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A New Airglow Photomeier* 
C. M. Purdy, L. R. Megill, and F. E. Roach 


(May 


The study of night airglow phenomena r 
the type depending upon the experimental 
analysis. It has recently become apparent 
simple, vet rugged and reliable, photometer 
intended to fill these needs. 


1. Introduction 


The photometer, which is shown in figure 1, is 
designed to monitor airglow intensities In one direc- 
tion only. The more likely choices for this direction 
are either toward the zenith or the pole. The present 
instrument is mounted on a base which gives zenith 
observations most readily, but for other directions 
a small modification in the base will give the desired 
angle. 

The photometer is designed to observe four colors 
in sequence. One and one-quarter minutes are 
required for each color observation. Of this time, 
25 sec are used in each of three operations as follows: 
(1) the instrument has an external shutter over the 
objective lens which gives a zero reading, 


(2) a 


FIGURE Photometer te lescope 


10, 1961) 


equires the use of various types of instruments, 
environment and the results desired for final 


to the authors that there is 


a need for a very 
The instrument 


which is described herein is 


standard light is moved in front of the objective lens 
giving a relative calibration of the instrument, (3) 
the objective lens is open to the sky giving a reading 
of the airglow intensity. A time sequence is shown 
in figure 2. The photomultiplier detector, the 
turret, the turret drive mechanism and the opties 
are contained in one unit—2 


1 in. high, 12 in. square. 
) 
o. 


Several views of the unit are shown in figure 


FIGURE sequence of operation of and lower 


u ppe / 


tele Scope turrets. 


2. Optics 


The optics used in this telescope consist of an 
objective and field lens as shown in figure 4. 

The optical design is the same as that used in an 
earlier model [St. Amand, 1953].! The field lens 
images the objective on the surface of the photo- 
cathode of the detector, which makes the instrument 
less sensitive to the position in the field of view of a 
point source such a star. This necessary 
because of the extremely nonuniform photocathode 
on the 1P21 photomultiplier used as a detector? 


as is 


3. Electronics 


The photocathode of the 1P21 has an S-4 spectral 
response, giving it some sensitivity as far as 6300 A. 
This enables one to observe all the prominent visible 

Figures in brackets 

lhe unit is readily 
cathode 

*Contribution from the 
Bureau of Standards, 
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airglow lines with only one detector. 
the tube is supplied by a commercial “ in high 
voltage supply, manufactured by the ¢ . Applegate 
Company. The addition of a 20K series resistance 
and a 10 mf 2,000 v d-c paper condenser in the output 
of the high volt power supply reduces the ripple at 
the output of the supply. The tube is operated at 
1,000 v 

The amplifying system consists of an Electro 
Instruments model A—12 d-c or A-14 d-c amplifier. 
This unit is an all-transistor amplifier that is chopper 
stabilized. Since, in general, the airglow intensities 
are independent of one another, four sensitivities 
will be needed for standard operation. The amplifier 
has a fixed gain and the individual gain adjustment 
for of the four airglow intensities is accom- 
plished by varying the size of resistor /?,. See figure 
5. The system is programed to change sensitivity 
each time a filter is moved into position in front of 
the telescope. 

The dark current of the photomultiplier presents 
a problem as it does in any d-c amplifier system. 
In order to eliminate its effect, a simple circuit Is 
set up at the input of the amplifier which may be 
adjusted to feed the amplifier an equal and opposite 
current as shown in figure 5. Since the dark current 
is a strong function of the temperature, this adjust- 
ment is made readily accessible and can be adjusted 
periodically by the observer. Each of the four 
sensitivity adjustments available for each filter has 
its own dark current adjustment. This allows the 
operator to set the dark current adjustment for each 
sensitivity. The voltage supply for dark current 
balance, shown in figure 5, a battery. Current 
drain is so extremely small that the shelf life of the 
battery is not shortened significantly. A galva- 
nometric strip chart recorder is used to record the 
data. In order to match the characteristic imped- 
ance of the recorder and accomplish optimum time 
response, a driver Was necessary as shown in 
figure 5. This for servo- 
tvpe recorders. 


each 


uses 


stage 


stage is not 


necessary 


4. Mechanical Drive 


The movements of both turrets are controlled by 
purely mechanical means. An exploded view of the 
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Voltage for 


is taken direct 


gear system is shown in figure 6. The rotating 
turret “b’’ which holds the filters and turret ‘a’ 
which holds the standard light, the sky opening, 


and a blank position for a zero, are driven by two 


concentric shafts. Each has its own gear system, 
but both are driven by the same motor. The motor 
runs continuously and a periodic motion of the 
turrets is obtained by a Geneva drive system. 


Descriptions of this type of oper: ation are available 
in the literature 's, 1930]. This type of motion 
combines high reliability with accurate positioning. 





5. Filter System 


The filters are mounted, as previously mentioned, 
in a disk which is periodically rotated into place. 
The particular line desired for study is selected by 
using an interference filter with a band pass of about 
50 A at half transmission. Colored glass filters may 
be used to eliminate side bands. Any wavelengths 
within the sensitivity range of the detector may, of 
course, be observed by changing filters. The 
present instrument is equipped with filters centered 
at 5300, 5577, 5893, and 6300 A. The last three are 
centered on airglow lines. The first is in a region 
where there is relatively little airglow and is used to 
help subtract out the integrated starlight and light 
from other extraterrestrial sources. 


6. Calibration and Standard Light 


The problem of calibration is one of the more 
difficult in the airglow field [Roach, 1958]. There 
is a well-calibrated photometer at the NBS Fritz 
Peak airglow station. The absolute calibration is 
obtained by comparison with this instrument. 

In order to maintain the calibration we are using, 
as previously mentioned, a standard light. 
phosphorescent source using a 
product as the activator. 

There are many phosphors on the market using 


This is a 
radioactive decay 


radioactive decay products for activation. An 
inexpensive source may be made from a radium 
activated phosphor which is applied as_ paint. 


These sources do exhibit some decay resulting from 
damage of the phosphor by the heavy ionization due 
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to the alpha particles from the radium. They are, 
however, quite cood for considerable periods and 
easily obtained. Most phosphors available have 
peak response in the green and little emission in the 
red. This is also true of some tritium 
phosphors which have been used. 

A phosphor which has been found to have good 
red emissions as well as the shorter wavelengths in 
the visible region is the U.S. Radium Corp. Isolite 
Standard Source. This contains a U.S. Radium 
Corp. red phosphor with carbon 14 as the activator. 
This phosphor is available with several types of 
radioactive activators. 


activated 


l 
mecnanisn 
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The authors recognize the contributions of Ralph 
C. Darr of the NBS Instrument Shop for excellent 
workmanship, particularly in the design of the 
details of the Geneva motion. 
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A Guide to the Use of the Modified Reflectometer 
Technique of VSWR Measurement * 
Wilbur J. Anson 


The theory of the modified refleectometer technique of measuring VSWR at microwave 
frequencies is discussed briefly, the operational procedure is outlined, and selected results are 


given of an unpublished error analysis. Much of the theory and procedure has been pub- 
lished in isolated papers. This paper unifies those details essential to the use of this tech- 
nique and includes procedural suggestions that have grown out of extensive experience with 
the technique. The error analysis provides the means to evaluate the accuracy of any 


particular measurement made with this system 


1. Introduction may be solved for the response in the form 


The modified reflectometer technique [Beatty and (1+KT,) 

Kerns, 1958; Beatty, 1959; Engen and Beatty, 1959] b;=b,k "taal ; 

recently developed at the Boulder Laboratory of the i: 

National Bureau of Standards is a method of making Ss. 

accurate VS R measurements in rectangular Wwave- k: Sar y Sur,’ 3) 

guide, This technique is capable of measuring low i song seak \ 

VSWR’s VSWR<2) to accuracies of 1 percent Ssil'e (I S33Dp) 

quite easily. Accuracies of 0.1 percent in VSWR are 

obtainable with extra care in the adjustments of the - 

modified reflectometer. kK S 
The purpose of this paper is to provide operating Sai 

knowledge of the modified reflectometer technique. | and 

Much of the material presented herein has been (i—Sul.¢) S Sal 

published elsewhere but this paper gathers in one SnT'< Ss, SpI'p 

place those details essential to the use ol this tech- S, (i-e7 (5) 

nique. The paper also offers suggestions on pro- rT’, ——— 

cedure which have grown out of extensive experience (I—Sule) Sisl'p 

with this technique, Also it provides selected Sailg (1—S33Ip)| 

results of an error analysis which are important in | |, 

evaluating the accuracy of any particular measure- Rearranged slightly, the equation becomes 


) (2) 


ment made with this svstem. 


2. An Outline of the Theory b, ii Sy Soo 


The modified reflectometer measures the magni- (I—Sule) Sislp I—PoP, 
tude of the reflection coefficient of a waveguide | Ssil’@ (1—S33I'p) 
termitiation at microwave frequencies. This method 
utilizes a directional coupler oriented to couple to the Now let S3,;=0 and T,,=0, then 
reflected wave with appropriately placed tuners 
adjusted so that the amplitude of the voltage wave h. , 
from the side arm is directly proportional to the (1— Sn P¢)A—S33P p) 
magnitude of the reflection coefficient. 

According to references 1 and 2, the seattering If b,, l;, and I» are constant, then 
equations for a three arm junction representing the 
combination of a directional coupler and two tuners 
(fig. 1) are 


b Sx SsoVr 


bs C Cecaatent yo (9) 


. . . ny “ , 1 a ) ase se SI j To 
a ee ee For a detector not phase sensitive 


h Soi + Soo t Ss nti b. CT, . (9) 
h Sod) + Salo + Sad (1) . . ‘ ie ‘ 
sie wiles Note that this last relationship is the desired response 
where a=), +-b,0¢, d2—boT,, and a,=b; Pp. These | Mentioned just before (1), i.e., the amplitude of the 
voltage wave from the sidearm is directly propor- 
*Contribution from the Radio Standards Laboratory, National Bureau of : } : . . ’ “ee 
Standards. Boulder, ¢ tional to the magnitude of the reflection coefficient. 
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Therefore, let [T,=T, (the unknown) so that 
lbs|,=|CT,|. Then let Tp=1, (a known quantity 
so that |6.!,=/CT,|. Now form the ratio 

bs, lg’ 10 
= —— { ) 
a | 


Thus, given the value of |T,!, and measuring the 
bs. 
ratio of the sidearm outputs h |? the value of 
) 
iT! may be calculated. 


3. Procedure 


Any practical use of the modified reflectometer 
technique depends upon an understanding of the 
physical conditions imposed by the theory. 

For |6;; to conform to (8) certain quantities as 
represented in (7) must remain constant. It re- 
quires a stable source (constant 6,) well isolated 
from load changes. It also requires apparent 
generator and detector impedances that are un- 
changing (constant T'g and ['p)—-requirements which 
can be essentially fulfilled by isolators placed at 
terminal planes 7; and 7; in figures 1 and 2. 

The remaining quantities of (8), namely 
Sy, Si, and S33, are constants of the 3-arm junction 
and therefore remain unchanging after adjustments 
are made for S;,;=0 and T.;=0. 


Soy 


Zl» 


3.1. Adjustment of S;,—0 and Ir;,=—0 


For a directional coupler with the arms labeled as 
in figure 2 the statement that S;;=0 is equivalent 
to saying that the coupler has infinite directivity 
since, in terms of scattering coefficients, directivity 


20 log 


Infinite directivity can be approached by ad- 
justing tuner A for no variation in the output ), 
when sliding a load of small reflection (/T,,|< 0.005) 
in the uniform waveguide. 

After making this adjustment, the condition 
r,,—0 can be approached by adjusting tuner B for no 
variation in the output (0; when sliding a load of 
large reflection (usually a shorting plunger) in the 
uniform waveguide. 

This adjustment sequence May need to be repeated 
several times since the two adjustments are not 


- - ~ 
NERATOR § [ NCTION ] TERMINATION 
T ‘ 
’ 7 3 
p* 03/bs en 
DETECTOR 
FIGURE 1. Three a unction representing the n odified 


eflectometer. 
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| 
I 


completely independent. However, if the reflection 
from the sliding load used in the directivity adjust- 
ment is small enough, then the two adjustments are 

more nearly independent and adjustment 
sequence will usually be enough. 


Nore: It is seldom necessary for the 
pletely The variation allowable 
slid in each adjustment depends upon the desired accuracy 
of the measurement to be made. This allowable variation is 
discussed in the error analysis and is displayed in figures 5, 
6, 7, and 8. 


one 


variation to be com- 


zero. while the loads are 


3.2. An Alternate Adjustment Procedure for S,,—0O 
and I.,—0 
Consider the definition for directivity in’ con- 
junction with (6) and observe that 6;;=0 wher 


r,=0 if the directivity is infinite. Therefore the 
adjustment could be made by attaching a flat load 
(T,=0) and adjusting tuner A for 0, =—0. But there 
is no commercial termination flat enough for the 
accuracy of adjustment often needed and therefore 
an adjustable sliding termination such as the one 
described by R. W. Beatty [1957] must be adjusted 
for T,;=0. 

The two adjustments for S,=—0 and T,=0, 
using an adjustable sliding termination, are carried 
out alternately. The adjustable termination is slid 
and adjusted so that the variations in |@;! are re- 
duced and then tuner A is adjusted so that the aver- 
age level of the variation lowered. When no 
variation in |b;| occurs as the adjustable termination 
is slid, and when |b, is zero, then S;,;=0 and T,=0. 
Note that after adjustment of tuner A the varia- 
tions will be greater than The greater 
variation results from increased — sensitivity 
infinite directivity is approached with the adjust- 
ments of tuner A. 

Adjustments for T,,=0 can be made after making 
the previous adjustment for infinite directivity 
The symbol [,; in (2), besides representing a group 
of svmbols as noted in (5), actually represents the 
reflection coefficient “looking in” arm 2 of the june- 
tion as indicated in figure 2. Therefore the condi- 


Is 


before. 


as 


tion T,,=0 is equivalent to saving that, at the 
terminal plane 7), this input reflection coefficient is 
Zero. 


> 


In figure 3 a scheme is illustrated for recognizing 
this matched condition. Power is coupled out of the 
main line to supply an auxiliary directional coupler 


cE a sone. | : 1] TERMINATION 


INIFORM WAVEG 
ECTION 


FIGURE 2. Diagran of modified reflectometer. 
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tuned for infinite directivity. This coupler is tem- 
porarily attached to the uniform waveguide with 
power shut off in the main line supplying the modified 
reflectometer. Tuner B is then adjusted for no 
output from the sidearm of the auxiliary coupler. 
The condition T,;=0 is thereby fulfilled. 

There must, of course, be enough isolation on the 
generator side of the modified reflectometer to main- 
tain essentially the same value of T,, during both 
this matching procedure and the normal operation 
of the reflectometer. Also, there must be enough 
isolation between the generator and whatever means 
is used for shutting off the main line power to insure 
no possibility of frequency pulling. 

The auxiliary directional coupler can be tuned for 
infinite directivity by utilizing the flat load previously 
This flat load is attached temporarily to 
the auxiliary coupler and tuner C adjusted for no 
output from the sidearm of the auxiliary coupler. 

This adjustment sequence may also need to be 
repeated several times since the two adjustments are 
not completely independent. However, if the adjust- 
ment for infinite directivity is carefully made, one 
adjustment sequence will be enough. 

Although this alternate procedure is sometimes 
easier to use, It does not provide sufficient informa- 
tion to evaluate the maximum errors due to imperfect 
adjustments of the tuners. Therefore, the adequacy 
of the adjustments must be checked by sliding the 
small and large reflection loads in the uniform wave- 
guide and by observing whether the variations in 
6; are within the allowable limits for the measure- 
ment accuracy needed. 

If the variations are not within the allowable 
limits it is important that the final adjustments be 
made with the small and large reflection sliding loads 
in the uniform waveguide. Otherwise, reflections from 
the waveguide joints will detune the infinite direc- 
tivity adjustment. 


discussed. 


3.3. Another Alternate Adjustment Procedure for 


Sar 


By adding an auxiliary channel as illustrated in 
figure 4, the adjustment for infinite directivity 
(S3;=0) made considerably easier and faster 
(Beatty, Engen, and Anson, 1960]. This auniliary 


Is 
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FIGURE 4. 


Auxiliary channel technique 
channel permits introduction of a signal of the proper 
phase and magnitude to null the signal component 
indicative of finite directivity. 

From (2) the response from the sidearm is 


, |1+ATP, 
b, b,K i r. r. Q 
If |T.;| is small, i.e., {T;{ <0.01, then a very good 
approximation is 
ef 
b bkK >+T;]- (11) 
K 


The following description of the procedure is best 
understood by referring to (11) and the vector 
diagrams in figure 4. 

1. With the shorting switch in the auxiliary arm 
closed, arbitrarily position the small reflection load 
in the uniform waveguide and adjust tuner A for 
b,—0. This corresponds to step (a) in figure 4 
wherein 1/A is oppositely directed and equal in 
magnitude to |T;). 

2. Slide the load for maximum |),!. This corre- 
sponds to step (b) in figure 4 wherein T, and 1/K 
are in phase. 

3. Open the waveguide switch in auxiliary channel 
A and adjust the phase shifter and attenuator in 
that channel for }6; =0 again. This corresponds to 
step (c) in figure 4 which illustrates a signal A having 


been introduced to null r,/. 


l 
| 
4. Introduce 6.02 db of attenuation into the aux 
iliarvy arm to reduce the amplitude of signal A by 
one-half as illustrated in step (d) in figure 4. This 
A/2 cancels the 1/A term in (11) and results in 
infinite directivity. Reeall that (T, assumed 

small so that (11) is a good approximation. 

The directivity obtained by this procedure can be 
checked by sliding the low reflection load and ob- 
serving the variations in |b;| as noted in section 3.1. 
If the directivity obtained thusly is not high enough, 
then trial and error adjustments using the phase 


IS 
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shifter and the attenuator in the auxiliary arm or 
tuner A will allow closer adjustment. However, 
adjustments using the components of the auxiliary 
arm are in general easier because these components 
usually present better resolution than tuner A. 

To avoid possible interaction between the two 
channels, two isolators should be included as shown 
in figure 4. Other tuning [Engen, 1960] could be em 
ployed to reduce interaction but these two isolators 
should suffice. 

3.4. Measurement of |T 


u 


Certain details vary from system to system and 
from operator to operator depending upon personal 
preferences, equipment available and the accuracies 
needed. Therefore, no attempt has been made to 
write this guide for a specific arrangement of equip- 
ment. 

However, regardless of the arrangement of the ap- 


paratus in the system, the ratio can be conven 


iently measured by some adaptation of the following 
sequence: 

A. Adjust S;;=0 and T,;=0 by some combination 
of methods as previously outlined. 

B. Check the adequacy of the adjustments to in- 
sure the achievement of the measurement accuracy 
required. , 

C. Attach the microwave termination of unknown 
reflection coefficient T,, to the uniform waveguide. 

D. Set \6,; at some convenient reference level on 
some suitable detector and note the setting of a cali 
brated attenuator. This attenuator, of course, must 
be properly matched [Beatty, 1954; Schafer and 
Rumfelt, 1959] to measure accurately the voltage 
ratios at the detector. 

E. Replace the unknown with some standard ter- 
mination whose reflection coefficient is accurately 
known. 

F. Return |), to the reference level by adjustment 
of the calibrated attenuator mentioned in step (D). 

G. Record the difference in settings of the eali 
brated attenuator between step (D) and step (F). 
Calculate the value of |T,,) from this difference, i.e., 


R=20 logio where F? is the difference in decibels. 


This calculation can be aided by reference to a table 
[Beatty and Anson, 1960] of reflection coefficient 
versus return loss. 


‘4: The Error in the Measurement of |I',,! and 
VSWR Due to S3, 40 and l;4+0 


As noted previously, the variation in 6;) allowable 
while the loads are slid in each adjustment depends 
upon the desired accuracy of the measurement made. 
This allowable variation can be computed from the 
following equations or taken from the graphs in- 
cluded for certain measurement situations. It i 
recommended that the adjustments be checked by 
this method both before and after measurements are 
made as an indication of changes created by accl- 
dental bumps or frequency drift. 
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4.1. The Error Due to S;,#0 


it can be shown that the maximum fractional error 
(assuming the worst phase combination) in the meas- 
urement of |T,|, where S;,40 but all other conditions 
are met, is given by the expression 


where be =reflection coefficient of the standard ter- 


mination, 


l,=reflection coefficient of the unknown ter- 
mination, 
and ‘ 101/294 -] 
K = , 
(10%1/2—])|P, 


in which P, is the db variation in |b, as the small re- 
flection load T,, is slid in the umform waveguide. 

The corresponding maximum error in the measure- 
ment of VSWR, when S,, +0, is obtained from (12) 
applying the following relation: 


de, 2\T dvr 


TI IP VSWR. 


where o,, 


4.2. The Errer Due to r.;40O 


The maximum fractional error in the measurement 
of T,!, when T,,—0 but all other conditions are met, 
is given by the expression 


d\T,|_— (Pal +IP.|) [Pe ais 
r. i—T,f 


where 


l,=reflection coefficient of the standard termination, 


lr, =reflection coefficient of the unknown termination, 
and I, is obtained from 
1+/T,,° 
R.=20 logiy9 — = 2 
re 


in which F, is the db variation in | 
large reflection load 
whose |T',|™1). 

The corresponding maximum error in the measure- 
ment of VSWR, when I,,;#0, is obtained from (14) 
by applying (13). 


observed as a 
(usually a shorting plunger 


4.3. The Error Due to Both S,,#0O ard I.,;+0 


Seldom, if ever, will Sj, or T,,; be adjusted per- 
fectly. However, it can be shown that the maximum 
error due to both Sj, #0 and T,,;40 simultaneously 
the sum of the individual maximum. errors. 
Therefore, these individual errors may be obtained 
from (12) through (14), or from figures 5 through 8, 
and added to form the maximum error of the meas- 
urement due to S.,+9 and P,;, 40. 
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Fr URE 5 The marimum erro 
modified reflectometer 


r 0 


in measuring VSWR by the 


technique due to S3;#0 and assuming 


Example: Suppose there is 0.5 db variation in 
bs! as a small reflection load \T, 0.005 is slid in the 
uniform waveguide and there is a 0.02 db variation 
as a short (\P 1) is slid. What is the total limit 
of error in the measurement of the VSWR of a 
termination whose reflection coefficient is |T,,! ~0.2? 
From figure 5 the maximum error in measuring 
VSWR due to S340 is 0.036 percent; from figure 6 
the maximum error due to T,,;#0 is 0.058 percent. 
So the maximum error (due to both S:,4#0 and 
r.,;~0) in measuring the VSWR of the termination 
is 0.094 percent. 


5. Procedural Suggestions 


1. Adjustment of the tuners for no variation in 
the sidearm response |b3| is easier if some systematic 
approach is used. One such approach is the follow- 
ing: 

A. Slide the load for either a maximum or mini- 
mum response. In this example assume that it was 
slid for maximum. 

B. Adjust a stub of the tuner in the direction 
which lowers the detector response. (Although this 
example uses a multiple-stub tuner the same _ pro- 
cedure is usually adaptable to other types of tuners.) 

C. Slide the load for a minimum. 

LD. Adjust the same stub in the same direction as 
before. If the increases, 1.e€., moves 
toward the average of the variation, then further 
adjustment of that tuning stub using steps (A) 
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FicureE 6. The maximum error in measuring VSWR by the 
modified reflectometer technique due to T2;#0 and assuming 
S3:=0 

through (D) repeatedly will systematically reduce 

the variation. If the response decreases, i.e., moves 

away from the average of the variation, adjustment 
of that stub will probably be useless and another 
stub should be tried. 

The adjustments for no variation as the load is slid 
are usually quicker if several effective stubs are 
located by this procedure and then alternately 
adjusted. 

If no effective stubs can be located, make some 
random adjustment of the tuner to raise or lower the 
detector response and repeat the previous procedure. 
This is a trial and error process. Experience with 
your own equipment will provide “feel’’ for the 
detector response level resulting when the sliding 
load variations have been reduced. 

2. An approximation to the detector response 
level resulting when the variation in the sidearm 
output has been reduced for the infinite directivity 
can be obtained as follows: 


A. Slide the low reflection termination for a 
minimum. 

B. Adjust tuner A for a null. 

C. Slide the low reflection termimation for a 


maximum. About 6 db down from the maximum 
is the detector response level resulting when infinite 
directivity is achieved. This level can serve as a 
guide when making the adjustments for no variation 
as the load is slid. 
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= +t eS Hay ea reegs This procedure avoids the one possible false 
citi * fo | adjustment which could result in an approximately 
POTS * AND Ik ; | ; eee is ; 
SO SK constant output with Kk? r,. This misleading 
‘ Miccc adjustment is seldom encountered except when 
1 AXON Ht NOTE THAT T,, IS ASSUMED 4 adjusting the equipment the first time. It some- 
BS KI NI] SS Oe times occurs at that time because the system lacks 
N NWN Ih ae EE sufficient power sensitivity for detection of the 
. NOS ag ee a 
mm i a8 ae ee ee A SS approximately constant output signifving K 7 
: Boe aw SSR iC 
3 IAG +4 WSO \ Th which occurs at a low power level. 
r ~% +4 =e SE ia i ‘ ; ° . 
3 Of IN | I \ ». It is recommended that the adjustments for 
a , % } : ty : r SK Nyse S;, ond T., be checked both before and aftera VSWR 
= +4 WO or reflection coefficient measurement is made to 
z Ni N\ Na N detect possible changes created by accidental bumps 
5 ok Be OS! Oh | | N NN or frequency drift a . . 
2 TY ‘ oo. SY 1. The reflection coefficient magnitude [,) of the 
SEN KC NST = NSN small reflection load need not be determined to any 
RAPE No tN TOT great accuracy. Computation of the maximum 
witH i <0025 Not TY Serre ri error requires only that the actual value of |T, is 
ata TNT ST RIT less than the value used in the computation. 
Nu : KING This reflection coefficient magnitude can be de- 
\N . , \ termined accurately enough as follows: 
NXT NUUTN A. Make the normal adjustments of the system, 
. repeating the adjustment sequence until adjustment 
, for T,;=0 does not require further adjustment of 
Figure 7. The maximum error in measuring \T,! by the S2,=0. 
modified reflectometer technique due to S3,4#0 and assuming B. Slide the small reflection load and observe the 
Pas=0. maximum and minimum (|b;. There should be 
te ~— little, if any, detectable difference between the 
minimum and the maximum. Therefore, it is usually 
sufficient to position the sliding load for the average 
. response and measure [,, for this position using 
the magnified reflectometer technique. For a truly 
conservative figure, however, measure the reflection 
coefficient magnitude of the load in the position of 
minimum |b 

If the standard termination being used is a “short” 

qr 1) the ratio R=20 log, . mav be too large 
for the calibrated attenuator of the system. Then 
it may be necessary to use a standard of smaller 
reflection coefficient in order to measure ([,,!. 

5. The adjustment of the tuners as outlined in 
ae 5.1 is an efficient technique if effective stubs can be 
‘2 i | found. The following technique aids in finding 

| | effective stubs for the infinite directivity adjustment, 
A. Adjust tuner A and the sliding load for \bs|=0. 
t+ B. Slide the load for maximum |} 
++ C. Now at least one of the stubs is an “effective” 
stub. Continue the adjustment as outlined in 5.1. 
05 6. When using the small reflection load (\T,)- 
0.005) the signal level yay disappear into the noise 
as the adjustments for infinite directivity are at- 
tempted. This difficulty indicates that the detection 
system does not have sufficient power sensitivity. 
Kather than invest in more sensitive equipment, a 
load withggreater reflection (perhaps || ~0.025) 














= - can be used for the infinite directivity adjustment 

provided the¥adjustment sequence outlined in sec 

db VARIATION AS SHORTING PLUNGER IS SLID . . , ae i : 

; tion 3.1 is’ repeated as indicated. However, this 

hiGURE &, The maximum error in measuring Ty hy the l ; : : . é : ; 4 

modified reflectometer technique due to T2;#0 and assuming a ternate procedure is) consideradiy more time 
S3,=0. consuming, 
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The author is indebted to Edward Niesen for his 
invaluable assistance im the operation and main- 
tenance of the equipment, to R. W. Beatty for 
his encouragement and helpful suggestions during 
the writing of this guide, and to David Wait and 
Michael Brady for their work on the error analysis. 
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An X-Ray Ditfractometer Cryostat Providing Temperature 
Control in the Range 4 to 300 ‘K’" 
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The diffractometer cryostat described earlier has been modified to provide a 


for controlling temperature ¢ 
Materials with vapor 
of insulating vacuum, 
range 4 to 20 °K 
thermal expansion of neon are 


pressures as great 


given 


1. Introduction 


In 1958 the 
in designing a ervostat 


authors collaborated with others 
for obtaining X-ray diffrac- 
polverystalline solids formed by 
deposition on a surface cooled by liquid helium [1]. 
It was used with success in studying 
temperature solids condensed from the 
phase, but the design imposed serious limitations 
on the choice of experimental conditions. The 
copper blow k on which samples were deposited Was 
soldered directly to the bottom of the helium 
reservoir, and gaseous samples were introduced 
into the insulating Vacuum through a tube 
directed at this block. Most of the material froze 
where it first impinged on the helium-cooled surface, 
and could be examined by X-ray diffraction methods 
through thin beryllium windows in the outer wall. 
Deposition and examination at the boiling tem- 
peratures of conventional crvogenic liquids could 
be accomplished readily, but, at intermediate 
temperatures, observations could be made only 
while the block was warming after depletion of 
the refrigerant. The normal rate of warming (1 
per min) was too high to permit adequate X-ray 
data to be obtained. It was not feasible to maintain 
au desired a by heating the block elec- 
trically while it was in contact with liquid helium, 
and when ¢ joke techniques for isolating and heating 
it were tried, the increase in vapor pressure of the 
sample usually interfered. When the vapor pressure 


tion patterns of 


low- 


FrxSeCOUS 


sOTLEe 


space 


exceeded 107° mm of Hg, diffusion of the sample 
into the insulating vacuum = spaces resulted in 
excessive heat transfer, while condensation on 
surfaces colder than the block eventually led to 


disappearance of the sample from the surface in 
the X-ray beam. 

These difficulties have been overcome in recent 
modifications which permit the copper block to 
be isolated from the helium reservoir. It can be 
heated electrically or cooled at any time to the 
med under the 


supported by 
ures in brackets indicat 


National Bureau of Standards Free 
the Department of Defense. 
e the literature references at the end of this paper. 


to demonstrate the 
maintaining a desired temperature while working with a volatile material. 


», 1961 


means 


and for confining the vapor in Srenretns with volatile samples. 
as Ll atmosphere 
and temperatures can be controlle dwith an 
decreasing to LOOK 4 


can be handled without loss 
accuracy of +0.1 °K in the 
Data on the lattice constant and 
suitability of the apparatus for 


“atm 


300 °K. 


temperature of the refrigerant. An electronic: con- 
troller maintains the desired temperature, and a 
Mylar window separates the specimen chamber 


from. the insulating vacuum spaces so that sample 
vapor pressures as high as 1 atm can be tolerated. 
Moreover, heat leaks have been reduced so that 2 
liters of liquid helium last for as long as 6 hr. 

The method of isolating the copper block from 
the helium reservoir is an adaptation of that used 
by Henshaw [2] in his neutron diffractometer 
cryostat. In most other designs, a low-con ductivity 
link is introduced between the sample and the 
liquid helium reservoir before an experiment is begun. 


Thereafter, the sample cannot be cooled to the 
temperature of the refrigerant and may require 
hours to reach its equilibrium temperature. Hen- 
shaw overcame these difficulties by placing’ his 


specimen 
stainless 
helium 


block at the lower end of a thin-walled 
steel tube joined to the bottom of the 
reservoir. Initially, liquid helium fills this 
tube and is in direct contact with the specimen 
block. When the block is to be heated, a plug is 
inserted in the top of the tube to exclude the liquid. 
Because of the low conductivity of the stainless 
steel used for the tube, the block can be heated to 
10 or 50 °K before the increase in the evaporation 
rate of the liquid helium becomes excessive. -igure 
1 shows how this design has been adapted to the 
diffractometer cryostat. 


2. Modification of the Cryostat 


By comparing figure 1 with the corresponding 
figure in the earlier paper [1] it can be seen that many 
features of the cryostat remain unchanged. It is a 
conventional double Dewar consisting of a 
liquid helium vessel (21) surrounded by a liquid 
nitrogen container (22) od radiation shields (19, 3) 
which screen the helium vessel from thermal 
radiation. 

In modifying the cryostat the original copper 
block was detached from the helium reservoir at 
the soft solder joint (7). The inner vacuum wall (10), 
which had separated the vacuum around the helium 


vessel 
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Figure 1. Vertical section through the X-ray diffractometer 
cryostat. 
1) Thermocouple, (2) heater bobbin, (3) radiation shield, (4) copper block, 
5) thin-walled support tube (stainless steel), (6 pper rin 7) soft-solder 
joint, (8) ring with centering pins, (9) flange with “‘o’’-ring seal, (10) inner 
vacuum wall, (11) heater leads, (12) vacuum valve, (13) radiation baffle 
14) helium-vessel neck stopper, (15) block vent tube, (16) liquid-helium depth 
probe, (17) helium vent, (18) nitrogen fill-tube and vent, (19) radiation shield, 
20) spacing spring, (21) helium vessel, (22) nitrogen vessel, (23) inlet-tube 
heater terminal, (24) bypass valve, (25) sample-inlet valve, (26) manifold block, 
(27) sample-inlet tube (stainless steel), (28) Mylar plastic window, (29) sampk 


inlet-tube heater, (30) sample-deposition surface, (31) aluminum foil window 


0.00035-in. thick), (32) opening for 0.025-in. thick bervllium X-ray window, 
33) translation-adjustment screw 34) rotational-alinement gear sector, 
35) plug, (36) internal siphor 


‘vessel from that in the specimen chamber, was then 
cut away from the helium vessel at the lower end to 
reduce heat conduction to the liquid helium. <A 
ring (8) with three centering pins was fitted to the 
lower end of the inner vessel. The pins contact the 


| 
| 
1 


| remaining portion of the inner wall and prevent 
lateral motion. They do, however, allow for a 
vertical component of relative motion so that the 
two vessels nay contract independently when filled 
with cold liquids. 

Replacing the original copper block is an assembly 
which makes it possible to introduce the thin-walled 
tube (5) between the bottom of the helium reservoir 
and the copper block (4) without changing the loca- 
tion of the deposition surface (30) relative to the 
windows and X-ray beam. The copper ring (6) is 
fitted with a %-in. o.d. tube that extends for 1%4-in. 
up inside the helium reservoir. At the top of this 
tube is a second ring which becomes, in effect, the 
new bottom of the liquid helium vessel. The %-in. 
o.d. tubing (5) of 304 stainless steel that supports 
the copper block is 2!s-in. long. The wall thickness, 
originally 0.035 in., has been reduced to 0.010 in, 
over a length of 2 in. to reduce the cross sectional 
area Of the heat path. Because of the low thermal 
conductivity of 304 stainless steel (~2* 107% watt 
em! °K~'), the heat flow along this path corre- 
sponding to a block temperature of 40 °K is only 
about 0.025 watt 

The plug (35) which is inserted through the neck 
of the crvostat is used to exclude liquid helium from 
the region of the block when it is to be heated. A 
diagonal passage through this plug connects the 
internal siphon (36) with the tube extending to the 


bottom of the Space below and permits liquid 
trapped in this region to return to the helium 
reservoir, Another passage connects this region 
with the vent tube (15). These passages may not 


be visible when the figure is reproduced. 

A new sample deposition block (4) was cut from 
a eylinder (1's-in. diam by %-in. high) of electrolytic 
tough pitch copper, selected for its high conductivits 


at low temperatures. A ‘s-in. wide cut was made 
halfway through the cylinder, leaving a_ plane 
| surface (30) on which samples are deposited. The 


block Was gold-plated to provide a chemically inert 
surface that elves an X-ray pattern suitable for use 
as an internal standard 

The uncut portions of the cylinder above and 
below the deposition surface provide circular flanges 
around which the Mylar window (28) is cemented 
using an epoxy resin [3]. The window is 0.0007 in. 
thick and attenuates the copper K, radiation 
ordinarily used by only 8.5 percent. [Tt 1s vacuum- 
tight and readily supports a pressure differential of 
latm. The one failure of this window occurred ata 
pressure difference of about 3 atm. 

A 23-ohm Manganin resistor wound on a copper 
bobbin (2) is provided for heating the block. The 
carbon (100 ohm, '%>o watt) used as the 
sensing element for the temperature controller, is 
mounted in a 'Me-in. hole drilled along the axis of 
the bobbin. The whole assembly is pressed into a 
\4-in. hole in the copper block. 

An inlet tube (27) connects the region enclosed by 
the Mylar window with the manifold (26) on the 
outside of the cryostat. Stainless steel tubing }s in. 
in diameter with a 0.005-in. wall thickness is used. 


resistor, 
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To keep samples from freezing in this tube, it js 
lined with Teflon tubing through which is passed an 
elongated helix of #30 Nichrome wire (29). The 
inner end is grounded to the copper block and the 
outer is connected to a terminal (23) in the manifold 
so that a current may be passed through the wire if 
blockage occurs. 


The bypass valve (24) may be opened while pump- 


ing to equalize pressures on the two sides of the 
Mylar window, and closed before the inlet valve 


(25) is opened to admit the sample gas. 
cryostat turns through 90° 
is scanned, 


Since the 
while the X-ray pattern 
flexible }-in. annealed copper tubing is 
used to connect this valve to the stationary specimen 
preparation train, flow 
auxilliary vacuum pump. 

Other features that 
the text are identified 


Fages, 


manometer, and 
have not been described in 
in the legend of figure 1. 


3. Temperature Controller 


Numerous controllers have been described in the 
literature [4], but the one used to control the tempera- 
ture of the specimen block differs from all of these 
and will be described briefly. It uses the 100-ohm 
carbon resistor in the heater bobbin as the tempera- 
ture sensing element (fig. 2). The resistor forms one 
arm of a Wheatstone bridge which is set to balance 
at a resistance corresponding to the desired tem- 
perature. The servo amplifier from a recording 
potentiometer amplifies the output of the bridge and 
drives a two-phase reversible motor to set a variable 
autotransformer. heater 


In this way, power input 
is automatically increased or decreased until the 
desired temperature is reached. To uuprove the 


stability of the System, a rate signal is obtained from 
a small d-c generator on the shaft of the balancing 
motor. While the system is approaching the con- 
trol point, this generator produces a voltage propor- 
tional to the velocity. A portion of this voltage is 
fed back to the amplifier input, phased so that, when 
amplified, it tends . drive the motor in the reverse 
direction. ‘Thus, a braking torque is provided that 
is very effective in reducing the amount of over- 
shooting at the control point 


FIGURI ae Temperature controller circuit. 


In order to reduce the time requiced for the resistor 


to sense a change of power input, the compact 
\vo-watt size Is used. It Is cemented into a '{¢6-in. 
hole in the heater bobbin. One lead is soldered 


directly to the block, and the other is cut very short 
to decrease its heat capacity. 

A separate resistor (56 ohm, 
in the block and used for measuring temperatures 
in the range below 30 °K. 


In addition, a thermo- 
couple (Au + 2.11 atomic % Co versus Cu) is used 
throughout the range 4.2 to 300 °K. 

The use of resistors in low temperature 
thermometry has been described by Clement and 
Quinell [5]. Two characteristics of these thermom- 
eters are their extreme sensitivity at low tempera- 
tures and the rapid decrease of sensitivity with 
increasing temperature. Thus, by measuring resist- 
ance to 1 part in 5,000, it js possible to detect tem- 
perature changes of 0.0001° at 2 °K, 0.01 °K at 
20 °K, or 1.0 °K at 300 °K. In the present applica- 
tion, the inherent sensitivity cannot be fully utilized, 
and below 20 °K the precision of control is limited to 
+0.1 °K. Above 20 °K the precision gradually de- 
creases to about +1°K at room temperature. 


1 watt) is mounted 


carbon 


4. Operation 


In order to fill the cryostat with liquid helium, the 
neck-stopper (14, fig. 1.) and the entire plug assem- 
bly are removed so that the transfer siphon can be 
inserted. When the inner vessel has been filled to 
approximately the level shown, the siphon is removed 
and the neck-stopper is replaced after drawing the 
plug (35) all the way up so that it remains in the 
neck tube to be cooled by the ese aping vapor. The 
sample is usually introduced at this time by opening 
the inlet valve (25) and allowing the gas to flow in at 
a rate of 20 to 30 em’/min (STP) until the thickness 
of the deposit is sufficient to reduce the intensity of 
the X-ray pattern from the gold substrate by a 
predetermined amount. When the block is to be 
heated, the plug is lowered and seated in the top of 
the tube (5) to exclude liquid helium. That trapped 
below the plug is returned to the reservoir through 
the internal siphon (36) by closing the vent (15) so 
that pressure develops in the well. 

With helium out of the well, the block may be 
heated electrically, using the controller to maintain 
a desired temperature. The setting of the bridge is 
determined by refe rring to curves based on the semi- 
empirical equation of Clement and Quinell for the 
calibration of carbon resistance thermometers. Ex- 
perience indicates the best settings for the 
sensitivity and damping controls. 

The relationship between power input and block 
temperature for three refrigerants is shown in figure 
3. Some heat reaches the block by conduction along 
the electrical leads and sample inlet tube even when 
the electrical power input is reduced to zero. As a 
result, the temperature of the block rises 10 to 20° 
above the temperature of the refrigerant. The in- 
ternal siphon (36, fig. 1.), which was not a part of 
the system described by Henshaw, was introduced 


soon 
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FIGURE 3. Variation of block temperature with heater power 
wmnpul for (a liq wid helium b) liquid h sdroqe n, and (« liq uid 
nitrogen. 
For the st refrigerant, the effect of re jual gas in the region below the plug 
een by paring cur btained with nitrogen gas at tn t we 
WwW u € , obt | ng the region with styrofoan { 
obtained by evacuating the region to eliminate both conduction and conve ) 


in an attempt to control temperatures in the interval 
between the boiling point of the refrigerant and the 
zero-power intercept. By pumping on the block 
vent (15) through a needle valve, helium (for ex- 
ample) can be drawn in at a controlled rate. If 
desired, the automatic temperature controller may 
be used simultaneously. This method is effective 
down to 10 to 12°K. At lower temperatures two- 
phase flow evidently commences and results in dis- 
continuous cooling. 

The arrangement that bas given the best results 
in the range 4.2 to 15 °K is a sliding-piston plug 
which is substituted for the plug shown (35, fig. 1.). 
[It was hoped that any desired temperature could be 
reached by sliding the plug up or down to vary the 
effective length of the heat path along the thin- 
walled tube (5). However, with the block vent (15) 
closed to create a positive pressure, keeping helium 
from leaking around the piston, the lowest tempera- 
ture reached with the plug all the way down is 12 
°K. To reach lower temperatures it is necessary to 
provide an escape to the atmosphere for helium gas. 
The rate of leakage around the plug is determined 
by the rate at which gas is allowed to escape and 
this is easily controlled by a needle valve (not shown). 
Each full turn changes the block temperature by 
very roughly 0.5 °K. The temperature remains 
constant (+0.1 °K) for periods of 10 min or more, 
and automatic control is available if desired. The 
valve is simply opened enough to give the lowest 
temperature required and, thereafter, the controller 
is used to adjust the heater input. 


5. Application 


Temperature was controlled in this way to obtain 
data on the lattice constant of neon at various tem- 
peratures in the range 4.2 to 22 °K (fig. 4) [6]. 
Values of the lattice constant were determined by 
the X-ray method except for that at the triple point, 
which was calculated from the density data of 
Clusius [7]. Both carbon resistors and the thermo- 
couple were calibrated in place at the boiling point 
of helium, the triple point of hydrogen and the 
boiling point of hydrogen. In addition, readings 
from both resistors and the thermocouple were com- 
pared with the vapor pressure of hydrogen through- 
out the temperature range 10 to 20 °K. Tempera- 
ture Measurements believed to be accurate to 

0.1 °K and lattice constants to +0.001 A. The 
neon data demonstrate the importance of accurate 
control at temperatures below 20 °K. There are, 
of course, many other problems requiring tempera- 
ture control that can be studied with this apparatus. 
For example, solid oxygen is known to exist in three 
forms, but, as far as we know, the structures of two 
of these and the effects of disorders that commonly 
occur, are still subjects for further study. It is hoped 
that this cryostat will contribute to the solution of 
these and other problems in the structure of low 
temperature solids, 
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FIGURE 4. The lattice constant of neon as a function of tem- 
perature, showing the reproducibility between three ex peri- 
me nts 

he value at 4.2 °K is the average of four measurements. That at the triple 

point is calculated from the density reported by Clusius. 


6. References [5] J. R. Clement and E. H. Quinell, The low temperature 
characteristics of carbon-composition thermometers, Rev. 
, ‘ Sci. Instr. 23, 213 (1952). 

fi] I. A. Black, L. H. Bolz, F. P. Brooks, F. A. Mauer, and [6] L. H. Bolz and F A. Mauer, Lattice constants and thermal 
H. S. Peiser, A liquid-helium cold cell for use with an X-ray expansion of solid Ne, A, Kr, and Xe, to be submitted to 
diffractometer, J. Research NBS 61, 367 (1958) RP2907. J. Research NBS. 

[2] D. G. Henshaw, Atomic distribution in liquid and solid [7] K. Clusius, Atomwirmen und schmelzwirmen von Neon, 
neon and solid argon by neutron diffraction, Phys. Rev. Argon, und Krypton, Z. phys. Chem. B31, 459(1936) 
111, 1471 (1958). 

[3] R. 5. Hickman, R. W. Kenny, R. C. Mathewson, and 
R. A. Perkins, Liquid hydrogen targets of adhesive- 
bonded Mylar plastic, Rev. Sci. Instr. 30, 983 (1959). 

[4] R. D. Goodwin and J. R. Purcell, Direct-coupled power 
amplifier for cryostat heating control, Rev. Sci. Instr. 28, 

581 (1957). (Twenty-six references are cited.) (Paper 65C4—75) 


229 








URNAL OF RESEAR( Bure 


~H of the National 
|. 65C, Ne 


4 Oct 


ot Research—Enginee 
196] 


ring and Instrumentati 


Jecember 


Apparatus for Determination of Pressure-Density-Temper- 
ature Relations and Specific Heats of Hydrogen to 
350 Atmospheres at Temperatures Above 14 °K’* 


Robert D. Goodwin 


June 


Method and apparatus ure designed for 


spaced, PVT 

pressure-temperature observations at nearly 
instalt method Temperatures of the 
electric heating under control of the 
ment 


fluid Calorimetric experimentation is 
calorimetric heat supply 
Details are given of the 


inde pendent data 


PVT 


calibrations, 


1. Introduction 


The high-density physical properties of hydrogen 
needed for technological applications [5, 10). 
Thermodynamic functions hes be computed from 
certain thermal data combined with a wide range of 
precise mechanical properties (PVT) and vapor- 
pressure data [1, 2,3, 9]. The purpose of this paper 
the essential features of apparatus 
used to obtain PVT data for parahvdrogen from 16 to 
100 °K and from 2 to 350 atm, and to obtain specific 
heats of compressed liquid and fluid in the same 
range. Advent high-speed digital computers 
renders practical the handling of a large quantity of 
data, produced by closely -spaced observations with 
accelerated methods to be described. 


ure 


is to describe 


of 


2. Experimental PVT Method 


The PVT method and its numerous calibrations 
emphasized in preference the better-known 
calorimetric techniques, outlined below. The popu- 


are to 


lar Burnett method for PVT determinations [17] is 
not suitable for compressed liquids [9]. That of 
Holborn and coworkers |20, 21, 22], termed the 


Reichsanstalt method [2], involves essentially direct 
measurement of variable. As employed by 
Michels and coworkers, a sample of fluid of known 
PVT behavior at normal temperature is confined in 
the piezometer and pressure system at that 
temperature. Pressure of the sample then decreases 
rapidly with tempe rature of the pesomoter [25]. 
In the Ohio State University modification for gaseous 
states, adjustments for obnoxious volumes of capil- 
lary and gage were eliminated by placing a valve in 
the ervostat. This necessitated an independent de- 
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20, 


more 
and specific heat data than realized by previous procedures. 
constant density 
essentially 
measuring 
and control are integrated with a high-pressure calorimeter for compressed liquid and 
accelerated by 
and of a d-c power regulator developed for automatic shiel 
adjustment computations, 
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1961 


rapid determination of accurate, closely- 
A sequence of 
a modified Reict 
are regulated by 
Instruments for measure- 


is made by 
adiabatic piezometer 
thermometer. 


is- 


of for the 
1 control. 


and comparisons W ith 


use an electronic battery 


termination of the amount of sample following every 
P-T determination [24]. For liquid states they em- 
ploved a boiling hvdrogen thermostat. 
small portions of the sample were released from the 
piezometer or ‘‘pipet’”’ to obtain a series of isother- 
mal P-v determinations from a given filling, each 
P-y point again requiring an independent volumet- 
ric determination [23]. 

Average time per determination 1s greatly reduced 
in the present method. Temperature of the nearly 
adiabatic pipet is varied to obtain a series of P-T 
determinations at nearly constant density. Tem- 


Successive 


peratures are automatically controlled at exact, in- 
tegral values, to permit handling the data as 
isotherms. Only one determination of the amount 


of sample is required following each experimental 

‘“pseudo-isochore’’. The portion of the total sample 
which was in the pipet during each P-T determination 
is calculated by subtracting the computed amounts 
residing simultaneously in capillary and diaphragm 
cell under those conditions. For this small adjust- 
ment, the known and estimated PVT behavior of 
normal hydrogen is employed. 

With reference to figure 1, a sample of the experi- 
mental fluid is confined by valve C—4 to the following 
system: the heavy-walled copper pipet, the stainless 
“transition” capillary tube in the cryostat, the capil- 
lary tubes and valves at room temperature, and the 
pressure-sensitive diaphragm cell (N. P. D.). Vol- 
umes of these elements are calibrated independently. 
The latter transmits pressure of the fluid to the oil of 
a piston-type, deadweight gage. With reference to 
figure 2, liquid hydrogen refrigerant resides in the 
tank. The pipet is cooled strongly as required by 
liquid hydrogen reflux action from hydrogen gas jn- 
troduced to the thin-walled, stainless steel ‘“reflux’’ 
tube supporting the pipet. for automatic tempera- 
ture regulation by electric heating, under control of 
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the platinum resistance measuring thermometer, 
pressure in the reflux tube is reduced to provide con- 
trolled cooling by CAS convection, alded by a copper 
rod suspended in this tube. 

An experimental run consists of measurement of a 
sequence of pressure versus temperature points, be- 
ginning at the lowest temperature. The total quan- 
tity of confined fluid then is determined by releasing 
it as gas into volumetric system D of figure 1 and 
measuring P, V, and T at about normal conditions, 
accounting also for the gas remaining in all obnoxious 
volumes and in the pipet at this pressure. K xperi- 
mental time required for fiiling and emptying the 
pipet is greatly reduced through capability of the 
commercial, nulJ-pressure diaphragm instrument to 
withstand high-pressure imbalance without damage. 
About 30 points on a pseudo-isochore are determined 
routinely in 8 hr by 2 or 3 men. 


Fas 


or 


3. Sample Preparation, Handling, and 
Analysis 


By mass-spectral analysis, electrolytic hydrogen 
in clean steel cylinders at 2,000 psig contains less 
than 5 ppm helium and 100 ppm air, respectively. 
The preparation system on the left of figure 1 has 
three functions. Impurities are absorbed at 76 °K 
in B-1 on 20 ml of 1.6 mm diam extrusions of ‘‘molec- 
ular sieve’’ silica [71]. The paramodification of hy- 
drogen is produced catalytically in B-2 at 19.6 °K 
on 10 ml of “30-100 mesh” particles of an iron oxide, 
batch 48-C, used in activity studies [72]. Small disks 
of sintered stainless steel retain the fine solids in their 
containers. Two stainless steel bombs, B—3, each of 
about 20 ml capacity, may be immersed in liquid 
nitrogen to provide a two-stage, thermal pump for 


boosting the parahydrogen pressure. Solidification 
of hydrogen at the catalyst is avoided by placing the 
catalyst ahead of these pumps in the train. The hy- 
drogen cylinders, preparation system, liquid hydro- 
gen transport Dewar, cryostats and vacuum pumps 
are separated from the instrument room by an eX- 
plosion-proof wall. Valve manifold C in the instru- 
ment bourdon rave as a doser for filling 
the pipet to selected densities and for controlling re- 
lease of the pipet sample to the gasometer system D). 

Systems B and C of figure 1 assembled with 
is-In. O. d. stainless steel tubing and high-pressure 
nidget valves with solid stems and polytetrafluoro- 
ethylene O-ring packings. Shop-fabricated gage con- 
nectors force a small, drilled, 
conical recess on the gage stem. 


< y 


roonh Uses a 


ure 


metal cone against a 
This special fitting 
avoid contamination of the gas. Bour- 
don gages for hydrogen must be phosphor-bronze, 
beryllium-copper, or type 316 stainless steel [75]. 

Analysis for parahydrogen is performed on gas 
bled from valve C-S at 100 ml NTP/min through a 
thermal conductivity instrument [58]. Calibration 
depends upon parahydrogen produced independently. 
In preliminary studies the rate of conversion of para- 
to orthohydrogen in the pipet at pressures up to 4,500 
psi was in the order of one percent per day at 76 CK 
and 25 to 30 percent per day at 275 to 300 °K. The 
vapor pressure of parahydrogen was measured at the 
beginning of a number of PVT runs [66]. Following 
one of the routine, 8-hr PVT determinations, the 
presence of parahydrogen again was confirmed by 
a Vapor pressure measureme nt. 


is necessary to 


4. Cryostat and Piezometer 


low 
14, 


Descriptions and summaries 


temperature 
techniques are available [11, 12, 


13, 38, 46]. 
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Figure 2 shows a cross section of the lower portion 
of the eryvostat. A 25-ml pipet sample cavity, 
5/8-in. diam, is bored into a two-in. diam, solid 
copper evlinder 8.5-in. long. To avoid oxidation, 
the plug is brazed with helium in the cavity. Metal 
stress is about one half the internal pressure. Helical 
grooves on the evlinder carry a 32-gage, 500 Q 
constantan wire heater. The platinum thermometer 
is tinned and cast into the pipet with Rose’s alloy. 
Above the pipet on the supporting, thin-walled, 
stainless steel tube, a guard ring thermally tempers 
electric wires leading to and wound on the pipet; 
the ring is automatically controlled at pipet tem- 
perature by a gold-cobalt' versus copper thermo- 
couple [40, 44] and a 36-gage, 100 Q constantan 
wire heater. A cold ring, integral with the refrig- 
erant tempers electric wires entering from 
room temperature. Before it was mounted and 
wound in the the bundle of 36-gage, 
low-level wires was constructed as a unit with all 
wires lightly varnished together parallel in a plane 
to insure good thermal contact on the rings. The 
levels of cooling through the reflux tube which do 
not disturb the temperature-versus-pressure rela- 
tionships are determined experimentally. About 10 
milliwatts cooling normally is used. 

An evacuable copper can, soldered to the tank 
with Rose’s alloy, encloses the pipet. A Dewar 
immersed in an 


tank, 


cryostat, 


vacuum jacket open Dewar of 
liquid nitrogen protects the can and the tank, 
A permanent liquid hydrogen transfer line with 


built-in valves enters the cryostat refrigerant tank 
through an O-ring seal, with the amount of refrig- 
erant in the tank shown by an electronic level- 
indicator instrument. Suitable valves and_ safety 
devices on the vent from the tank permit the use 
of a rotary oil pump to attain the triple-point 
temperature of the refrigerant. 
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5. Measurement of Sample Volume 


5.1. Normal Pipet Volume 


This volume is required for computation of 
density. By ignoring the adjustments contained 
in eqs (10.1-1) and (10.3-1) there is obtained the 
rough approximation for sample density, 

ot(PoIRT,)CV | Ve); (5.1-1) 
where P,,, T, and V, refer to the gasometry system 
and V, is pipet volume. Whereas the ratio V,/V, 
dominates the computation of eq (10.3-1), it is a 
practical convenience to perform absolute 
brations. 

The pressure and temperature-dependent volume 
of the pipet is required. Connection to the pipet 
cavity is made with a 3-in. length of 0.35 mm 
stainless capillary brazed to the pipet drill hole, 
figure 2. The pipet was briefly tested to 8,000 psi. 
\ midget valve was attached to the capillary, 
the pipet evacuated, and freshly boiled water 
admitted to the pipet. The three determinations 
yielded an adjusted volume of 25.83+0.016 cm* 
at 25°C Volume at room temperature after 
mounting in the cryostat was determined by ex- 


cali- 


panding hydrogen from the pipet at about 1,000 
psi into the gasometer system of section 6. Seven 
determinations yielded adjusted volume 25.854 


+ 0.037 em’ at 25 °C. Following two series of high- 

density experiments for PVT data to 5,400 psi, 
low-density experimental results violated general, 
low-density, limiting gas behavior. Six redetermi- 
nations of pipet volume by gas expansion vielded 
normal volume 25.913+0.009 cm* relative to the 
gasometer volumes, in agreement with deductions 
from low-density behavior. The first series of 
PVT experiments is rejected, and the increase of 
molal volume by factor 1.00263 is applied to all 
other data. The annealed pipet is thus assumed to 
have deformed in With this correction, the 
estimate of error in normal pipet volume is 0.5 
parts per thousand, 


use. 


5.2. Elastic Stretching of Pipet 


Circumferential and longitudinal stresses in a 
closed-end cylinder are, respectively, 

S, kP/(R 1) and S, Ig (R? 
where P is excess internal pressure, /? is the ratio 
of external to internal diameters, and k incorporates 
both end-effects and relative behavior of thick 
walls. With modulus of elasticity /, relative 
diameters, lengths, and volumes are, respectively, 


D/D,=1+kP/E(R—1) (5.2-2) 
L/L,=1+P/E(R?—1) (5.2-3) 
V/V.—1=(P/E)[((?—1)714-2k(R—1)74]. (6.2-4) 








The temperature-dependence of EF 
rs. 8 


from [73] in 


E=1.15 - 10°{1—4.35 - 10-*7], 


atm, (5.2-5) 
and the behavior of k by comparison with [70, 74, 79] 
Is 

I: a(R l | R25 1)/(R? 1)+m| 1. 
by introducing R=3.2, a 
m=0.3 for Poisson’s ratio. 
above estimates, therefore, 


1/3 for end-effect and 


A=V/V,=141- 10-144.35-10-'7]-Paim. (5.2-7 


5.3. Thermal Contraction of Pipet 


Thermal expansivity of copper [78] in the form 
a=(L,—L)/L, as a function of temperature is used 
to compute a table of V/V temperature. 
For use with computing machines, this is fitted by a 
quadratic polynomial on the Kelvin temperature 
SsCcaLC, 


versus 


V,=25.913 em’: a=0.990069: b=1.667 - 107°: 


e=Lsl7 - 107’. (5.3-1) 


6. Measurement of Amount of Sample 


6.1. The Gasometer 


The fluid contained in the pipet, the capillary, and 
the diaphragm cell is released into the gasometer 
system PD (fig. 1). This is a set of spherical glass 
standard volumes in a precision water-bath thermo- 
stat, and a valve manifold with precision mercury 
manometer [48] in a circulating-air cabinet.  Ele- 
vation of each arm of the leveled manometer is read 
to +0.05 mm by a highly reproducible lightbeam 
and photocell arrangement. A standard volume is 
selected such that the pressure is in the upper range 
of the manometer; i.e., 500 mm or greater. The 
instrument room temperature is controlled to about 
+0.5 °C. To compromise between low obnoxious 
ion and adequate pumping speeds, the mani- 
folds are \-in. copper tubing of 0.19-in. diam bore. 
For leak-free behavior, bellows-sealed valve 
and gaskets are a plasticized vinyl chloride polymer. 
Bellows volume increase per single turn is 0.15 em*, 
This adds to standard volumes of 10 liters or more 
used for most determinations. Jonization 
serve for vacuum leak testing. 


seats 


gages 


6.2. Gasometer Calibration 
Volumes of the 1- and 2-liter spherical glass flasks 
are determined by weighing them with water on a 
3-kg analytical balance against class Q weights. 
The 6-liter flask is weighed on a 20-kg single-beam 
platform balance against calibrated weight — slugs. 
Adjustment is made for air buoyancy, water density, 


Within the accuracy of 


and thermal expansion of the borosilicate glass [76]. 
Volumes at 300 °K 
of table 1 


are given by the second column 


TABLE 1 a ation of qasometer flask 1 ( 
Nor \\ Gi \ ned 
vo ‘ 
11 4.4 4. 45 4.1 
21 Ot 
61 120 i 6, 427 
I ¢ 0 21, 29 
I 2 Lil a. 
After determination of all manifold) volumes in- 


cluding the manometer, the 2-liter flask is calibrated 
by Yas expansion relative to the 1-liter flask; each 
successive flask then is calibrated relative to the sum 
of all smaller flasks With all adjustinents the 
results shown in the third column are obtained. 
Absolute volume standard is selected aus the sum of 
volumes of the 1- and 2-liter flasks by water weigh- 
ing, and the ratio of two volumes by 
expansion is selected for assigning the values shown 
in the last column.  Pressure-dependence of the 
flask volumes was computed from physical dimen- 
sions and mass of the individual flasks and the modu- 
lus of elasticity [76]. The nearly uniform § results 
were confirmed experimentally upon a flask filled 
with water extending into a capillary neck, 


these vas 


(1/V)(dV /dP)=1-107* atm 


6.3. Adjustment for Capillary and Cell 


The amount of sample 
amount behind valve 4 of 
gasometry, less amounts in 
diaphragm cell under conditions of measurement 
foreach P-T point. The rapid experimental method 
is justified provided these latter amounts can be 
estimated with sufficient accuracy. This is attained 
with a diaphragm cell volume which is small relative 
to the pipet, and by utilizing known and estimated 
gas-imperfection behavior of normal hydrogen. 

A volume of for capillary tubing and 
diaphragm-cell exterior to the top of the cryostat 
and up to the pipet valve C-4 is obtained by weighing 
mercury out of the tubing and by finding differences 
with suitable gas-expansion methods. The quantity 
of hydrogen in these obnoxious volumes is computed 
from their temperatures, the pressure, and the 
compressibility factor for normal hydrogen [65]. 

The volume of the 71-cm-long, transition capillary 
from the top of the ervostat to the guard ring is 
V.=0.0672 cm*, or 0.26 percent of pipet volume. 
The amount of fluid in the capillary is computed by 
calculating the temperature 7, corresponding to any 
position at a fraction + of capillary length from the 
ring at pipet temperature 7, to the top of the 
cryostat at room temperature 7), from thermal 
conductivity A of stainless steel [77] by the relation, 


in the pipet is the total 
figure determined by 
capillaries and in the 


0.655 em 
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(6.3-1) 


The number of gram moles of fluid in the capillary is 


N.=(PV «RT .)0 (6.3--2) 
where 
6=T dz/T,Z (6.3-3) 


Indices o refer to the rime or pipet, r to relative 
position along the capillary, and Z=Pv/RT where 
v molal volume. By known [23, 64, 68, 69] and 
extrapolated [3] data for Z, the capillary factor @ for 
hvdrogen is determined at a sufficient number of 
sensitive points for smoothing, taking note of liquid- 
vapor discontinuities; results are summarized by 
figure 3. The ratio of amount of fluid in the tran 
sition capillary to amount in the pipet is necessary 
for an estimate of required accuracy. This ratio is 
(V/V) - @, where ¢=Z,@ is summarized by figure 4. 
Due to the small value of V./V,, a ten percent error 
in @ would correspond to a relative error well under 











0.03 percent in density of the pipet sample at 
temperatures <100 °K. 
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7. Measurement of Pressure 


7.1. Piston Gage and Diaphragm Cell 


The commercial, dual-range, precision, dead -weight 
and null-detector operate in the oil-filled 
f-manifold, figure 1 [52, 53, 57,59]. The diaphragm, 
with adjustable electric contacts to indicate null 
position, is centered by use of an open-tube manom- 
eter. The oil pressure is adjusted to the experi- 
mental gas pressure with the screw-type oil press P. 
Prior to placing weights on the piston gage, the 
approximate pressure is read from bourdon gages. 

Piston diameters of about 0.182 and 0.407 in. are 
measured with a light-wave micrometer to +5 
millionths of an inch. Tolerance on the weights, 
determined by comparison with class 5 standards 
[54], decreases from 0.05 percent on the 1 psi weight 
to 0.002 percent on the 1000 psi weights. Adjust- 
ments to gage readings are made for temperature, 
pressure, acceleration of gravity, and barometric 
pressure. Observed sensitivity of pressure measure- 
ment with the small piston varies from under 

0.05 psi below 500 psi to under +0.5 psi at 5000 
psi. 

As compared with known diaphragm-detectors 
(49, 61], the commercial instrument used here 
withstands extreme pressure imbalance without 
damage, but exhibits temporary hysteresis effects 
up to about +0.05 psi. A measured shift of null- 


cage 


point with absolute pressure is appliel by the 
relation 
Po. Fipces T (1.00002) P94). ec 1) 


Table 2 summarizes estimates of accumulated errors 
at different pressures. 


TABLE 2. Estimates of accumulated pressure 
errors 
2 A) M1) 
(rage pressure, atm 
Error, atm X10 

Barometer 0. 50 0.50 0.50 
Diaphragm cell 1. 00 4. 00 1. 00 
Piston-cyvlinder diameter 0. 06 0. 6 6. 00 
Calibration of weights 67 1. 52 1. O7 

23 6. 62 14. 57 
Relative error (6P/P)-10 2.6 0. 33 0. 07 


Thermal equilibrium of the fluid sample, following 
a step-increase of temperature, is found by pressure 
observations to be established nearly as rapidly as 
the heating rate of the massive pipet. Occasional 
pressure balances are taken an hour apart to confirm 
this inference from the steady behavior of single 
observations. 


7.2. Adjustment for Capillary Column 


The column of cold fluid standing in the transition 
capillary tube produces an additional pressure P, at 
the pipet, 


eA gi 


IMLA/RT, 








using symbols of section 6.3, pressure conversion 
factor J, molecular weight M, and capillary length 
L.=71 em. Since JML,/R=1.68 - 10~* °K for hy- 
drogen, the adjustment for capillary column pressure 
is negligible. 


8. Measurement of Temperature 


8.1. Method and Instruments 


The potentiometric method of figure 5 is employed 
with a 25-ohm platinum resistance thermometer cali- 
brated by the NBS Temperature Physics section and 
a six-dial microvolt potentiometer [41, 42, 46]. The 
low-impedance, unsaturated standard cells are pro- 
tected from occasional extreme temperatures ol the 
laboratory by two concentric boxes of 1/2-in. alumi- 
num in a plywood case [43]. The upper potentiom- 
eter current supply is stabilized by placing a lead-acid 
battery across an isolated, rectified, and filtered a-c 
source [60]. A 25-v, 200 amp-hr, thermally insulated 
and shielded battery, of a similar, low internal- 
discharging type, supplies the stable thermometer 
currents which flow in the calibrated precision stand- 
ard resistors. The microvolt switch is wired so that 
the thermometer potential connections to the poten- 
tiometer can be reversed when the thermometer 
current is reversed, Microvolt galvanometer ampli- 
fiers with floating input circuits are coupled to 50-ya 
panel galvanometers in suitable damping and range- 
changing circuits. To minimize spurious effects in 
the thermometer and thermocouple potential circuits, 
the continuous copper wires from the pipet in the 
cryostat are joined by a special solder [45] to un- 
tinned, copper cable wires near the cryostat on copper 


| 





are observed by use of an internally-shorted emf 
position on the potentiometer and are bucked out 
at the galvanometer amplifier. When the cryostat 
is cooled with liquid hydrogen, a 0.35-uv spurious 
potential, which is nearly independent of pipet tem- 
peratures up 410 °K and increases with higher 
temperatures, appears across the thermometer po- 
tential leads in the absence of thermometer current. 
Potential measurements with reversed thermometer 
currents show the effect to be independent of current. 

Thermometer currents are adjusted to exact in- 
tegral milliampere values. To a corresponding, 
tabulated, thermometer potential is added the spuri- 
ous thermometer potential to obtain a potentiometer 
dial setting yielding an automatically-controlled, 
integral-valued temperature. A thermometer current 
of 5 ma is used at temperatures below 35 °K; 2 ma, 
from 35 K; and 1] 


to 


ma, at 75° and above. 
Thermometer heating by the 5 ma current produces 
an error of 0.0007 °K at 40 °K. All measurements 
occur in the low range of the potentiometer. 

The specified error limit of the potentiometer low 
range is 0.01 percent of reading plus 0.02-uv. The 
calibrated value of the l-ohm standard resistor is 
0.99999+-0.00001 ohms. Apparently satisfactory 
agreement was found between the platinum ther- 
mometer calibration tables and measurements with 
above procedures for the ice-point and for vapor 
pressures of high-purity nitrogen and of parahydro- 
gen, as described in section 11. 

The accuracy of temperature measurements is 
summarized by table 3. It is based on the ther- 
mometer characteristic and currents, the potentiom- 
eter specifications, and a latitude of 0.05 yv for 
bucking out spurious potentials. Thermometer 


to 75 


lugs in an iron-shielded, heavy aluminum box. The | Calibration uncertainties are not included. 
crvostat wires and shielded cables run in iron con- TABLE 3 Estimates of accumulated ten perature errors 
duits. The cryostat and all electrical instruments 
are interconnected and grounded to the earth with — —— Et r, °K - 10 Ter 
heavy copper conductors. lr, °K aes 
7 i R External 
+005 pv 
8.2. Procedure and Calibration 
0.049 “ 04 1.2 135 170 
- ° ° — " a 11 ( 7 RO () BR 74 
When the instrument room air conditioner is 10 2 107 234 fi2 1: 19 
started, spurious potentials up to +0.1 microvolt ka 4.6 109 “+ 2 09 3 
arise in the galvanometer circuit. These potentials 
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236 


9. Regulation of Temperature V,=vol of lines, gasometer to C-valves, including 
gage, 22.37 cm’, 


Precisely controlled electric heating compensates 


for crudely adjusted, spontaneous cooling (section 4). V.=vol of glass bulbs used, em*, 

The pipet temperature can be controlled [29] from 

a thermocouple, relative to the temperature of the V, vol of water-bath manifold, 13.00 em*, 
refrigerant tank (upper left of figure 6), by a wide 

range microvolt bucking circuit, +7 | microvolt nam- V vol of pipet, cm. 

plifier [56] and a d-c power-regulator [50]. For all 

PVT measurements the pipet temperature is con- | V,/7,=25.913[0.990069 T,.+1.667.10~' 

trolled from the platinum resistance thermometer 131730" « TJ, 10.1-5) 
and potentiometer. There is no perceptible noise at 

0.0002 °K sensitivity, The accuracy of controlled | 7) (at Ty and P,,)=1.000044+0.000555P,,, 
temperature is equivalent to the described accuracy 10.1-6) 
of measurement. ‘To maintain a steady temperature 

during other operations with the potentiometer, the Z... (at JT, and P..)=1.0007 


power regulator input signal switch is opened and 
the steady-state heating carefully adjusted by the 
power regulator control before undertaking such op- 


1055/7” for T,>40 °K. (10.1-7) 


erations If required, the above thermocouple con- 10.2. Pressure 
trol can be used. The guard ring temperature is 
controlled with similar instruments. The pipet fluid pressure P in international atmos- 


pheres is 


10. Computation of PVT Observations 


P=0.0680457P,+-P, (10.2-1) 
10.1. Gasometry 
a where 
The number, NV, of gram moles of hvdrogen con- 
tained behind pipet valve C—4 is calculated by the P,=deadweight oil gage pressure in psig, 
relation, P,—barometric pressure in atm, 
RN/P (EV.+V.)/7.4+ V./7 P,=k(qg/gs)[1—2.4 - 107-°(t,—20 °C)] [1 
, ey , T) y, VT Z (10.1-1) 4.2-10°-SP,]-[Nominal wt psi]. 10.2-2) 
where P,, is mercury manometer pressure in atmos- k=1.00002, eq (7.1-1), 
pheres, (q/g 979.615/980.665, 
‘ t,—deadweight gage temperature, °C, 
R=82.057 cm* atm/g mol deg K, P,=KgD - Ly as in eq (10.1-2), 10.2—3) 
P.—KeDI (10.1-2 L,—barometer height, em, 


t,== barometer temp, °C. 
K=0.986923.10~° atm em?/dyne [4], 


10.3. Density 


g=979.615 dvne/gram in this laboratory room, : 
; : Density, o,¢g moles em®, is calculated by 


D 13.5948 — 0.00245 .- f ,giem [4], (10.1 a A (P R){I 9 T +1} JZ, 1 10.3-1) 
t, manometer mercury temperature, °C, 

where 
L,—=Manometer mercury height, em (brass seale), T,=temperature NPD, °K, 


I-=glass volume/pressure coefficient T.,=pipet temperature, °K, 
1+10-'*(P,,—0.825), (10.1-4) ay 
25.913 A[0.990069 + 1.667 - 107°7 
T,=air-bath manifold temp, °K, L3i7 - 107° 73k (10.3-2) 


ow . y } 0672 e 
l’,=NPD temp, Ix, \ sarcaiaiaaa soll 
Va, volume of N PD and warm capillary 0.655 em’, 


T ipet temp, CK, . ° a 
PH Z,, for gas at P and room temperature in V4, is 


T..=water-bath temp, °K, Z, =0.99990 + 0.000623447 P, (10.3-3) 
. Fr 9 
V vol of air-bathed manifold and manometer, hy OF LO.8— 8 
59.66 em’. 6, eq (6.3-3) and figure So. 








11. Checks on PVT Calibrations 


11.1. Thermometry Fixed Points 


Vapor pressures of a commercial, high-purity 
nitrogen and of parahydrogen were measured prior 
to the PVT experiments. Examination of experi- 
mental technique later revealed failure to compensate 
some of the spurious potentials. Determination of 
these potentials provided adjustments, estimated 
accurate within 0.005 °C, which have been applied 
in preparing table 4. 


TABLE 4. Thermometry fixe d 
point observations* 


P,mm . Tehe Tiit AT 


Vapor pressure of nitrogen 





625. 05 75. 732 0. 002 
625. 42 75. 741 —(). 006 
626. 03 75. 744 +0. 001 
631.75 75, 825 +0). 009 
631. 68 75. 826 +0. 011 
Vap ire of parahydrogen ¢ 
626. 33 19. 639 19. 632 +0. 007 
626. 72 19. 645 19. 634 +0. 011 
626. 9O 19. 650 19. 635 +. O15 
628. 18 19. 656 19. 642 +-0.014 
629. 12 19. 656 19. 646 +. O10 
*All temperatures in deg K on NBS 
1955 seale 
* Mercury column height reduced to 0 
C at standard ¢ \ 
bG. T. Armstro 





¢ Hloge and Ar 


11.2. Elevated Vapor Pressures 


Initial PVT points of several runs were taken at 
saturation conditions to check instrumentation. If 
the substance is nearly pure parahydrogen, the 
observed pressure-temperature relations should agree 
with those of Hoge and Arnold. Table 5 presents 


} 


TABLE 5 Elevated va por-pressures Obse¢ rvations 


Pods, atn Data from Hoge 
und Arnold 4 

> 
r°K I] Ill P, atm dP/dt 
22 1.612 1.614 0. 417 
23 2. 069 2.071 44S 
24 2.614 587 
25 3. 246 3. 249 685 
26 3 O78 , O83 3. OR5 700 
27 4. 832 4.829 4.8352 904 
28 5. 785 5. 792 5. 795 1. 026 
29 6. S86 6. SSY 1. 162 
30 8. 121 1. 304 
31 9. 499 9. 502 9. 5OO 1. 459 
32 11. 051 11. 047 1. 643 
32 11. G51 11. 047 1. 643 
32 11. 052 11. 047 1. 643 

* Data of Hoge and Arnold interpolated to N BS 1955 temperature scale 


observed vapor pressures of parahydrogen from the 
series II compressed liquid runs and from the series 
III low-density runs in the vapor and gaseous region. 


Data of Hoge and Arnold are interpolated to the 
NBS 1955 temperature scale of the present work. 


Temperature deviations corresponding to pressure 
deviations in 
28 °K to 0.0030 deg at 32 °K. 


series IL range from 0.0098 deg at 


The series [1] range 
is from 0.005 deg at 22 °K to 0.003 deg at 32 °K. 


11.3. PVT Data for Normal Hydrogen 


The first PVT run was made with normal hydrogen 


. ‘ - . mn 
for comparison with known data, table 6. The 
: | ABLE 6. PVT ohse rvations on normal h ydroge n 

iy 8 P, atm Vv, Cl Z 1.0026+ Zexy Z 
28 30. 869 30. 443 0. 4090 0.4101 0. 40934 
30 45. 357 30. 460 wi) 4 O27 56168 
32 59. 738 30. 476 6933 6952 6945" 
36 SS. 443 30. SOS 9134 9160 914* 
10 116. 969 30. 538 1. ORS83 1. 0911 1. 0902 
45 151. 884 30. 574 1. 257¢ 1, 2609 1. 259° 
0 186. 213 $0. 607 1. 3892 1. 3928 1. 3928) 
* Friedman and Hilsenrath, private communication 


Johnston and White, Tran 
*Nonlinear interpolation 
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known data have been interpolated linearly except 
where indicated. Present results for Z in column 4 
are multiplied by 1.0026 in column 5 as an adjust- 
ment for the pipet deformation which occurred at 
an unknown time during the experiments. Devia- 
tions of the former results are within 0.07 to 2.6 
parts per thousand; of the latter within 0.00 to 2.2 
parts per thousand 


11.4. Virial Coefficients 


Corresponding to the virial equation, 


Z Pv RT | B v C'/v" t a a (11.4 1) 
the data function, 
o=RT(Z—1)v=RTB+RTC/v ~ 2+. (11.42) 


exhibits satisfactory linear dependence on density 
at each temperature when the recalibrated pipet 
volume, V)=25.913 em’ is used for computation of 
experimental densities below critical. Unsmoothed 
values of RTB and RTC, determined by least squares, 
are given by table 7. These new data are equivalent 
to an equation of state for parahydrogen at densities 
up to 9 g mol/liter. For comparison are given 
values for normal hydrogen calculated from the 
smoothed table 19 of H. W. Woolley et al. [69], and 
values for RTB only, smoothed by us from the 
“smoothed” table CIV of A. S. Friedman [19]. 
Negative third virial coefficients in the latter table 
below 50 °K suggest that corresponding second virial 
coefficients are small in absolute value. Our 
closely-spaced data suggest a so-called lambda- 
point maximum in the third virial coefficient at the 
critical temperature The criteria of its validity 
have not been investigated. Concerning comparison 
of results, it recently has been shown that second 
virial coefficients for parahydrogen may be about 1 
percent smaller in absolute value than for normal 
hydrogen at liquid hydrogen temperatures, in 


too 
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agreement with 
study of the 
correction of 


these results [63]. Quantitative 
¢g-density behavior hefore and after 
the pipet volume suggests accuracy of 
volumetric calibrations within one part per thousand. 


rABLE 7 Virial coc flicients, parahydroge n 
Units of em itm, deg K 
RTB-10 RTC.-10 
je I Er] 
i 229 2. 247 2. aw 3. 053" 
Is 3. 432° 
13Y} 2. 159 2. 11 3. 273% 
OOS 3.519 
ss, Osu 2.07 (V2 3. 047 
; (3 3. 826 
( s 1, Gut ) 4. YAS 
| V52 4. (21 
1.917 1. ¥21 1. 8S 4. O06 
1. 8S 4.124 
$4 1. S48 1. 849 1.81 4.110 4.171 
$. 1. S14 4. 103 
M 1. 780 1. 781 4.073 4. OSY 
5 1. 747 4.043 
$5 714 1.71 1.70 4. 020 4, 020 
$4 1. 683 4 GOS 
4( 1. 651 1. 653 1.64 3. 966 3. W62 
12 1. 591 1, 591 1. 59 3. 934 3.914 
14 l $3 4. 920) 
it 1. 478 1. 473 1.49 3. 922 3. 846 
1s 1. 423 3. ¥2Y 
‘ 1. 369 1. 360 1. 40 3. 926 3.811 
1. 239 1. 224 2s 3. 971 3. 804 
60 1.11 1. lf 4. 074 
‘ 0. 99 4.143 
( 8 “4 4. 232 
7 4.414 
st 4 t 1.514 
s $1 4. 587 
" 117 1) 4. OSS 
) $14 1.849 
1K f | On 
Ii. W. Woolley et J. Research NBS 41, 379 (1948) RP1932, normal I 
A. S. Friedman, Ohio State dissertation (1950), normal HH 
*Based on only two data points 
Basedf{on only three data point 


12. Comments on Overall PVT Accuracy 


Requests for a statement of overall accuracy of 
data in the form Z=Pr/RT are commonplace. 
This accuracy depends not only on calibrations, but 
also on all the adjustment computations and the 
partic ‘ular conditions defining the state of the fluid. 
The relative error in Z is the sum of relative errors in 
experimental variables, 


6Z bP ; ov éT" 
FA P T y T 


Estimation of pressure errors is indicated by section 
7.1. and table 2. The complexity of error estimates 
for molal volume, v, is appreciated by reference to 
eq (10.3-1). In view of sections 5 and table 1, 
and the unpublished details of calibration and ad- 
justment of obnoxious volumes, we venture the 
experimentalist’s opinion that error in » generally is 
under two parts per thousand. Temperature errors, 
in view of tables 3, 4, and 5, are within 0.02 °C 
except near 14 °K, comparable to uncertainty in 
the low temperature scale [38]. Value of the gas 


constant used for present computations is R=82.057 
em? atin/g mol deg. It rey be noted that at high 
densities, anv isotherms of pressure as a function of 


density are highly sensitive to errors in density, due 
to the relatively incompressible state of the fluid, 
a situation giving spurious indications of low pre- 
cision in the measurement of pressure. 
13. Calorimetric Apparatus 
13.1. Introduction 
Measurements of the volume, mass, and tem- 


perature of a fluid sample are necessary to define the 
state corresponding to each specific heat determina- 
tion. For high pressures the PVT pipet must be 
massive so that its volume will be calculable. The 
PVT data are obtained first because they are used 
for ealibration of the volume of the thin-walled 
calorimeter as a function of temperature and _ pres- 
sure. Calorimetric methods have been reviewed 
[29] and the calorimetry of some compressed fluids 
described [28, 30, 33, 34, 36]. The purpose of this 
section is to sami a gts description of some 
improved instrumentation and of the integration of 
calorimetric with PVT apparatus. Calorimetric 
calibrations will accompany the publication of results. 


13.2. Cryostat and Calorimeter 

The calorimeter resides in an additional cryostat 
of basic construction identical with that for the PVT 
pipet, with the high-pressure hydrogen capillary, 
vacuum, and liquid-transfer plumbing in parallel 
through suitable valves. Alternation of experiments 
requires no dismantling. Figure 6 is the calorim- 
eter region of the crvostat, and figure 7 is the 
calorimeter section. All parts are copper except for 
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Figure 7. Calorimeter section 


sphere, capillary tube, and the guard-ring support 
tube. The 2-in. diam, type 316 stainless steel. 
spherical calorimeter has a nominal wall thickness of 
0.06 in. It is fabricated from spun hemispheres, 
welded together, with addition of a 1-in.-long neck of 
\4-in. tubing. The spherical shape is improved by 
expansion with water pressure. Yielding to a per- 
manent volume increase of 10 percent was observed 
with a pycnometer surrounding the calorimeter. — It 
occurred at from 5,000 to 6,400. psi, 
corresponding to metal stress roughly 8.3 times these 


pressures 


values. Elastic behavior following this deformation 
was AV/VAP=1.59-10 atm~', corresponding 


roughly to a circumferential modulus (3d/4?) 
VAP/AV=23-10° psi, where d is diameter and ¢ is 
wall-thickness. After electroplating copper inside 
and outside to about 0.2 mm thickness, voiding 
sulfate electrolytes, the normal volume is 72.35 
em’, determined by gas expansion from the cal- 
orimeter in the cryostat. The 100-ohm, constantan 
wire heater is varnished directly onto the sphere, 
and shielded by the lightweight calorimeter case, 
which serves also as the wire spool for thermally 
anchoring thermocouple and thermometer lead- 
wires. The '{s-in. o.d. stainless steel capillary has 
a bore diameter of 0.0345 em. This relatively 
large bore is essential for rapid filling and emptying 
of the calorimeter. Volume of this capillary, 
including a valve at the head of the cryostat, is 
0.0SO0 em’, 


Adiabatic shielding is fn two parts: the guard 
ring, and attached, lightweight shield can. The 
ring tempers all electric wires, pressed into indi- 
vidual, longitudinally machined grooves. These 
wires lead from similar grooves in the cold ring and 
lead to the calorimeter case. Separate heaters are 
wound on the guard ring, the evlindrical surface of 
the shield and as a plane spiral on the bottom of the 
shield can. The two shield heaters operate in a 
divider circuit from a single regulator. These 
heaters are automatically controlled from thermo- 
pile 7 and thermocouple 2 respectively of figure 6. 
Thermocouple 3 serves for manual adjustment of the 
side-to-bottom heat divider. Thermocouple 4 
serves for automatic temperature control of shields 
and calorimeter, during gas-expansion calibrations 
of the calorimeter volume, with reference to the re- 
frigerant tank temperature. 

To facilitate control of the calorimeter tempera- 
ture during some calibrations, a small heat leak is 
introduced to the present apparatus. Thermal 
contact between electric wires and tempering ring 
is reduced to obtain the calorimeter cooling drift de- 
scribed in 13.4 The calorimetric method 
therefore may be described as quasi-isothermatl. 


section 


13.3. Control and Measurement of Heat 


The control circuits for calorimetric heating are 
indicated by figure 8. Although they are con- 
ventional in principle, a useful simplification is 
application of an electronic battery: a 0.5 amp, 
transistorized, d-c power supply, continuously ad- 
justable to 0.02 v from 0 to 32 v, with regulation 
and stability better than 0.002 v. The rms ripple is 
0.001 v. The battery source impedance, PR, of the 
figure, is readily adjusted equal to the calorimeter 
heater resistance when the batters potential in- 
dicated by P-—1 is twice the heater potential indicated 
by P-2 [31]. To provide time for potentiometric 
measurements of the heater power, and for shield 
temperature adjustments, a temperature rise of 
about 2° in a fixed interval of 10 min at all temper- 
atures is realized by adjusting the battery potential. 
This adjustment does not affect R,, in contrast to 
behavior with conventional batteries and rheostats. 
Relay with direct current to 
reduce induction of spurious a-c effects in the low- 
level circuits. 

Potentials are measured with the 6-dial microvolt 
potentiometer and with a type K—3 universal poten- 
tiometer, each indicating through an electronic 
galvanometer [56], as in figure 5. Potentials to be 
measured are connected to each potentiometer by 
high-quality, commercial, microvolt selector switches. 
Four measurements of heater potential and current 
are made in 10 min. Heater power is constant to 
better than 1 part in 10,000. The standard electric, 
clutch-type, synchronous stopclock utilizes local, 
60 c/s, powerline frequency for the time standard. 
The NBS Time Standards Section finds that this 
frequency usually is accurate and constant to well 
within 0.03 percent 


coils are operated 
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FIGURE 8. 
13.4. Measurement of Temperatures 


The calorimeter platinum resistance thermometer 
is calibrated and the thermometric technique checked 
described for PVT determinations. Imperfec- 
tions of adiabatic shielding vield empty calorimeter 
drift rates from —0.0009°/min at 20 °K to —0.0019 
min at 100 °K. These are linear in time. 

Five temperature observations in 5 min precede 
a heating interval. Nonlinear equilibration fol- 
lowing the interval requires a few minutes, and is 
observed with a recorder. five similar temperature 
measurements follow. Temperature-time behaviors 
are extrapolated mean time the heating in- 
terval, by a least-squares computer program, at 
which time tuken the calorimetric temperature 
Increment 


as 


to of 


is 


13.5. Control of Shield Temperatures 


The different methods, isothermal and adiabatic 


shielding, are useful for different purposes [87]. The 
adiabatic method is selected here for rapid adjust- 
ments in the range 20 to 100 °K where thermal 
transport ts almost entirely conductive. Gold- 


cobalt versus copper thermocouples indicate temper- 
ature differences between the calorimeter and the 
guard ring and shield respectively. An individual 
regulating system is emploved for each of the latter. 
An improvement of present instrumentation, in- 
dicated by the bottom row of figure 5, is its simplicity 
relative current practices [27]. Control instru- 
ments are the microvolt amplifier [56] and the power 


Lo 


regulator [50]. Two factors contribute to satis- 
factory operation. (1) Thermocouples are so placed 
on shields as to give rapid but average response. 
The heater windings and powers are designed to 
match the areas and masses of metal to be con- 
trolled. (2) The power regulator contains a solenoid 
relay, peat: by the heating interval signal, 
figure 8, which automatically boosts the shield 


heating to a manually pre-set level. This eliminates 
the delay otherwise required for development of the 
temperature difference normally utilized for auto- 
matic control. The booster circuit does not affect 
the low output impedance of this regulator, thereby 


Calorimetric 


60 C/S, 0.03 % 


he al control Mie thods. 


avoiding the loss of power-gain and temperature- 
regulation otherwise suffered when resistances or 
other impedances are inserted in the heater circuit 
for boosting control. Because automatic resetting 
instruments are not necessary for the present work, 
the normal steady-state heating also manually 
adjusted at the power regulator. For calibration 
work, the shield versus calorimeter temperature dif- 
ferences may be recorded from the microvolt ampli- 
fier output signals by means of simplified strip- 
chart instruments. A typical chart shows that, 
with routine error in adjustment of the boost heat 
level, the guard ring-versus-calorimeter temperature 
difference peaks at 0.001 °C at the start of an interval 
and decays to within 0.0001° in The same 
behavior with opposite sign follows the ending of 
an interval. Relative change in the steady heat-leak 
over the copper wire bundle between gu: ard ring and 
calorimeter is negligible during a heating interval. 


is 


20 sec. 


The numerous machined parts for 
were made in the NBS, Boulder shop under S. 
Landis. The skillful assembly and silver-soldering 
of these parts as well as the construction of the liquid 
transfer and hard vacuum systems by Cryogenic 
Division instrumentalist W. R. Bjorklund greatly 
advanced the project. D. E. Diller placed all of the 
small electric wiring in the cryostats, intercalibrated 
the gasometer volumes, performed the gas- 
expansion calibrations of the pipet, checked fixed- 
points of the platinum thermometer, and measured 
null-point shift of the diaphragm cell in addition to 


the crvostats 


ols ass 


other important contributions. L. A. Weber as- 
sisted with the many details of placing the PVT 
apparatus in final operating condition. B.A. 


Younglove determined mechanical properties of the 
calorimeter and with D. E. Diller placed the calo- 
rimetry in successful operation. Hans M. Roder, in 
programing the computer work, devised subsidiary 
equations presented in sections 5.3 and 10. Con- 
tinuing encouragement from Robert J. Corruceini 
and Russell B. Scott have been essential in bringing 
this long task to completion. 
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The adsorption-desorption phenomenon which occurs when a multicomponent liquid is 


forced through a column of solid adsorbent is accompanied by an evolution of heat. 


The 


resulting temperature change occurring within the column may be monitored and correlated 


with the compositional changes occurring in the liquid phase. 


thermistor sensing cells is described. 


The construction and use of 


These cells are an integral part of the column, and 


they are maintained at appropriate temperatures by either constant temperature, liquid- 


circulating baths or by transistorized temperature controllers. 


Both methods of temperature 


control are described, and various applications of the cells are discussed. 


1. Introduction 


As part of a study of liquid-solid chromatography, 
it was desirable to develop means which would allow 
the process to be observed at various positions down 
the column length. The adsorption of a liquid upon 
the solid adsorbent causes an evolution of heat, and 
the displacement of this adsorbed liquid by a more 
tenaciously adsorbed liquid produces a further evolu- 
tion of heat. The heats of preferential sorption 
upon a wide variety of solid sorbents have been 
determined using a calorimeter of simple design.’ 
An application of this principle to observe the heat 
from the adsorption-desorption process as it takes 
place in the chromatographic column is described 
herein. 

A small thermistor, embedded in the column 
packing, is used to detect the temperature changes 
produced by the adsorption-desorption process. 
Since the magnitude of this temperature change may 
vary from, at most, a few degrees to as little as a 
tenth of a degree, the portion of the column contain- 
ing the thermistor must be kept at nearly constant 
temperature if good sensitivity is desired. This is 
achieved by a constant temperature, circulating 
bath or by electrical resistance heating with a suit- 
able controller. The thermistor, the brass block 
which houses and supports the thermistor, and the 
necessary heating jackets are referred to herein as 
the “cell”. These cells prov ide an in situ analytical 
tool for sensing and recording the changes in column 
temperature which result from the heat generated 
by adsorption. The inner diameter of the cell ts 
made the same as the inner diameter of the column; 
the packing is continuous through both the cells 
and the column; and there are no large obstructions 
to the liquid flow. The cells, described in detail 
later, are connected to the glass column by means 
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of a kovar-glass seal. As many cells as are desired 
may be placed at predetermined positions down the 
column length, and the column behavior may be 
monitored without loss of material and without 
disturbing the column continuity. Not only can 
the progress of the separation, as it occurs on the 
column, be observed, but the collection of fractions 
can be accomplished automatically, and qualitative 
estimates of the differential heats of adsorption can 
be obtained. 


2. Apparatus 


Several different types of cells have been con- 
structed. While they differ, depending on their 
intended use, they are similar to the one shown in 
figure 1. This cell is made in two sections, which 
are electrically insulated from each other by means 
of a Teflon gasket. A glass-coated bead thermistor 
(2,000 ohms at 20 °C) is centered between the two 
sections, as shown (D). It has been found that the 
2,000-ohm thermistors give adequate response when 
used one arm of a balanced-arm Wheatstone 
bridge. The platinum-iridium leads of the ther- 
mistor electrically bridge the two brass cell mem- 
bers. These sections are held together by spring- 
loading, which maintains a liquid-tight seal despite 
temperature changes. The column is packed with 
solid adsorbent so that it is continuous throughout 
the total length of the column, including the cell 
sections. The inner diameter is the same as that 
of the column so that the chromatographic separa- 
tion may be recorded without disturbing the ad- 
vancing liquid fronts. Figure 2 is a schematic 
representation of a typical column showing two 
different types of cells in a single column. The 
upper cell is the same as that shown in figure 1. 
In the upper cell, figure 2, water maintained at a 
constant temperature by an auxiliary water bath is 
circulated through a water jacket. For the lower 

design, an electrically- 


cell, which is of one-piece 
heated copper block is substituted for the water 


as 
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Ficure l. Te mperature sensing cell. 
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F. Kel-F insulator 
G. Brass washer 
H. Spring 
I, Screws ho halves together 
= ry . . . 
jacket. The block is maintained at constant 


temperature by the transistorized proportionating 


controller? shown in figure 3. 


3. Experimental Procedure 


An example of the chromatogram obtained with 
this apparatus is given in figure 4. The test mixture 
contained equal volumes of n-heptane, carbon 
tetrachloride, and benzene. The mixture was 
chosen to illustrate different degrees of compound 
separation. Ten milliliters of this mixture were 
poured onto the top of a silica gel (28-200 mesh) 
packed column. The column was 6 ft long and * 
in. in diam. The cell was maintained at 30 °C 
Ethyl alcohol was used to elute the test mixture. 


2 The author is indebt 
for designing this low-cost 
paper. The instrument 
operating temperature of 


» Malcolm Morse of the Instrumentation Division 

ontroller and for the privilege of including it in this 
provides control of +0.02 °C for several hours at 
60 °C, 
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FiGurgE 2. Typical packed column. 


I'wo cells are showt Che top one is a two-piece cell being thermostated by 





onstant-temperature circulating liquid. The lower cell is a one-piecefcell beir 
thermostated by an electrically heated jacket which is maintainedj{at constant 
temperature 

A. Kovar to Pyrex graded seal 

B. Wire serving as electrical contact 

C. Jacket containing thermostating liquid 

D. Cell as show fig. 1 

E. Electrical contact 

F. Kovar t ex vraded se 

G,. Brass tu ~ 

H. Control tor 

I. One-piece br ‘ 

J. Electrical } ing wire 

K. Sensing thermistor i on Os diam, in. length 

L. Glass insu 


The cell was placed at the exit end of the column 
that the refractive indices of the eluent could 
be compared directly with the recorded temperature 
Variations. However, it could have been placed 
in any position down the length of the column, 
or several cells could have been used in series so us LO 
follow the separation process. When a signal at- 
tenuation of 30 times was used, the record shown in 
figure 4 was obtained by using a 10-inv, poten 
tiometric strip-chart recorder. The refractive index 
curve (dotted) was obtained by measuring succes- 
sive 0.5-ml portions of the sample as it emerged 
from the bottom of the cell. Accumulative error in 
measuring the volume of these cuts is probably 
responsible for the phase shift between the two 
curves. 


sO 








tr 4. Discussion 


o ‘4 
AMA >—+ a eo 
a, ft | | 35 This cell has been used with several different test 
Te + i -.., | mixtures. In all cases distinct peaks were obtained 
+ be ti | «when the refractive indices showed that a separation 
ray be LY | _) of the components had taken place in the column. 
' fat SR In general, a broad peak is indicative of incomplete 
i Ce "Slut r- - separation, and the amount of tailing appears to be 
% ow i + Fi inversely proportional to the degree of separation. 


The height and broadness of the peaks are influ- 
enced by many factors such as degree of resolution, 
flow rate, thermal conductivity, thermal diffusion, 


Figure 3. Wiring diagram for transistorized proportionating 


controller. cell temperature, ete. Under constant column 
head-pressure, the flow rate is nearly constant, 
Resistance (Rj), 1 watt unless otherwise specified and the height and area under the peaks give val- 
1 10K 15 100 ohm uable information concerning the relative differential] 

_ ood aie heats of adsorption and quantity of material. 
a4 Ay ae The cell may be used to actuate an automatic 
6 800 ohm 20 510 ohm fraction collector. This can be done by using an 
: oo = to adjustable limit-switch on the recorder to rotate the 
bE coal a+" fraction collecting table when a new peak appears. 
11 500 ohm 25 3.5K For those materials which do not separate sharply, 
2. oo oS oe the heading and tailing of the peaks may be used to 
M4 1K 2 51K actuate the fraction collector. In this case, the 
Capacitors (( rate of peak decay or buildup could be used as the 
ante intl ore controlling signal. This method would have the 
wie advantage afforded by the taking of fraction as a 
function of changing purity, but a combination of the 
Sean’ ahi . first method together with the usual volume or 
4 ol 6 aa time-actuated mechanism of the fraction collector 
would achieve nearly the same result with much 

See ees eee simpler control. 
Hav +9 When either the sample or the adsorbent is elec- 
A trically conducting, the cell consisting of the single 
sas ° block should be used. The thermistors used in this 
1,2 800-ohm thermistors at operating temperature cell are of the glass-enclosed probe type and are 
Relay inserted through the side of the cell. The probe 
Sigma, 7AOY 3,500 SD-SIL. 3,500-ohm dual coils, polarized ean be sealed in the block with the usual stuffing 
box arrangement or by making a liquid-tight seal 
with epoxy resin or soft solder. Such a thermistor 
has a much higher heat capacity than the bead type 
and the recorded peaks tend to be less high. The 
, larger size of this probe does not seem to disrupt the 
NIEN ‘—- a flow pattern of the liquid. The problem of increased 
: \ F heat capacity can be minimized by using a low-gain 
en | d-c amplifier prior to recording the signal or by using 
; | thermistors of higher resistance. 


| \ The use of cells of this and similar types in con- 
junction with other analytical tools should make a 














| = ‘ z i. _ P 
systematic study: of liquid-solid chromatographic 
separation more feasible. 

TAN t | 
' ) 
WiLLLITER ENI The author is indebted to Miss Naney Beach for 
her assistance in the performance of early experi- 
FiGure 4. Recorded chromatograph. ments and to Herbert B. Lowey for his assistance 
in the design and fabrication of the cells. 
Refractive index plot is shown as dashed curve. Refractive indices of the pure 

components are shown at left of figure. ‘A’’ is the peak due to the heat of adsorp 

tion of n-heptane; “‘B”’ the peak due to carbon tetrachloride, “C’”’ that due to / 7-4 @) Ae f 

benzene and “D”’ that due to ethanol ’ { Paper H5¢ 4 4d) 
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f radiation shielding problems involve calculations of the response of an 


isotropic detector to radiation of arbitrary angular distribution from uniform finite plane 


sources 


Series expansion methods previously used for treating the rectangular source are 
here extended to include the circular disk source with detector off-axis. 
involve the family of integrals /’s (cos 6 dS/r?) P; (cos 6) and the integral f’s (dS/r?) exp ( 


These methods 
br) 


where @ is obliquity with respect to the normal to the disk surface S, r is the distance from 


an element of source area, dS, to the detector, and yu is the attenuation coefficient. 


Tabula- 


tions of the first type of integral facilitate use of Legendre expansion representations of 


radiation directional distributions. 


The second integral relates primarily to exponentially 


attenuated radiation from a plane isotropic disk source, but the expansion coefficient solution 
can be readily adapted to take into account a point isotropic buildup factor of the polynomial 


Ny 
orm >* 
form p 

i=0 
coefficients previously 
are presented 


1. Introduction 


In radiation shielding, photometry and related 
studies, a frequently occurring problem is that of 
calculating the radiation field resulting from a finite 
plane source of radiation with arbitrary polar angular 
distribution g(@). A number of authors have in- 
vestigated the response of surface-type [1]? detectors 
to rectangular and circular Lambert’s-law radiators. 
Smith and Storm [2] have derived a generalized 
series expansion method for calculating the response 
of an azimuthally symmetric detector off-axis from 
a disk source with arbitrary surface angular dis- 


tribution, in which evaluation of the series terms 
involves successive differentiation of the detector 
response to a differential element of the source 


area. The resulting formulas usually consist of 
triple or quadruple summations in which one sum 
can be obtained by recursion methods. 

In-a previous paper [3], a Legendre expansion 
method suggested by Berger and Lamkin [4] was 
applied to a rectangular source as part of a general 
study of the penetration of radiation into  struc- 
tures [5]. In the Berger-Lamkin method, the geo- 
metric relationship of detector to source is character- 
ized by a set of Legendre coefficients, p;,, over which 
corresponding penetration-dependent coefficients, gj, 
of a Legendre expansion of g(@) may be summed, 
This has some advantage over the results in reference 
{2} in that not only are the geometrical and penetra- 
tion parameters separated for economy of data 
tabulation, but use of the present Spencer-Fano [6| 
moments method provides most directly the gq, 


coefficients, rather than the actual angular dis- 
tributions g(@ 
Work supported jointly by the Office of Civil and Defense Mobilization 


the Navy Bureau of Yards 


ind Docks, and the Defense Atomic Support Agency 
2 Figures in brackets 


indicate the literature references at the end of this paper. 


This adaptation applies equally 


given for the rectangular source. 


well to the corresponding expansion 


Formulas and numerical results 


For small values of ut, where yp is the attenuation 
coefficient and t is the barrier thickness, the Legendre 
series solution does not converge adequately for the 
particular angular distribution for the unscattered 
radiation from a plane isotropic source ° 

9(0)=9°(0) =(0/4) exp (—pt/cos 6)/cos 6, (1) 
in which o/4x is the source strength in units of 
radiation (e.g., photons) emitted per unit time from 
unit surface area into unit solid angle. The reason 
for this poor convergence is that, as pt—-0, g°(@)— 
(¢/4)/cos @, and a function of the form 1/cos @ {is 
not representable by the positive powers of cos 6 
appearing in the Legendre polynomials. An alterna- 
tive series expansion in powers of ut with geometry 
dependent coefficients, gn, was provided for calcu- 
lating the unscattered radiation component when 
O<utsSl. 

The present work extends the two expansion 
methods in reference [3] to include the case of a 
circular disk source with detector off-axis. A typical 
example of this situation would be the estimation 
of the neutron flux around the face of a cylindrical 
thermal column of a research reactor [7]. In addition, 
it will be shown that if the scattering properties of 
the medium are given as a point isotropic buildup 
factor [6] in polynomial form [8, 9] 


N 
B(ur)~ >> Bi (ur), 


i=0 


(2) 


rather than as angular harmonics, g;,, at the exit 
surface of the barrier, the g, expansion coefficients 
here and in reference [3] can often be used to evaluate 
the scattered as well as unscattered contribution to 
the radiation flux. 


3 Here and in what follows, g° and J° will be used to denote the angular distribu- 
tion and intensity of direct radiation, i.e., radiation which arrives at the detector 
from the source after having undergone no scattering interactions with the in- 
tervening medium. 
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2. Arbitrary Angular Distribution; Legendre 
Expansion of the Angular Response Function 


2.1. Isotropic Detector (Flux) 


If the emission from each source element, dS, 
of the disk, is independent of azimuth and uniform 
over the source surface S, the radiation intensity at 
any point in space, as measured by a small spherically 
symmetric detector, may be written as the integral 


J Sg(9) -vs(0)d(cos 6). f< 


Here @ is the polar angle as indicated in figure 1. 
The angular distribution function g(@) of radiation 
crossing the source-plane will now be defined as the 
radiation flux * (e.g., photons) traveling in direction 
@ per unit solid angle, unit time, and unit area normal 
to direction 6. wWs(@), called the 
function” in reference [4], is the azimuthal angle 
subtended by the source area from the detector-to- 
source-plane vertical for a given @-cone around this 
vertical. 

There are a variety of ways in which the integral 
(3) may be evaluated, depending on the nature of 
g(@) and S. The Legendre expansion method, es- 
pecially advantageous if g,; coefficients happen to be 
already available, has been discussed in detail in 
references [3] and [4], but for convenience a_ brief 
summary is given as follows: 

The two functions, g(@) and ¥s(@), in (3) which are 
parametrically independent except in 6, may (if free 
of singularities) be expanded in Legendre harmonics 


as 


“angular response 


q(é 3 (l+-4)9,-P (cos 0 (4) 
and 
Ws(0 pe + )p, -Pi(cos 8). (5) 


‘In this notation the liation flux in direction @ per unit solid angle, unit 


time and wnit area of the rce would be g(@)cos @, i.e., g(@) may be identified with 
the angular distribution functions ¢(2,9) and @(m) of references [4] and [7], re 
spectively, but differs from G(F) of reference [2] according to 9(#) =G(P)/cos 6 
< 
< 
; 
«. 
a 
- 
s > 
a 
> 
> ‘ 
ee 
FIGURE 1. Disk source and off-axis detector, showina polar 
angle 6 and “angular response function,’ Wg (0). 
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When this is done, the integral 


3) may be written 
as the series 


which converges rapidly for many cases of interest. 
The p; coefficients may be evaluated in the usual 
way according to 


p P, (cos 0) - ¥s(0)d(cos 6), (7 
or since ¥s(@)d(cos @) is a differential solid angle 
dQs(0), 

Pr P;, cos 8)dQ.(4) (S) 


As in reference [3], the Ws(@) and related prs ure 
dependent on the relative, rather than on the abso- 
lute dimensions, so that we may select one dimension, 
here the disk radius, as the unit of length. The 
quantities p and /, as indicated in figure 2, are then 
the distances off axis and off disk-plane, respectively, 
measured in disk radii. 

Since the poly nomials P (cos 7) are linear combina- 
tions of even or odd powers of (cos @) it is convenient 
to break up (dS) into components I (p,h ), each over a 
particular power of cos @, i.e. 


I ph COS 6d2(6,p,h), (9 
from which we reconstruct. the desired functions 
Pilpyh according to 

Pi ph) > a j;F (p,h), (10) 


where ay are the coefficients of x 
polynomial 


in the Legendre 


P I) Ja a’ 
of order lL. 
The differential solid angle, dQ(6,p,h), can be 


written in cylindrical coordinates h, R, and ¢@ around 





A 
r ¢ 
pi. 
o ? 
<< 
‘es . 
ANE OF 
a ~ 
> 
s 
FIGURE 2 Source-detector geometry, showing relevant param- 
elers p, h, 0, o, and R, in which linear dimensions are 
measured in radii of the source disk. 


the perpendicular from detector to source-plane as 


d2(0,p,h)=cos 6 dS/r?=(h? +R?) "AR dR dod, (11) 


h yh i, dS R dR dd, and , h? T Re’. 
The integral in (9) may now be written explicitly as 
the double integral 


since Cos @ 


F'i(p,h) 


h | 


R?) 2h d Rdo, 


(12) 


in which the limits hWiO.p,h) of the azimuthal 
variable @ are indicated in figures 1 and 2. The 
angular response function ¥(@,p,h) is the quantity 
¥<(0) defined following eq (3), in which the general- 
ized source detector geometry svmbol S has been 
replaced by the parameters p and A which character- 
ize the disk source situation. 


Now, since 


W(0,p.h 2 cos p?+R2—1)/2pR } . 13) 
we may write (12) as the single integral 
I (ph h ) 2RdR(h? I?) 
JR 1 
COS p” R? l } 2oh ° | }) 
Integration by parts and use of the substitution 
R2—= p? l 2p cos t (15) 


transforms (14) to the form 


h 
F'i(p,h 
jv 
> : dt 
p*—1) | 
p 1—2p cos t)(h?+-p?+1—2p cos t) 
|" dt ‘ 
. a b (16) 
(h? T p" T l 2p cos ft) . 


J 


The first of the two integrals in (16) may be further 
reduced by use of the identity 


l 
: ¢ fr 4 


wh? “=a h2/-2" (h? +a?) 


> 9 9 ’ Ga 
a?(h? +-a?) ri 


in which a?=p?+1—2p cos t. 

Since at this point the evaluation of the integral 
for the case of odd 7 (and thus odd /) is quite different 
from that for even 7, these cases will be treated 
separately in subsections A and B, as follows: 


rhe limits of integration p—1, p+1 in (12) refer to the case 1<p< For 
0<p<1 the integral should be split into the ranges 0< R<1—p and l—p< R<1+ p 
However, the resulting expressions are identical] 
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A, Even j,l 


Combining (16) and (17) we have for the F- 
integrals for even ); denoted as F;, where q Q, 
Ls 

h? 
I (p,h) : 
20 


ia | a ~ dt 


. p~ l 2p cost)y h? p" +-] 2p cos f 
I l [" dt 
) 
=p h?/-2* Jy (h?+-p?+1—2p cos ft) a 
(" dt ; 
. { Ss“) 
0 (h? p” l 2p cos t) ; 


The first integral in (18) can be evaluated in terms 
of complete and incomplete elliptic integrals of the 


first and second kinds. However, in combining 
these F’,,’s according to (10) to give p.;, it can be 
shown that for />0 this integral cancels out, 1.e., 


the polynomial > \a,,(27+-1) has the properties 
1 for /=0 

(19) 
0 for />0,/ even 


The zeroth coefficient Po( ph (which is the same 


as Fy(p,h is the response of an isotropic detector, 
at an off-axis position p,h, to radiation from a 
“perfectly diffuse” circular disk source, i.e., g(@) 

constant Fortunately, Po(p,h) is also the solid 


angle subtended by a circular disk from the point 
ph, and in the latter role has been discussed and 
evaluated by a number of authors [10, 11]. If we 
use one of these tabulations * for po, we need only 
evaluate the integral appearing in (18) in the sum- 
mation and also in the final term. Factoring out 


the quantity (A?+-p?+1)~“*!” from this integral 
reduces it to the form 
S, (kk) (1—k? cos t)~ +! dt, (20) 


in which n corresponds to the 7 and 7 in (18) and 
k?=2p/(h?+-p?+1)<1. This family of integrals may 
be evaluated in terms of complete elliptic integrals 
of the first and second kinds [12]. 


2 , A ied 
vi+k . (, za) 


Si(k) qd Byli | J.2 E(y 7B) 


and for n> 2, 


l 2n—-3 . 2n—3 4 
l p(t i) Se Qn [> }. 


(23) 


S,(k) 


S,,(k) 


For tabulations of po see, for example, references [10] and [11 
eters p and h in this paper were explicitly selected in order to use 
extensive table in reference [11] 


The param- 
directly the 








Extensive tabulations [13] of K and £ exist in addi- 
tion to an excellent polynomial approximation by 
Hastings [14] suitable for use on a_ high-speed 


computer. 
Substituting (20) into (18) we have 
F,,(p,h) cancellation terms in p2,(p,h) for 1>0 


l h? a SP 
254 Ai -p*-+-} ) Sk) 


aw 
eS 2 is 


ae ges 71) ‘sul }- (24 


The even-/ response function coefficients p2;(p,h) may 
then be built up from the F;(p,h)’s in (24) as indi- 
cated in (10). Equation (10) may be written more 
explicitly for even / as 


: l 21—7\ (4l—2) 
) >: = os 1 ( ) 
Poi(p,h — > ( j )( 9] j ) F; 2(p,h), 
(25 
in which the quantity 
! 
i tes 
oe 2%) 
r riin—r)! 
is the rth coefficient of an nth degree binomial 
expansion. 


B. Odd j,l 


Denoting the odd-j F-integrals in (16) as F, 
where 7=9, 1, 2, , we have, using (17), 


l _ dt 
F g+-1 wh) 5 (p? 1) 
z i( p,ht 2) 9 p a (p? } ] 2p cos 1) 
PUTS aoe” dt 
h2 |= . (h?4 p- +] 2p cos t)'* 


— nes | : . (27) 
. (h?+ p?+-1—2p cos ¢t)/*! ie 


the integrand denominators have integers, 
rather than half-integers, as exponents, closed-form 
algebraic expressions result, rather than quantities 
in terms of tabulated higher functions or infinite 
series expansions. Use of standard integral tables 
[15] to evaluate (27) gives 


Fy 541(p,h) 4 {, 
sided 2) +-2 L 
—(- | h? 7 i: ( h?+-p?+1 ) 
~ (h? p” -1)*+ 47h? \ (h2 p” +-1)?-+ 47h? 
_ : l —s ( h? ) ; 
h? = \ (h? p*+1)?-4 4p*h? 


<P, h'+e" } (2S) 
\ (h2 p” +1)?-+ =a) 7 


in which (cos 6)” is replaced by 


(ph 


Since 


{+96 - = ~ = —s)m 
\(h? -p°+) y (h? -p*+1)?-4 4 ph? 
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in the usual Legendre polynomials 7,(cos 6) to give 
the quantities ?; and ?; in (28). 

Substituting (28) into (10), we have for the first 
two odd-/ coefficients 


p (p.h 5( ] h re ) (29) 
- \ h? p” L)? T 4p°h? 
and 
j 7 f 6h? p~ l 
jp »./t) l 
Bs\t S 4 \ h p-+1 24 17h? 
5h? hh? p 


+— 1) (h? +p? +1) 
° (30) 
((h? p" +-1)?-+ 1 p*h?| . 


The higher odd-/ terms can be written in a form com- 
parable to (25) as 


Prsi(e,h)=2) “Sp 


21—j+1\ (4l/—2j+2\ ,, 
( ; ey, wy) Far-ay+1(0,h) (31) 


using the F; p.h)’s from (28). 

The p,(p,h) coefficient in (29) is equivalent to the 
response of a plane detector parallel to the plane of a 
perfectly diffuse (Lambert’s law) uniform. circular 
disk source, and has been derived and discussed by a 
number of authors {16}. 

The F3;(p,h) and F;4,(p,h 


28), and the expressions in 


expressions in (24) and 
(25) and (31) for their 
respective Legendre combinations were programed 
for evaluation of the pi(p,d)’s on an IBM 704 
computer. The resulting numerical values for the 
ppd) are presented in table 1 for 0<p<10, 0.1 
<h<10 and 1</<13. The entries for Po(p,h) were 
taken from reference {11 For p or hk greater than 
10, the source mia be treated for all practical pur- 
poses as a point anisotropic source with the same 
angular distribution, g(@) cos @, as each surface element 
of the source. The Pilp,h) coefficients in table 1 
may now be combined with g, coefficients as in 
equation (6) to give the radiation flux intercepted by 
a small spherical detector at coordinates p,h. 

As a check on these pilp,h re. they may be inserted 
in (5) to reconstruct the angular response function 
W(0,p,h Three detector positions (o=05, h==6:5- 
p=1, h=1; p=2,h=—1) are indicated in figure 3, and 
the exact angular response functions for these positions 
are given by the solid lines in figure 4. Summing 
over the first 10 terms in (5) gives us the approxi- 
mations shown by the dashed lines in figure 4. 


2.2. Cosine Detector 


As has been discussed in reference [3], a similar 
analysis of a finite uniform plane source geometry 
may be performed in the case of a surface-element 
detector with response proportional to the cosine 
of the angle of incidence of the radiation. For the 
case of a surface-element detector parallel to the 


TABLE 1.* Coefficients pi(p,h) of the Legendre expunsion (5) of the angular response function y(0,p,h 




















Pilp, h 
h=0.1 l 0 0.2 0.5 OS 1.0 
0 0 5. 6389 (0 5. 092 (0 4.9511 (0 2. 7036 (0 
l 0 3.1079 (0 3. (0 2.9346 (0 1. 4924 (0 
2 l | $ l ». 9042 (—1 1. 6880 (—1 
3 l l 7 (—1 1185 (—1 3. 2419 (—1 
1 l 1 2 (—1 3. 3127 l 9.7153 (—2 
3. 1 l 2 1 1. 2095 1 1.4395 (—1 
( l l ( l 2 l 1. 7067 ] 6.9956 (—2 
7 h 1 | l l 1. 8423 2 7.7337 (—2 
s l l | l 1 7. 5020 2 ». 2967 (—2 
v 1. OO68 l 2 6 2 5. S812 3 4.6161 (—2 
10 1. 2734 1 l ] l $1191 (—2 4. 2324 2 
11 $333 2 2 l 2 1. 7343 2 2. 8446 2 
12 1. 0734 1 1 v 2 2.2105 (—2 3. 4397 2 
13 2.2831 (—2 2 6 5 2. 4535 4 1. 7852 2 
1.2 2 ( 10.0 
6. O61 1 1. 7518 1 S816 2 6296 $ 3. 1768 4 
1. 3539 l 1.9422 2 $. 4561 $ ». 4493 ) 3. 2047 t 
2 2. 482 l 8. 3056 2 2. 6553 2 1. 3131 3 1. 5879 4 
; 1. 7525 l 2. 8306 2 5.1412 3 &. 1680 5 4. 8063 6 
4 1.0623 (—1 5. OS15 2 l 2 9.8181 (—4 1.1901 (—4 
1. 6643 l 3. 3581 2 6 3 4) 6.0060 (—6 
t 1.4981 2 3. 9768 2 1 2 $ 9. 9063 5 
7 1. 2740 l 3.6289 (—2 7 5 4 7.0039 (—6 
s 3. 5553 2 2.6539 4 1 2 4 8.6547 | 5 
; 7. 8309 ) 3. H85S 2 7 3 4 7.8748 6 
{ 1249 2 1. 5553 2 1 2 4 7.7741 
3. 6546 2 3. 5609 2 s 3 4 8.6560 (—6 
2 4. S885 4 6. 3576 $ s 3 4 7. 1095 5 
l 1. 1524 4 3. 2871 2 s $ (—4 9. 3693 6 
h=02 0 0.2 ( 0.8 1.0 
5.0510 (0 5.0156 (0 4.7828 (0 3.9593 (0 2.4052 (0 
1 , OBO (0 3.0113 (0 9394 (( 2.5521 (0 1.4145 (0 
2 Y242 l 6. 0653 l 6.9042 l S. 1628 1 2.6714 | 
6. O00 l ». 725 l 0554 l 1.7057 (—1 . HYO (—1 
i 4.0444 l 4. 0939 1 4. 2399 l 2. 4971 1 1. 4043 ] 
1. S601 | 1.7310 1 Y. 0334 2 1080 (—3 (—2 
( 2. SS45 l 2. SO48 ] 2. 5204 ] . 7523 (—2 2 
2. 4333 1. 3696 4. 0103 3. 7853 2 2 
S 2. 24 l 1.9326 l 1. 3335 l 5.9687 (—2 2 
iF) 1.972 2 ». 6315 2 7. 1927 2 2.6729 (—2 (—2 
10 1. 2879 l 1. IS95 | 5. 438 2 7. 2910 2 2 
ll S. 1495 2 8. 3167 2 6. 2859 2 1.8096 (—2 3 
12 7. 0789 2 6. O437 2 1. 2489 2 4. 1686 2 2 
13 Ss. SAUZS y4 §. 486 » 4, OUU0 2 2 SUIS 2 $ 
1.2 15 2 () 5.0 10.0 
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2. 4386 l 1. 3978 1 0043 4 2. 6108 3 3 4 
; 3. 4860 l 9. 4767 2 1.9485 2 3. 2513 4 l 5 
{ 5. H872 2 6. 59. 3. 2928 2 1.9343 3 2 (—4 
1.7291 1 2 2, 2940 4 4. (426 4 2 5 
: 1 2953 1 9 2 1798 2 1. 58090 3 1 4 
7 3. LISS 2 2 2 44% 4 4. 6802 4 2 5 
s 6. 7218 2 2 1. 2973 2 1. 3466 (—3 l 4 
” 2.8126 $ 2 2. 4463 2 5. 2123 4 3 5 
10 2. 3823 2 3 2 ». AG21 $ 1.1705 (—3 1. 5 4 
11 1. 2139 3 2 2. 3072 2 5.6621 (—4 3 5 
12 2. HH54 2 4.1470 4 5. 4447 4 1. 0276 3 1 4 
13 4. USH4 $ 1. 0522 2 ? 0585 2 6. 0426 $ 3 5 
h 0 0 0.2 0 OS 1.0 
0 $ 0 3.4184 (0 ( 2.4123 (0 l 0 
l 2 0 » 4802 (0 0 1.7426 (0 l 0 
2 ] ) 1.12382 (0 0 8. 2414 ] 4 l 
) 0 2. 3036, B | 1. 5252 l 7 2 
i ' 1 2304 l 1 6.8444 (—2 l 1 
3 l 2.0631 l l 6. 2302 $ 2 2 
f 4 2 2. 5241 4 2 9. 8306 2 7 2 
7 2 1 1.9522 1 2 7. 1346 2 4 
s ] 1 1. 2808 1 2 1. A869 2 2 2 
” 2 3. 2603 2 3] ». 7460 2 3 9 
10 1 1 1. 1162 1 ; » &509 2) 1 4 
ll 6 | 6. 26. 2 &. 6621 3 2 
12 5 2 3. 055 2 2 QI85 3 1. 2729 2 
13 Q 2 6. 8208 2 | ) 0240 2 1. 4838 2 
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0 l 2. 3242 l | 2 1. SNE ; 
] 1 6.04 » $330 
2 l 8. 3311 2 6). 2644 4 
3 l S. H034 2 wid i 
1 2 SSH7 2 1.34 { 
5 710 2 2 { 1.4771 } 
6 2.671 tid 4. 6834 { 
7 4.431 2 t ’ 1. 7041 ' 
8 4. 5078 9 99 O3t { 
4 2 1501 2 2.779 1. &kOD j 
10 4 4. 2228 2 is } 3618 j 
ll 2 162 ) » ®1Y 2» OA { 
12 i 2. 48 2 8. W073 \ 2. 8842 1 
13 3 2. 4048 2 } » 1626 { 
1.0 ) ( 8 1. 
0 1 0 l 0 6371 ( | 1S ({ 122 
l l 0 l 0 T6010 1.0921 s 1 
2 l 0 l 0 ¥. 4773 l 6.90 } ] 
; 5 F l 1.9529 l 03 1 1 
4 1 l ] ] 1. 4522 l 40s 2 2 
5 l l l l 4. 2344 2 SOAS s 2 
t 2 l 2 l 8.4731 1423 1 2 
7 1 l l l 4. SOHS2 2 4410 y. ] é 
Ss 5 2 5 2 4. 4288 3 26S84 4 
9 6 2 $ 2 G, 2186 3 4. TARO 1 
10 l l Ss 4 6. 4612 } . Howe l y. 
Bi 9 2 7 2 2. 7229 2 1, O68 >. 2 
12 3.3 2 2 2 3. 1179 2 sHitid 
13 3 2 l 2 1. 4585 2 i l 2 
l 4 { If ) 
0 l ( l , ney ] 9 9 ) 
l l } l 658 l f 107 j 
2 s 9 § 1552 9 ~ | 
; 2 1 1 ! 6. 97¢ 1.6313 { 
! l l 1 s 2 ( l “4 ; 
| | 1 1 2 6054 H12 i 
3 ; 7 2 1. 8398 tis4 1 
7 ; 2 2 7 r is { 
s 7. 7352 } 2 f } 1.1179 lye } 
Y 2.1418 2 $ $ } y 10S 6. 6240 { 
iT 4. 8204 2 7 3 5) s 4 Om { 
1] 2. 6571 2 ( 3 3 4. OOK 6. 7448 i 
12 2.4071 1 2 2 4 uu 175 4 
13 1.021 2 l 2 l 2 1.9591 6. 6418 i 
p 
h=2.0 l 0 0.2 0 0.8 1.0 
0 6 1 6. 5663 | 6, 2286 | 659 OM4 l 
l 6 l l RDF 1 1359 1. (OOS 
2 5 l l l 4. 200; 3. 55US l 
} 4 1 l l , 0454 2 3208 1 
J 3 l l l 1.89 1. 1540 | 
5 2 l l l 4, 28k 2. 7138 2 
t l l l 2 2.7404 2 Hs2 2 
5 2 2 2 4. 933 3 O35 
& 1 ‘ 2 ; } 1. { 2 404 - 
4 6 2 2 2 1. 1624 4.4048 
10 s 2 2 2 682] 1, 102 9 
11 s 2 6. 8555 2 3 7O00 y. 1.8444 2 
12 7 2 SIRS 2 $ 022 1.4479 
13 4.6553 > 4 4.4432 2 $ 4418 4.7874 } 
1.2 l 0 10.0 
0 4. 6970 l $. 9438 l 2.8271 1 4. 1466 2 847 
l 4.0376 l $9134 l OU 1 ] 7 2 1. ISS2 
2 2.9065 l > 10014 l &. GOR4 2 1. 1646 » 640 3 
1.6107 | 6. 9906, 2 2 17h 2 1. 7831 1. 6583 
} 4. 6708 2 2 s MS 2 3. 1704 1. 405¢ 
3.0151 2 ) &. O892 2 } 1. S25¢ ; 
6 6. 2584 2 2 0978 2 1. O89 5. OSSO j 
7 2 » 190 2. 4108 ; 1. 7496 3 
Ss 2 2 01 2 1. 1062 2 2. 0353 i 
4 2 3 6461 2 689 1. 4767 3 
10 3 2.1731 2 909 7. 3707 4 
11 $ ; 1477 3 9.1154 1.0619 3 
12 3 3 6.8639 3 1.5814 1. 0761 ; 
13 - 3 GOR4 oor 6923 { 
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Pile 
0 0 0 
1. 2202 l l 
1. 2083 l ] 
1. 1848 i l 
1. 1502 1 l 
} 1. 1051 | l 
5 1. 0503 1 l 
‘ YU STOO 2 2 
4. 1627 2 2 
s S. 3046 2 2 
4 7. 5798 2 2 
( 6.7330 (—2 2 
SOUu0 4 2 
y 4 » OFS v4 2 
4. 1487 2 3.877 2 
1.2 ] 
1277 1 OSU l 
SAU l 1.0201 1 
2 O140 l 4 2074 2 
; 1 0832 y. 9469 2 
4 933 2 6. 3284 2 
(). 3572 2 4. 5910 y 4 
‘ ‘ 2 2.8750 2 
7 w 4 1. S078 4 
s 2 8.7710 
; ] 1. OOM 
1 1. 6621 2 
] 1. USI 2 
2 2 ? W9S 2 
2 1. 812 
oo ’ ov 
3. LINZ 2 3. 1164 2 
l I1f 2 3. LONI 2 
OSI 2 3.0914 2 
$. OF 20 ) 3%. 0666 2 
j s O414 2 3. 0337 2 
0035 2 2.9930 2 
( 2. O584 2 2 9446 2 
2 GO4 2 2. SSS 2 
s 2 8476 9 ? $261 2 
Q 2. 7R2 2 2. 7A 2 
10 2.7113 2 ? HSOS 2 
11 2. 6343 4 ). SOY 2 
12 2. 5520 2 2. 5119 2 
W47 2 2. 4198 2 
0 ; 2 3. O176 2 
3 2 2. O77¢ 2 
> O77 g 2 ROS 2 
SASS 2 » TR32 2 
t 74 2 2. 6336 9 
422 2 2. 4540 2 
» 4800 ? 2. 2487 2 
2. 3181 2 2. 0229 2 
S 2. 1327 2 1.7821 2 
9 1. 9360 l 18 2 
0 1. 7318 2 1. 2780 4 
1. 523 2 1. 0261 
] 1. 3148 7.818 3 
13 1. 1089 ) 5 4916 3 
I f parenthese t , 2, and 3 indicate the power of 10 by wl 


source-plane, table 1 may be used in two different 
Wiys: 

a) If the angular distribution data is already in the 
form of g, coefficients of a Legendre expansion as 
given in (4), a new set of coefficients p:(p,h) may be 
obtained from the entries in table 1 according 
Lo 


p,h ) 


which can be combined with the g; data as in (6) 


Pi(pyh lp, (L+-1) pisi(p,h) } /(214-1) (32) 


wh 








>) of the angular re sponse function ¥ 6,p,h Continued 
oS 1o 
1, 2033 l 1.1776 l 1. 1547 
1. 1863 l 1. 1520 l 1. 1235 | 
1. 1528 l 1. 1050 l 1. 0628 
1. 1041 l 1 0358 1 4. 7041 
1.0417 l 9. 4882 8. H923 2 
9. 6759 &. 4810 2 1741 2 
8.8421 2 7. 3820 2 6.1761 2 
7.9414 2 6, 2393 2 4. 8664 2 
7. 0609 2 5. 1001 2 §. GOR 2 
6. 0477 2 4. OOR6 2 2. 4596 
L080 2 3. 33 (—2 1. 4624 2 
4. 2057 2 2.1148 2 6. 4728 5 
4. SO1Y 1. 3644 2 2. 92. } 
939 2 6574 $ 3. Y175 ; 
0 10.0 
Y. 88903 2 $. 4590 2 1. 1307 2 
9. 1512 d 3.1724 N32 ; 
7719 2 1. 1644 2 2. 220 
G208 2 7. 1461 3 0443 $ 
$. S554 2 1. 7407 2 2. 3029 5 
1. 7058 2 1. 647 2 2. 1603 $ 
2. 3434 4 7.113 3 6874 $ 
1. 42909 ) 3 l OD 
2. 3258 2 2 2.1739 3 
2. 8458 2 2 3 
2. 6000 2 3 4 
2. 1555 2 3 » 
1. 5101 2 3 2 5 
®. 1402 3 i 
0 0.8 1.0 
3 1068 ? 3. OR90 2 3. OF 28 2 
OY 2 3. O718 2 3. 0504 2 
4. O72 2 3. O3876 2 , O50 y 
4. 04384 2 2 USE 2 4400 2 
OO3t 2 2.9202 2 2. 8539 4 
Y382 2 2.8385 (—2 2. 7488 2 
2.8728 2 2. 7429 2 2. 6266 2 
2. 7981 2 2.6344 2 2. 4880 2 
2.7145 2 2. 5145 2 2. 3381 2 
2. 6228 ) 2. 3848 2 2. 1764 2 
2. 5239 2 2. 2466 2 0063 2 
2.4184 2 2.1018 2 1. 830 2 
2. 3074 2 1.9520 2 1. 6506 
2. 1916 2 1. 7990 2 1. 4702 2 
o4 0 10.0 
2.9429 2 2. 2412 4 1.1118 2 
&RO] 2 0042 2 7. 8735 3 
7508 2 1. SASS 2 3 
1. 0051 2 3 
? 4. 0196 3 $ 
2 1.4972 3 3 
2 5. 7256 3 3 
?. §. LAU 3 3 
2 8. 6369 3 3 
$ 7. 3437 3 3 
3 4.758 $ 3 
4 1. 4090 3 4 
} 1. 6059 3 3 
$ 4. 0706 3 3 
h the adjacent entry is to be multiplied; e.g., 2.2831 2) =0.022831. 


If the source data are given in terms of the 
actual angular distribution g(@), which could be 
either a mathematical expression or experimental 
data, it would be simpler to modify the source 
angular distribution coefficients, that is, to generate 
Legendre coefficients g; of the quantity 
cos 6q(@). These g; coefficients could then be com- 
bined directly with p, values from table 1 in the 
same way as has been specified in (6) for an iso- 
tropic detector. 


q' (0)= 
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DISK SOURCE ,S 

Figure 3. Three examples of detector positions in which the 

distance 8 off-axis and off source-plane, measured in 

radii, are (A) p=0.5, h=0.5, (B) p=1, h=1 
h=1. 


























Figure 4. Reconstruction of the response da 9g oH 

according to (5) for (A) p=.5, h=.5, (B) p h=1 and (C 
p=2, h=1 using table 1 aa summing over 7 basen 0<1<9. 
The reconstructed functions are shown as dotted lines and the 
exact functions as solid lines. Note that the scale for (B) 
is on the right-hand side. 


y (0,p,h 


2.3. Limitation of the Legendre Expansion Method 


The utility of this method is dependent on the 
oscillatory nature of P,(cosé@). Within any finite 
region between cos 6 and cos @+ Acos @ the increasing 
number of these oscillations, for increasing J, yields 
decreasing values of p; and g; as a result of higher 
and higher order cancellations. If both functions 
¥s(0) and g(@) are concentrated mostly within very 
narrow regions of cos 6, adequate convergence of (6) 
may not be achieved without going to exorbitantly 
high /’s. 

‘Evaluation of the unscattered ieidiatinn com- 
ponent provides an example of the above limita- 
tion. This problem has been qualitatively discussed 
in section 1, and will be treated in some detail 
the following section. 


disk 
and (C) p=2, 


3. Unscattered Radiation Component 


3.1. r>1: Legendre Expansion 


As in the rectangular source case [3], the Legendre 
sum (6) converges adequately for an angular dis- 
tribution of the form (1) when the barrier thickness 
ryt is of the order of one or greater. In this 
region the behavior of g(@) at the barrier exit surface 
lies between that of a “perfectly diffuse” plane 
source (g(@)=constant! and that of a plane normal 
source {g(@) od(@—0) The g? coefficients are 
easily obtained from appropriate combinations of 
tabulated expone ntial integrals [17], and a sample 
tabulation of g?’s is given in reference [3] for 
0</<13 over the range 0.001 <7 < 100.0. 

3.2. O<r<¢1; h>O: Power Series Expansion 

Since one of the chief reasons for separate treat- 
ment of the unseattered radiation component is the 
inadequate convergence of (6) for +: l, an al- 
ternative means of evaluating the integral 


[°(r, p,h) (33) 


(dS /r?) exp (—7/cos @) 


is desirable. Again, as in [3], no success was had 
in obtaining a closed-form result for /°(7,p,h) in (33). 
However, Sievert [18] has shown that by expanding 


exp (—r/cos@) as a (7/cos 6)" power series, a solution 
corresponding to eq (10) in reference [8] may 
be written 
[°(r,p,h) >> In(p.h)r", (34) 
4 fa 


in which coefficients 


the qn(p,h) are evaluations 

of the integrals 
( 1)* (° yo\ 57h ( p? he? 
Yn(p,h) nh it +-F2?)2, 2 ¢os ( Dl? )RaR 
(35) 

for 1<p<o@ and 

[}* i cai 
Gn(p,h) To (h?+R?)? 2aehdR 

i J0 


+ (h?+-R?2) = 2cos 
l—p 


for O<p<l. As in the p,(p,h) evaluation, three 
distinct expressions arise from (385) and (36) de- 
pending on whether n is zero, even or odd, but by 


proper manipulation the distinctions between re- 
sults from (35) and (36) were removed. In reference 
[18], an(p,h) coefficients are given only for n=0, 1, 
and 2, of which ¢,(p,h) is not exact. 

The zeroth coefficient qo(p,h), when multiplied by 
the source strength factor o/4n, is the radiation 
field from a disk source in a nonattenuating medium 
(i.e., 7=0). This can be written, for h>0, as [18] 
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An/o—a In\(1+h?— p? 


ry (1 +f? 


Jo(p, h) i. 


o\9 912\ /91,9\ ~— 
ad ad 4p*h?) 2h2;. (37) 

Higher even-n terms (1.e., n(e)=2, 4, 6,...) are 
closed-form purely algebraic expressions, the next two 
of which are 


/ 7 3 
(2(p, lr) ON? (38) 
and 
T 
: 2 ae 2 
gs (ph) - (h?+ p?+3). (39) 
fs 1')3 , e 


These are included in the general expression for the 
even-n terms: 


, h) J l 
Vnce\(P a 


h?+-p?+-1 2 f tmaxgni4 27 /2 2 
(eS CCE) axa) 
2 l ‘ 4 1/2 2i\ . 2 h2 n/2 
te l = ( ; (Ge i) [Gas i) 


~ * 
as Os 
ition 


The odd-n 


’ ' ) are all 
included in the series-expansion sol 


terms (1.e., n(o)=1, 


q (p,h 


h? p* L\* ‘ 27 n/2 y . 
+( he ) {35 (7) (2?) axa) 


=  (ni2\ / 24] m 
2 | I) (Vain) |} } (41) 


which strongly resembles (40), except that here the 
quantities (’¥7) and (”/?) now coefficients of the 
infinite series expansion of half-integral powers of 
a binomial [19], ie., 


are 


. (ni2 
n/2 ” 9 
(1+2) ("ha (42) 
m 0 m 
These are readily generated using the recursion 
n/2 > 
(er 
0 
) 
LTP 4 
( ) n/2, and 
l 
n/2 2(m—1)—n/ n/2 = 
( ) ( ); (43) 
m 2m m—1/ J 





which is also valid when n/2 is an integer, 
using (43), is zero for m>n/2. 
By setting p 


since {"/?) | 
0 in (37), (40) and (41) and substi- 
tuting the resulting values for q,(p,h) in (34), we 
find the on-axis unscattered flux to be 


[°(0,h) =(0/4r) 2x [in vhi+ 
h 


1)” h?+1\ 
> Oe \ )- r” ) 
Toa nl { ( a i} | ae 


This series can be identified with the familiar on- 
axis solution in terms of exponential integrals [20] 


yre+1e7 . 
2ardr 


h 


T°(0. A) 


(o/47) 


9 





Ey (rvh?-+1/h)}. (45) 


(o/4m) 2m} Ey (7) 

Equations (37), (40) and (41) were programed 

for evaluation on an IBM 704 computer and a 

resulting set of g,(p,h) coefficients for 0<p<10, 
0.1<h<10 and 0<n <9 is given in table 2. 


3.3. Comparison of Legendre and Power Series 
Expansions 


As a partial check on tables 1 and 2, and because 
the unscattered radiation component is of some gen- 
eral interest, the quantities J°(7,p,h)42/o as defined 
in (33) were evaluated using the methods of the 
preceding two sections, and are given in table 3. 
Numbers preceded by asterisks were evaluated using 
the Legendre expansion approach (6) with g{ values 
from reference [8] and 0</<13; the others are 
evaluated using the 7" series (34) for 0<n<20. 

Some notion of the convergence for each method 
may be obtained from table 4, in which partial sums 
over 0</<13 and 0<n<13 for 7=1.0 are given for 
the three situations shown in figure 3. A qualitative 
comment about the Legendre sums for these three 
cases is that in each case the first two terms con- 
tribute over 80 percent of the result. 

For scattered radiation the tendency is toward 
even faster convergence. For example, the Fermi 
approximation [21] 


g(0)= (14+ v3 cos 6)/(14+ 73/2), (46) 


used in reference |7| for the angular distribution of 
neutrons from a reactor thermal-column face, has 
been shown by Placzek [22] to be in error by less 
than one percent for 0.1<.cos 6<1, with a maximum 
error of 7.2 percent at cos 6=0. The Legendre ex- 
pansion coefficients for the g(@) in (46) are: go= 
1.072, g,=0.6188 and g, 0, which may be used 
with the p; values in table 1 according to (6). 
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TABLE 2. Coefficients q'n(p,h) [=h"qn(p,h)| of the (r/h)™ power series expansion (34) of the unscattered component (4x/0)1°(r7,p,h 
of the radiation flux, from a circular plane isotropic source, inc ident on a detector at position p,h separated from the source- 
plane by an attenuating layer of thickness r=yt. Each term contains the factor h" in order to include meaningful tabulations 
an the region p>1,h~0O as in (47 


h=0 p 0 0.2 0.5 s 1.0 
n 
0 
1 
> 
3 
} No physical solution wit! the r 1 rl 




















SN 
y 
p 1.2 1.5 2.0 , 10.0 
0 3.6157 (0 1 9. 0378 l 1. 282 l 3.1574 2 
! 2.9533 (0 2 1.6252 (0 6.3151 l 3.1455 1 
2 1. 5708 (0 1 1. 5708 (0 TOS (0 1. 5708 (0 
$ 6. 8470 1 bad l 1.0803 (0 2 6311 (0 5 
4 2.5395 (—1 3. { 1 ». 8905 l 3. 3379 (f l 
5 8.1621 (—2 1.3 1 ? H8965 (—1 $ 4200 (0 2 
6 2. 3069 2 4.° 2 1. 0617 1 2.9467 (0 4 
7 5. 8004 3 1.2 2 3. 6999 2 2 1951 (0 6 
8 1. 3158 3 3.2 3 1.1551 2 1.4425 (0 s 
9 2.7095 4 2.8 4 3. 2677 3 &. 400) l 9 
h=0.1 p 0 0.2 0 s 1.0 
0 1.4499 (1 l 1.3619 (1 1.1476 (1 7. 3908 (0 
l 5. 6862 (0 5. 2804 (0 1.5455 (1 3.7113 (0 
2 1. 5708 (0 l 1. 5708 (0 1. STOR (O 1. 5708 (O 
3 3. 53897 (—1 3.6 1 4. 1803 l 1342 1 9576 l 
4 6. 6759 2 7 2 9. 9484 2 1. 5054 l 1. 9766 ! 
5 1. 0736 2 & 4 2, 1058 2 + 840 2 ». 740 4 
6 1. 4985 3 l 5 3. 9747 3 8. 9266 $ 1. 467¢ 2 
7 1.8441 4 2. 5 4 6. 7332 4 1. S28 ; s. S511 8] 
8 2.0270 (—5 3 1.0316 4 3. 3TRS 4 6. 8960 i 
9 2.0120 (—6 3.¢ 6 1. 4403 5 680. 1. 2003 4 
p 1.2 1.5 2.0 5.0 10.0 
0 3.6033 (0 l l 1. 2810 l 3. 1571 2 
1 4 (0 2 (0 6. 3138 1 3. 1454 l 
2 (0 1 0 1. 5708 (0 1. 5708 (0 
3 l 8.¢ 0 2.6316 (0 5. 2428 (0 
4 1 3 1 3. 3393 1.3157 (1 
5 - 2 1.é (—] $.4220 (0 2.6479 (1 
6 2 4.3 l 2.9489 (0 4.4516 G1 
7 3 1.2 2 2.1971 (0 6.4304 (1 
‘ 3 3 2 1.4440 (0 8.1470 G 
4 - 4 7.é 3 8. 5000 9.1963 (1 
h=0.2 p 0 0.2 0 0.8 1.0 
0 a 1 al 9.4210 (0 7.6461 (0 5. 3698 (0 
l 0 0 ) (0 41.7664 (0 4.1033 (0 4.4555 (0 
2 0 1.§ (0 1. 5708 (0 1 7O8 (0 1.5708 (0 
3 (—1 3 (-—1 4. 3050 1) 5. 241 l 6. 0468 1 
4 (—2 7.¢ (—2 1. 0341 l 1. 544 1 2.0159 (—1 
5 2 1. < (—2 2. 2013 2 4. O11¢ 2 ». RBRS 2 
f 3 2 (—3 4.1776 3 +, Wt $ 1. HOTS 2 
7 4 2 4 7.1178 4 1. 8901 ; $. 4547 ; 
8 (—5 3. 3 5 1.0971 (—4 $. 5247 4 7. 1350 4 
9 > (-6 3.7 6 1.5412 (—5 5. 9423 1. 3399 4 
p 1.2 1.5 20 5 | 10.0 
at 
0 3.2477 (0 1 l 1. 2803 1 2 
1 2.8313 (0 2 (0 j l 1 
2 1. 5708 (0 1 0 0 0 
3 6.9432 1 8 0 { 0 
4 2 1 3 | 0 | 
5 8.¢ 2 L.é (—1 { l 
6 2 2 4 (—] ( ] 
7 6 3 1.3 2 () l 
8 lL. 3 3. < 2 0 l 
9 2.8243 (—4 7.7 3 1 1 





TABLE 2. Coefficients q'» p,h) [=hqn(p,h)] of the (7 h)" power series expansion (34 of the unscattered component (4x a) I°(r7,p,h) 
of the adiation flus, from a circular plane isotropic source, incident on a detector at position p,h S¢ parated from the source- 
plane hy an attenuating layer of thickness r= ut. Each term contains the facvor hn in order to include meaningful tabulations 
in the region p 1,.h~0O as in (47 Continued. : 














h=0 r 0 { 0 0.5 1.0 
0 5.0562 (0 4.9753 (0 4.5353 (0 3.6855 (0 2.9550 (0 
1 4, RRS? (0 48383 (0 3, 6032 (0 $1826 (0 2. 8363 (0 
2 1. 5708 (0 1. 5708 (0 1 TOR (0 1. 57C8 (0 1. 5708 (0 
} 4.4420 1 4.5354 (—1 5. 0155 1 5. 8610 (—1 6. 5368 (—1 
4 Y. 8175 2 1.0341 (—1 1, 3090 l 1.8195 l 2.2907 (—1 
1.7967 (-—2 1.9734 2 2. 9378 2 4.8879 (—2 6. 6468 2 
t 2 SISO 3 3. 2613 3 ». SI78 3 1.1587 (—2 1.8090 (—2 
$. 8751 4 1.7714 4 1. 0318 3 2.4543 (—3 4.2497 (-—3 
s 4.7480 5 6.2749 (—5 1.6557 } 4.6921 (—4 8.9969 (—4 
4 ». 2476 6 7. 5023 6 2. 4234 (—5 8. 1650 | 1.7321 (—4 
p 1.2 1.5 2.0 5.0 10.0 
0 2550 (0 1 0 &. 2602 1. 2690 ( 1 3 ( | 
l 2.4916 (0 2 0 1. 5625 (0 6 25 ( l 3. ( l 
2 1. 5708 (0 1 (0 1. 5708 (0 1. 5708 (0 1 (0 
5 7. 4054 1 8.7 l 1.1135 (0 2.6442 (0 5.2 (0 
t 2.8667 (—1 3 l 6.2177 (-—1 3. 3707 (0 ] 6 
5 9.4997 (—2 l 1 2.8953 1 3.4695 (0 2 8 
6 2.7574 (—2 4.§ 2 1. 1626 l 3. 0026 (0 4 1 
7 7.1192 3 1 2 4.1153 4 2.2463 (0 6.4683 (1 
s 1.6543 3 3 3 1. 3040 2 1.4823 (0 8. 2044 (1 
9 3. 4929 4 9.1603 (—4 3. 7422 3 8.7601 (—1 9.2716 (1 
h=1.0 p 0 0.2 0 0.8 1.0 
0 2.1776 (0 2.1463 (0 1.9878 (0 1.7219 (0 l 
l 2. 6026 (0 2.5805 (0 2.4687 (0 2.2805 (0 2 
2 1. 5708 (0 1.5708 (0 1.5708 (0 1.5708 (0) 1 
; 6. 3824 l 6.4562 (—1 6. 8363 (—1 7.5100 (—1 8 ] 
4 1. 9635 1 2.0159 (—1 2. 2907 1 2.8013 (0 3 | 
5 4. 8767 2) 5.0998 (—2 6. 3027 2 8.6736 (—2 1 1 
6 1.0181 (—2 1.0886 (—2 1, 4817 2 2. 3138 2 3.1998 (—2 
7 1. 8368 (—3 2.0161 (—3 3. 0566 3 5.4471 (-—3 8. 2091 3 
s 2.9219 (—4 3.3053 (—4 5. 6368 i 1.1504 (—3 1.8895 (—3 
v 4+. 1608 5 4. 8702 9. 4190 2. 2055 4 3. 9469 4 
f 1.2 1.5 2.0 5.0 10.0 
0 1.290903 (0 1.0123 (0 6. 6582 1 1. 2302 ( 1 3.1257 ( 2 
| 1. 9698 (0 1. 7394 (0 1.4173 (0 6. 1877 1 3.1298 (—1 
2 1. 5708 (0 1. S708 (0 1. 5708 (0 1.5708 (0 1.5708 (0 
3 8.7743 (—1 9. 9202 l 1.2063 (0 2.6832 (0 5. 2687 (0 
1 3. 8485 ! 4. 9087 l 7.1995 (—1 3. 4688 (0) 1.3286 (1 
1. 4065 1 2.0193 (—1 3. 5482 l 3.6193 (0 » 6869 (1 
t 4. 4562 2 7. 1540 2 1.4981 (—1 3.1736 (0 4.5389 (1 
7 1.2351 (--2 2.2333 (—-2 5. 5516 2 2.4045 (0 6. 5876 (1 
s 3. OS30 3) 6.2418 (—3 1. 8369 2 1. 6064 (0 8.3854 (1 
i) 6.9802 (—4 1. 5806 (—3 5.4956 (—3 9.6077 (-1 9.5095 (1 
h=2.0 { 0 0.2 0 0.8 1.0 
0 7.0103 l 6. 9603 l 6. 7068 l 1 5. ¢ l 
1.4833 (0 1.4777 (0 1.4491 (0 l.é (0 
2 1. 5708 (0 1. 5708 (0 1. 5708 (0 J (0 
5 1.1102 ( 1.1148 (0 1.1391 (O 1 0 
t 5. 8905 1 5. 9428 1 6. 2177 l l 7 l 
5 2. 5030 1 2. 5382 1 2. 7251 l | 3 (—] 
t &. 8721 2 9.0473 2 9. 9902 2 1 l ( l 
7 )_ G9S83 2 2. 7684 (—2 3. 1511 2 3.9278 (—2 4.7 2 
S 7. 1878 3 7. 4236 $ 8. 7303 3 1.1475 (—2) ! 2 
’ 1. 7037 3 1.7722 (-3 2. 1583 3 2.9996 (—3 3 3 
r 1.2 l 2.0 5.0 10.0 
0 5.5179 1 4. 80905 1 3. 9168 l 1.0970 ( 1) 3. 0342 2 
l 1.3098 (0 1. 2327 (O 1.1012 (0 5.8485 (—1 3. OR39 1 
2 1.5708 (0 1 1.5708 (0 1.5708 (0 1.5708 (0 
; 1.2 0 1. 36 1.5157 (0 2.8335 (0 5. 3464 (0 
} 7 l & 1.1127 (0 3.8615 (0 1.3679 G 
5 3.8510 (—1 4 6. 6249 1 4.2401 (O 2.8063 (1 
6 1.6060 (—1 ? 3. 3307 l 3. 9066 (0 4.8085 CU 
7 5. 71 2 Ss 1. 4532 1 3. 1058 (0 7.0780 
s 1. 8501 ) 2 5. 6119 2.1742 (0 9.1364 (1 
9 §. 2954 (—3 8 1. 0471 2 1.3612 (0 1.0506 (2 





TABLE 2. Coefficients q’ n(p,h) [=h"qn(p,h)] of the (r/h)" power series expansion (34) of the unscattered component (4x/a) 1° (7,p,h) 

of the radiation flux, from a circular plane isotropic source, incident on a detector at position p,h se parated from the source- 

plane by an attenuating layer of thickness r=pt. Each term contains the factor h" in order to include meaningful tabulations 
in the regular p>1,h=0O as in (47)—Continued 


Qn’ (p,h hrqn(p,h 











h=5.0 m p 0 0.2 0.5 0.8 1.0 
n 
0 ] 1. 2206 l 1. 2031 l 1. 1873 ] 
] i 6. 1922 1 6.1471 ] 6. 1063 l 
2 1. 5708 (0 1. 5708 (0 1. 5708 (0 
3 2.6568 (0 2.6767 (0 2.6049 (0 
4 3.3707 (0) 3.4217 (0 3. 4688 (0 
5 3.4215 (0 3.5003 (0 3. 5735 (0 
6 2.8946 (0 2. US4t 0 3. 0689 (0 
7 2.0993 (0 2.1820 (0 2.2599 (0 
& 1. 3324 (0 1. 3962 (0 1.4568 (0 
9 1 7.5174 1 7. 9429 1 8.3502 (—1 
p 1.2 15 2.0 5.0 10.0 
n 
0 1.1686 (—1 1.1355 (—1 1.0699 (—1 2 2. 5193 2 
1 6.0574 (—1 5. 97 1 5.7943 (—1 1 2.8111 l 
2 1. 5708 (0 l 1.5708 (0 0 1.5708 (0 
3 2.7170 (0 —2 2.8422 (0 0 S610 (0 
4 3. 5264 (0 3 3.8615 (0 0 1.6428 (1 
5 3 3 4.2021 (0 ( 3.6895 (1 
6 3 3. 37 3.8150 (0 l 6.9160 (1 
7 2. 2. 54: 2.9723 (0 l 1.1130 (2 
® 1.5 1. 6802 (0 2.0286 (0 l 1. 5695 (2 
G 8 1 ¥. 8760 l 1. 2322 (0 0 1.9705 (2 
h=10.0 p 0 0.2 0.5 OS 0 
nv 
0 3. 1260 2 3. 1248 2 2 3. 1064 2 3. 0955 2 
l 3.1338 (—1 3. 133% ] 3 ] $. 1239 l 3. 1184 l 
2 1. 5708 (0 l (0 l 0 1. 5708 (0 
3 5. 2491 (0 5.2501 (0 5 (0 0 ». 2750 (0 
4 1.3155 (1 1.3161 (1 ] l l l 1.3286 (1 
5 2.6377 (1 2.6392 (1 2. l 2.6630 (1 2.6772 (1 
6 4.4071 (1 4.4106 (1) i l 4.4637 (1 4.4957 (1 
7 6.3116 (1 6.3180 (1 6.2 l 6.4133 (1 6.4710 (1 
8s 7.9093 (1 7.9188 (1 7.§ a 8.0629 (1 8.1501 (1 
9 8.8102 (1 8.8226 (1 s (1 9.0106 (1 9.1247 (1 
p 1.2 1.5 20 50 10.0 
n 
0 3. OR23 2 2 2 2. 072 2 2 
l 3.1117 (—1 1 l 2. 8060 ] l 
2 1.5708 (0 0 1. 5708 (0 0 
$ >. 2864 (0 0 5.8645 (0 0 
4 1.3344 (1 l 1.6428 (1 l 
i) 2.6947 (1 ] 3.6829 (1 ] 
6 4.5349 (1 l 6.8833 (1 2 
7 6.5419 (1 l 1. 1031 2 
8 8.2576 (1 l 1. 547 2 2 
y 9. 2656 (1 4 1. 930, 2 





TaBLE 3. Unscattered component (4x/0)I°(r,p,h) of the radiation flux calculated (a) using table 1 of this work and table 2 of 
reference [3] (numbers preceded by asterisks) and (b) using table 


of this work (no asterisks). 


An/o) I rT, p,h 
h=0.1 0 0.2 0.5 0.8 0 12 5 20 5.0 10.0 

0 1.45 (1 1.44 (1 1.36 01 115 (1 7.39 (0 3.60 (0 1.83 (0 0.00 l 12 (—} 3.16 
0.01 1.39 G1 1.38 (1 1.31 (1 1.10 (1 7.03 (0 3.33 (0 1.62 (0 7.53 l 7.84 2 1.17 2 
0. 02 1.34 (1 1.33 (1 1.26 (1 1.06 (1 6.71 (0 3.08 (0 1.44 (0 6. 31 1 4.81 (-—2 4.31 $ 
0.0 1.20 (1 1.19 G 113 01 9.54 (0 5.86 (0 2.47 (0 1.0 0 3. 7¢ l 1.13 2 
0.1 1.01 01 1.00 (1 9.57 (0 8.11 (0 4.81 (0 1.76 ( 5.99 l 1. 66 

"1.00 (1 Y.98 (0 "9.56 (0) *8.02 (0 *4.72 (0 *1.65 (0 
0.2 7.38 (0 7.34 (0) 7.10 @) 6.12 (0 3.45 (0 9.94 1 

*7.27 (0 "7.24 (0 *7.04 (0 *6.07 (0 *3.41 (0 *Y. 38 ] *2.09 *4.565 (—2 *1.36 (-—3 *1.34 (-—4 
0.5 *3.50 (0 *3.50 (0 *3.46 (0 (0 *2. 51 "2.74 2 *3. 45 $ *2.17 (-—5 "1.72 
1.0 *1.38 (0) *1. 38 (0 "1.38 (0 l *4. 25 2 "1.03 3 *2. 49 4 
2.0 *3.08 (—1 *3.08 (—1) *3.08 (—1 l *2. 11 $ 
5.0 *7.22 (—3 *7.22 (-—3 *7.21 (—3 3 


10.0 *2.61 (—5) | *2.61 (—5 *2.61 5 





TABLE 3 Unscattered component (42/0) 1° (7,p,h) of the radiation flux calculated (a) using table 1 of this work and table 2 of 
reference |3] (numbers preceded by asterisks) and (b) using table 2 of this work (no asterisks Continued 


Anja) I° (7, p,h 




















p 
h=0.2 0 0.2 0.5 OS 1.0 1.2 1.4 2.0 5.0 10.0 
0 1.02 (1 1.01 (1 9.42 (0 7.¢ (0 37 (0 3295 (0 1.77 (0 &. 90 l 1. 28 l 3.16 2 
0.01 9.98 (0 9.87 (0 9.19 (0 7.44 (0 5.20 (0 3.11 (0 1.67 (0 8.13 l 1.00 l 1.92 (—2 
0. 02 9.74 (0 9.62 (0 &.96 (0) 7.25 (0 5.04 (0 2.98 (0 1.57 (0 7.43 (-—1 7.83 2 1.17 (-—2 
0.05 9.04 (0 8.94 (0 &.32 (0 6.71 (0 1.60 (0 » 63 (0 1.31 (0 5.70 (—1 3. 76 2 1.63 ( 3 
0.1 8.01 (0 7.92 (0 7.38 (0 5.93 (0 3.97 (0 2.15 (0 9.78 (—1 3.70 1 1.12 2 1.01 (-—3 
0.2 6.35 (0 6.28 (0 5.88 (0 4.71 (0 3.17 (0 1.49 (0 5.66 (—1 1.62 ( l 
*6.39 (0 *6. 36 (0 *5 8S (0) *4.70 (0 *3.02 (0 *1.46 (0 *5. 36 ( 1 “1.38 l *3.48 3 "3. 36 4 
0.5 3.37 (0 3.35 (0 3.19 (0 2.60 (0 1.55 (0 
*3.36 (0 *3.34 (0 *3.19 (0) *2 5's (0 *1.55 (O "5.79 ( l *1.37 1 *2.19 2 *1.62 rt *4.90 (—6 
1.0 *1.37 (0 *1.37 (0 *1.34 (0 *1.14 (0 "6. 35 l *1.67 1 *1.73 (-—2 "1.02 3 
2 0 *3. OS 1 *3. ON 1 "3.06 ( l ww ] *1.42 1 *2.06 2 *5.93 ( 4 
( "7.20 3 “7.20 3 “7.21 (—3 *7. O08 $ *3. 48 3 "1.04 (—4 *7. 67 7 
10.0 *2. 61 *2. 61 5 *2.61 (—5 *2.61 ( "1.27 "2.71 (-—8 
p 
h=0 0.2 0.5 OS 1.0 1.2 15 2.0 0 10.0 
( 06 (O 4.98 (0 4.53 (0 3.69 (0 2 6 (0 2.26 (0 1.49 (0 &. 26 1 1.27 l 3.15 (—2 
0.01 1408 (0 1.90 (0 4.46 (0 3.62 (0 2.90 (0 2.21 (0 1.45 (0 7.95 1 1.15 l 2.58 (—2 
0. 02 4.00 (0 4182 (0 1.349 (0 $55 (0 2 84 (0 2.16 (0 1.41 (0 7.66 1 1. 04 1 aan @ 2 
0.05 1.68 (0 1.6 0 4.19 (0 3.38 (0 2 69 (0 2.02 (0 1.30 (0 6.84 l 7.74 2 1.16 (—2 
0.1 4.34 (0 4.27 (0 3.87 (0 3.11 (0 2.45 (0 1.81 (0 1.13 (0 5. 68 1 4.74 (-—2 4. 30 3 
2 3.73 (0 4, Ot 0 3.32 (0 2.63 (0 2.03 (0 1.47 (0 8.69 ( 1) 3.94 1 1.79 2 
( 2 38 (0 2.34 (0 2.11 (0 1.63 (0 1. 21 0 &. 07 1 1.09 1 1. 38 1 
"2.38 (0 *2. 34 (0 2 11 (© *1.63 (0 *1.21 (0 *8.07 ( 1 *4.08 ( 1 *1. 36 1 of $ 6. 69 5 
1.{ 1.16 ( 1.14 (0 1.03 (0 7.78 ] 5.43 l 3.26 (—1 
1.16 (0 “ae oo *1.03 (0 l "5.4 ] *3. 26 1 *1. 30 1 *2.73 2 1.8 
0 *2 93 2 OO | *2 69 l *2 03 1 0 i *6). 64 ? "1.73 (—2 *1. 62 3 
( ° 20 ; *7.04 $ *5. 76 $ 29 ; *1.09 $ *9.00 (—! *6. 03 7 
10.0 *2. 61 *2. 61 *2.61 ( *2. 31 22 5 *2.29 (—6 *2. 52 8 
} { F { o.2 0 os if 2 l 2 ) { oo 
2.18 (0 » 1 0 1.99 (0 1.72 (0 1.51 (0 1.30 (0 1.01 (0 6.66 (—1 1 3.13 (—2 
0.0 »15 (0 2.12 ( 1.96 (0 1.70 (0 1.49 (0 1.26 (0) 9.95 (—1 6. 52 1 l 2.83 (-—2 
(2 2.13 (0 2.10 (0 1.94 (0 1. & 0 1.47 (0 1.20 (0 9.78 (—1 6 1 1 2.56 (—2 
0.05 2 05 (0 202 L887 (0 1.61 (0 1.41 (O 1.20 (0 9.29 l 5 l 2 1.90 2 
( 1.93 (0 1.90 (0 1.76 (0 1.51 (0 1.31 (0 1.12 (0 8.53 (—1 5 1 2 1.15 (—2 
0.2 1.71 (0 1.69 (0 1.55 (0 1.32 (0 1.14 (0 9.62 (—1 7. 20 l 4.: | 2 4.24 (-—3 
( 1.20 « 1.18 (0 1.07. (0 8.96 (—1 7.56 l 6.18 (—1 4.37 1 2.35 (—1 1.02 2 
*1.21 (0 “1.18 (0 *1.07 (0 *8.904 (—] *7. 55 l 6.17 (—1 *4.37 (-—1 *2. 42 l *9.31 (—3 "3.96 (—4 
1.0 6.04 I 6.51 (—1 5.85 l 4.73 1 3. 87 ] 3.02 (—1 1.95 (—1 8.65 2 
"i. 04 l 6.51 1 "5.84 (—1 *4.73 l *3.87 (—1 *3.04 (-—1 "1.95 (—1 *8.65 (—2 *1.06 3 *2.27 (—6) 
2.0 2.0, l Ol | 1.77 (—1 1. 37 1 
2. OF l 2.01 (—1 "1.77 (-1 *1. 38 l *1.06 1 *7.67 (—2 "4.22 2 "1.29 2 *6. 65 ( 
0 6. 54 3 7 (-—3 *5.70 (—3 *4.21 (-—3 *2.95 (—3 "1.79 (-—3 "6. 51 ( 4) *7.40 (—5 
10.0 *2. 58 5 *2. 56 (—5 *2.41 (—5 *1.78 5 *1.14 5 *5.57 (—6 *1.13 (—6 *2.66 (—8 
f 
h=2.0 0 0.2 0.5 O.8 1.0 1.2 1.5 2.0 5.0 10.0 
) 7.01 l 6.96 (—1 6.71 (—1 6 (—} . 91 (—1 5 | 4.90 (—1 3. 92 l 1.10 (—1 3.03 (—2 
0.01 6. 94 l 6.89 | 6.63 (—1 6.20 (—1 ) l 5 l 4.84 (—1) 3. 86 (—1 1. 07 l 2.88 (—2 
0.02 6.86 l 6.81 1 6.56 (—1 6.14 (—1 l l 4.78 (--1) 3.81 (—1 1. 04 I 2.74 (—2 
0. 05 6. 65 l 6. 60 1 6.35 ( 1) 5.93 ( l ( l 5.2 ] $60 ( l 3. 65 1 9. 60 2 2.35 ( 4 
0.1 6. 31 l 6. 26 l 6.02 (—1) 5.61 (—1 (—1 i 1 4.32 (—1 3.40 (—1 8.40 (—2 1.83 (—2 
0.2 9. 67 l ». 63 l 5.40 (—1 5.02 (—1 (—]} 4.¢ ] 3. 81 (—1 2.96 (—1 6.44 (-—2 1.10 (—2 
0 4.13 l 1.09 1) 3.91 (—1 3.60 (—1 3.34 (—1 3.05 (—1) 2.62 (—1 1.95 (—1 2.91 (—2 2.41 (—3) 
*4.14 | *4.10 (—1 *3.90 (—1) *3.59 (—1 *3.34 (—1 *3.07 (—1 *2.62 (—1 *1.95 (—1 *3.00 (—2 *2.18 (—3 
1.0 2. 44 l 2.41 (—] 2.28 (—1) 2.06 (—1 1.89 (—1 1.70 (—1 1.41 (—1 9.74 (—2 7.83 (—3 
*2. 43 l *2.41 (—1 *2.28 (—1 *2.07 (—1 *1.89 (—1 "1.70 (—1 "1.40 (—1 "9. 74 2 *7. 69 3 *2. 35 4 
2.0 8. 47 2 8.36 (—2 7.79 (—2) 6.85 (—2 6.09 (—2 5.29 (—2 4.12 (—2 
*8. 44 2 *8.33 (—2 *7.79 (-—2 *6.85 (—2 *6.09 (—2 *5 2 *4.12 (—2 *) 50 2 *6. 16 4 
0 *3. 59 3 *3. 52 3 *3.17 (-—3 *2.60 (—3 *2.16 (—3 *1.72 (—3 "1.14 (-—3 *4.77 4 *2. 28 7 
10.0 *1. 89 F *1.84 (—5 *1.61 (—8 "hae 5 *9. 55 (—6 *6.86 (—6 *3.45 6 *9.02 (—7 





TABLE 3. Unscattered component (4x/o)I°(r,p,h) of the radiation flux calculated (a) using table 1 ef this work and table 2 of 
reference [3] (numbers preceded by asterisks) and (b) using table 2 of this work (no asterisks) —Continued 
tra) T° (p, h 
h=5.0 0 0.2 5 OS 1 l 20 50 10.0 
0 1. 23 l 1. 23 ] 1. 22 l 1. 20 l 1.19 7 1 ] 1.0 6.28 (—2 2. 52 2 
0.01 1.22 (—1 1. 22 l 1.21 1 1.19 Is ( 1. OF l 6.2 2 2. 51 2 
0. 02 1. 21 l 1. 21 l 1. 20 ! 118 ( 1. lt 1 1.0 G.11 (—2 2.41 4 
0. 05 1.17 l 1.17 1 1.16 l 1.14 l 3 ()S Lo l 5. 85 2 2. 25 4 
0.1 1.11 ( 1 1.11 l 1.10 l 0O9 l O7 0 7d 9.60 2 44 (—2? 2. 02 2 
0.2 1.01 (—1 1.01 (—1 9. 96 2 9.8] 2 9. 66 2 9, 50 2 9, 20 2 8.62 (—2 4.74 (—2 1.61 2 
0 7. 44 2 7 2 7.3 2 4.00 2 7.10 , 6. Yt 2 ( l 2 6, 23 2 3.10 (—2 &. 26 ] 
7.24 2 7.2 2 7. 24 2 7. 20 2 7.14 2 *7. 03 2 6. 80 2 6. 24 2 *3.14 (-—2 8.00 ( 
( 4.49 2 4.48 2 4.43 2 4.33 (—2 $. 20 2 4.14 2 3. 97 2 3. 62 2 1.53 (—2 2.71 3 
*4.50 (-—2 "4.49 2 *4.44 2 *4.33 2 *4. 25 2 *4.14 2 *3. 96 2 *3. 62 2 *1.53 (—2 *2. 72 3 
2. ( 1. 64 2 1.63 2 1. 60 2 1.56 (—2 1. 42 2 1. 47 2 39 2 23 2 3.76 (—3 
*1.66 (-—2 "1.66 (-—2 *1. 62 2 1. 57 2 *1. 52 Z *1.47 2 7 ta} 2 "1.22 2 *3.77 ( $ *2. 94 4 
0 7.91 { 7. 87 4 fit { 9 { 6.97 { 6. 60 { 7 { 1.83 } 71 *3.77 7 
10.0 *5 09 6 *5.05 (—6 *4,.83 6 *4.45 (—6 { } i *3.78 ( *3. 20 6 *) 27 t ay S 
i 1>4 0 0.2 ( O& ) ] ) 
0 3.13 2 3. 12 2 3.12 2 3.11 310 2 3. 08 2 rt 2 3.01 2 2.51 2 1. 57 2 
| 3.09 2 3.09 2 3.09 3. OS 2 3. OF 2 4. 08 2 $03 2 2. OS 2 2. 48 2 1.55 ” 
0. 02 3. 06 2 3. 06 2 3. 06 2 3. 04 2 3. 03 2 3. 2 2 aL 2 2. 95 2 2.4 1. 53 
0.0 2 97 2 2. 97 2 2. 97 2 2. 95 2 2.94 2 2. 93 } 2 2. sf 2 >. 37 2 1. 46 2 
0 2.83 2 2. 83 2 » &2 4 2.81 2 2.80 2 » 79 y. 2 2.72 2 2.24 2 1.36 2 
0.2 56 (—2 2. 5 2 2.55 (—2 2. 54 2 2.53 (—2 2.52 (—2 2 2 2. 45 2 2.00 (—2 1.18 (—2 
0 gg 2 1,89 2 1.89 2 1.88 2 87 2 St 2 & 2 x0 2 1.43 2 792 (9 
7¢ 2 7¢ 4 1.77 2 1.79 2 be 2 1.8 2 82 sd 2 1. 40 2 7. Be 5 
1.0 2 1.15 2 1.14 (—2 1.14 2 1.13 2 1.12 1.11 2 1.08 (—2 8. 20 3 3. 82 3 
2 "1.15 2 1.15 2 1.14 (—2 at BE 2 *1.13 2 a 2 "1. 08 2 *8.20 (-—3 *3. 80 
2.0 4.21 3 4.21 (-—3 4.19 3 4.16 (—3 4.13 3 4,09 4.03 3 3. 90 3 2 68 3 9.31 4 
4.33 3 *4. 32 3 4. 30 $ 4.25 $ 4.2 3 *4.16 *4.08 3 *3. 91 3 *2. 67 3 *9. 33 ( 4 
5.0 *2.08 (—4 *2.08 (—4 *2.06 (—4 *2.04 (—4 2. 04 t 2.01 4 13 i 1. 82 4 *9.37 (—5 *1.36 ( 
10.0 39 6 “1.38 ft 1. 37 ( *1.34 ( 1. 31 t *1.28 f »”) 6 10 t 3.53 7 *1.26 (—8 
TABLE 4. Comparison cf converaence in calculating the un- A physical interpretation of the rather wild 


scattered component (4x/a)1°(r.p,h) in the case r=ul=1, for | oscillations in the first few terms of the 7” series is 
the three situations (.A4), (B), and (C) indicated in figure 


For each situation the first column contains partial sums that the zeroth term would be the result for a bare 
disk source (r=0), from which the radiation ab- 


Le sorbed in the barrier of thickness 7+ must be sub- 
> (l+4)gi( 4) pil p,h tracted. This subtraction is accomplished by a 
. succession of undershoots and overshoots which are 
and the second colun contains sums quickly damped oul if rT 1S small enough. 
- 
>> an(p,h)r*. 3.4. h=O, p>1: Power Series Solution In or Neor 


0 Plane of Disk 


(A): p=0.5, h=0.5 R): p=1, h=1 C): p=2, h=1 If both source and detector are embedded in a 
x t- | homogeneous attenuating medium, the series (34) is 
“| Legendre | r=series | Legendre | rm series | Legendre | r=series | valid except for h=0, where the concept of “barrier 
sum sul sum sum sum sur 


thickness” + becomes meaningless. As h—>0, g,—> 
as 1/h”", but at the same time 7”"->0 as h”, while 








0 0. 34031 4. 5353 0. 12316 1. 5118 ) 66582 ’ a - ‘ . ; . 

$4660 2.6711 31658 0. 6170 0.75149 | qo(p,0) 1s a finite but indeterminate quantity of the 

2 QO567 3 38273 0. 9538 0. 81932 » } . 

3 99211 -() 37916 1441 0. 38701 form oO O;r, 

i. aoe tC rock, aa pees The series (34) can still be used if the indeter- 
234 U. 4" ONG dOlz - W218, . A in ns 

f 1. 02365 1. 1264 38665 3932 0.12798 | minacy is removed from the product q,,7" by multiply - 

7 1. 02628 0. 9943 38793 3850 07242 . - ET. " ; ee rn 

8 1. 02590 1. 0367 38753 | 3869 09079 Ing and div iding WT» and , respectively, by h 5 

9 1. 02615 1. 0243 38718 | 3865 OR520 ones ine in a serie ; 7 

10 | Loaei9 ose sara | aor pss’ | resulting in a series of the form 

11 1. 02655 1. 0268 BR688 3865 OR6408  ORG41 

12 1. 02677 1. 0270 BRORS 3865 O86602 08650 

13 1. 02685 1. 0269 38672 3865 086496 O86 48 ' 7 ’ : i 

/ (r,p,h (c/4r) > q n(p sh) (u-1)”, (47) 


n=0 





in which d'n(p,h) h"qr(p,h) and pl t/h=upt/h 
uh/h==the radius of the disk source measured in 
mean-free-paths of the attenuating medium, and 
q'o p,h )= ol p,h ). 


For h=0, p>1 the expression in (37) reduces to 


9 


qo(p, 9) In ( =. ): (48) 


p 
4. Buildup Factor 


When both source and detector ure embedded in 
a homogeneous medium, the secondary, or scattered 
radiation contribution J* is often taken into account 
by means of a buildup factor [6] 


B(ur)=(U°+-T*)/ 1°. 
Such a buildup factor may then be used to evaluate 
the total flux J from a point isotropic source ac- 
cording to 
T=1°+T°=I°B(ur)=(0/4x)B (ur) exp(—ur)/r?. (49) 


In references [S] and [9] B (ur) is expressed as a set 
| 
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of polynomial coefficients 8;, in the notation of 
reference [9], as defined in equation (2), with N=3. 

Now, since the integration of (33) was performed 
term by term over the expansion 


exp (—7/cos 6)=exp (—ur) >5(—1)"(ur)"/n!, (50) 


no new integrands appear if these series terms are 
multiplied by terms 8;(ur)' of the polynomial repre- 
sentation of Bur). if this multiplication is per- 
formed, the total radiation flux as defined in (49) 
may be written as the series 


, 


(0/4) >5B' n Unt (51) 


in which r=yph. The modified buildup coefficients 
B’, are generated from 8,’s, such as are given in 
[S] or [9], as 


n\/(n a)!} B, (52) 


in which n!/(n—7)!=0 for 7>n. For N=3 the 
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Fraure 5. Some circular aperture or disk source situations to which results are applicable 
(a): [1°+]* calculable using pi(p,h), qi(r) coefficients. (b): 1°+-1* calculable using 8; ,q;(p,h 
coefficients. (a-f): I° calculable using either pi(p,h) and qi’, or qi(p,h) coefficients. 
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first five modified buildup coefficients are 


By Bo= l 7 


, 





B; Bo— Bi 
8; =By—28, +2 - 18, , a 
(9) 
Bz = Bp — 361 +3 -28.—3-2- 1B, 
and 
B,=Bo—46: +4 - 3Bo—4-3-265. J 
Such a set of 8; coefficients, once generated, 


may be applied as in (51) to either the disk-source 
dn(p,h)’s in table 3 of this work or the rectangular 
source qg,(a,6)’s in table 3 of reference [3]. 


5. Discussion 


Situations to which results in this paper may be 
applied are schematized in figure 5. 

The Legendre method of section 2 can be used 
to calculate J=J°+J* in such cases as figure 5a 
in which the source is distributed in a uniform layer 
(or lavers) parallel to the exit surface of the medium. 
If g; angular harmonic coefficients are used which 
have been calculated assuming an infinite medium, 
the result will be an overestimate because of the 
depression, in the realistic case, of flux near a 
boundary. This flux depression results from the 
one-way nature of the radiation traversal across 
the boundary. For gamma rays the infinite medium 
gis may be a rather good approximation, whereas 
for neutrons more realistic g,’s including the boundary 
effect may have to be used. 

The unscattered component J’ may in principle 
be calculated for all situations in figure 5a-f using 
either the Legendre p,(p,h)’s of section 2 or the 
dn(p,h)’s of section 3, but in practice the Legendre 
series (6) converges better for barriers of the order 
of one mean free path or thicker while the q,(p,h) 
series (34) converges better for barriers of less than 
a mean free path. 

Adaptation of polynomial representations of point 
isotropic source buildup factors, as suggested in 
section 4 for use with the q¢,(p,h)’s of section 3 for 
computation of J=J°+J*, is strictly valid only 
for the situation in figure 5b, although some qualita- 
tive results may be obtained in the other situations 
depending on how little the physical situation 
departs from figure 5b. 
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The Bauschinger effect in alpha brass, specifically the lowering of the yield strength in 
the direction opposite to the preceding plastic deformation, was studied by an examination 
of both tensile and compressive elastic limits measured at the quarter-cycle stages throughout 


a complete cycle of uniaxial plastic strain of one percent amplitude. 


X-ray measurements of 


axial residual directed microstresses indicated that the latter could be correlated with the 
decrease of the elastic limit measured in the direction opposite to the preceding deformation 


only at the first quarter-cycle of deformation. 


After 3 quarter-cycles the axial residual 


stress was always tensile regardless of the direction of previous deformation; whereas the 
elastic limits continue to show strong directionality. 


1. Introduction 


Although it is nearly S80 vears since Bausch- 
inger [1]' first reported the effect which may be 
described very broadly as the anisotropic modifica- 
tion of the plastic properties of a metal specimen 
that has undergone a preceding plastic deformation, 
a simple quantitative explanation for the phenome- 
non has not yet been established. Perhaps the 
most common aspect of this effect is the raising of the 
yield strength of a plastically deformed metallic 
specimen in the direction of previous deformation 
and the lowering of the vield strength in the Oppo- 
site or reverse direction. The basis of an explanation 
often advanced for this phenomenon is the dif- 
ference in degree of susceptibility to slip in differ- 
ently oriented crystal grains or regions of differing 
hardness within the specimen [2]. This inhomog- 
eniety of slip is thought to give rise to a system of 
residual internal upon release of the de- 
forming load; and, upon reapplication of load, these 
residual stresses presumably modify the external 
stress required for vielding. 

Residual stresses in metals can be calculated from 
measurements obtained by means of X-ray diffrac- 
tion [3]. If the stresses are random in sign and 
magnitude in the diffracting material, the X-ray 
diffraction lines are observed to be broadened; if 
the stresses are nonrandom, the peak of the line 
will be shifted. This latter effect, the change in 
Bragg angle of diffraction after plastic deformation 
of a specimen, implies a nonzero net stress in the 
diffracting material, which is usually observed to 
be of such a sign as to oppose the prior deformation. 
S. L. Smith and W. A. Wood [4] observed such 
residual stresses by means of X-rays and believed 
them to be related to the Bauschinger effect. 

Recently a detailed quantitative study was made 
of this type of residual ‘‘oriented microstress”’ and 
its relation to the Bauschinger effect by B. M. 
Rovinskiy and V. M. Sinayskiy [5]. These experi- 
menters worked with three groups of specimens of 


stresses 


' Figures in brackets indicate the literature references at the end of this paper. 


40X steel (Russian designation), each group having 
had a different thermal or strain-rate history. The 
specimens, after having been annealed, were com- 
pressed uniaxially by different amounts, ranging 
from 0.30 percent to 1.46 percent. Although the 
authors do not give details of their X-ray procedures, 
they state that back-reflection X-ray patterns of the 
(310) spacing (using cobalt radiation) after these 
prestrains enabled them to calculate the residual 
microstress, gg, In each case. The direction of op 
is not specified, but it is presumed to be axial. The 
conventional 0.2 percent yield strengths in com- 
pression and in tension were then determined, and 
the difference, Ac, between them was obtained for 
each degree of prestrain. Both the yield strength 
difference and the residual stress were found to go 
through a maximum at a deformetion near | percent. 
The value of prestrain that resulted in these simul- 
taneous maximums near 1 percent depended on the 
thermal and strain-rate history of the specimen. 
The ratio of cp to Ao was found to differ from one 
group of specimens to another, but to vary by only 
a few percent within a group. These authors 
observed that this ratio, in most cases somewhat 
greater than unity, was comparable to that observed 
for the ratio of the X-ray stress, measured on a 
specimen under load, to the applied mechanical 
They concluded that, if one keeps in mind 
the coefficient of correspondence between stresses 
determined by X-ray and mechanical methods, one 
can predict the magnitude of the reduction of the 
vield strength during reversed loading of a specimen 
from the size of the residual oriented microstress 
observed with X-rays. Thus the Bauschinger effect 
after uniaxial plastic deformation is, according to 
this picture, the result of the establishment of a 
residual oriented microstress that can be summed 
with the macrostress created by the external load. 
Or, in equation form, 


stress. 
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where o, is the yield strength of the annealed mate- 


265 








rial, oo is the yield strength observed in the reverse 
direction after a plastic deformation of the speci 
men, @ is constant nearly equal to unity, and gp is 
the oriented axial residual microstress messured by 
X-ray diffraction. Or, from another viewpoint, the 
reduction of the yield strength is 


Ag=64—6)—acp. 


Since the residual stresses in a specimen are balanced 
internally by stresses of the opposite sign, there is a 
further corollary implicit in this model, that is, that 
the region or crystals that support the additive stress 
observed by means of X-rays is that same region in 
which plastic flow first begins, thus determining the 
yield strength. 

This article on the Bauschinger effect by Ro- 
vinskiy and Sinayskiy was of sufficient interest that 
it was decided to investigate the Bauschinger effect 
and the state of residual stress at appropriate stages 
in specimens that had undergone a complete cycle of 
plastic deformation, unlike those of Rovinskiy and 
Sinayskiy, which had been plastically deformed in 
one direction only. 


2. Experimental Material, Procedure, and 


Results 

When a metal specimen containing a_ residual 
internal stress system is sectioned, it is assumed that 
at the exposed surface the stress component norma! 
to the surface vanishes and the components of stress 
paratle! to the surface are not affected. A few 
vears ago an X-ray study of residual stresses on 
various sections of plastically deformed iron by 
C. J. Newton and H. C. Vacher [7] vielded results 
consistent with this long-held hypothesis. The 
procedure in that study emploved an X-ray beam 
normally incident on the sections of the specimen 
with the diffraction pattern recorded on photographic 
film. The observed lattice spacing was compared 
with that of the annealed material. A similar pro- 
cedure was used in this study of sections of alpha 
brass specimens, except that the diffracted X-rays 
were detected by means of a diffractometer (counter) 
technique instead of by film. 

The alpha (Cu=70, Zn=30, Pb<0.07 percent) 
brass specimens which served both for elastic limit 
measurement and for X-ray residual stress measure- 
ments, were machined from 1 by *% in. bars, to be 714 
in. long with a gage length of 1% in. and a reduced 
section measuring 0.625 by 0.373 in. X-ray examina- 
tion indicated no appreciable preferred orientation 
in the material. The plastic extension was _per- 
formed on a hydraulic testing machine using Templin 
grips, the plastic compression was made possible by 
careful machining of the ends of the specimens and 
by the use of massive steel guides around the speci- 
men undergoing compression. During the course of 
interrupted straining, resistance wire strain gages 
were used to follow the amount of strain. The 
routine of the mechanical testing and strain measure- 
ments, along with a discussion of the difficulties that 


sometimes arise when one attempts to measure very 
small permanent sets, and hence elastic limits, with 
resistance gages, has been previously reported [8]. 
The elastic limit was taken to be the stress required 
to give the specimen a permanent set of 21075, us 
indicated by gages not previously subjected to large 
plastic strains. 

The values of the limit for various conditions of the 
brass specimens are listed in table 1. It is difficult 
to specify precisely the reliability of these values, 
most of which are the average of several observations. 
The precision of the testing machine is high, but 
variations from one specimen to another may be 
large. A reasonable estimate of average uncertainty 
might be about + 1 10° psi. 


TABLE 1. Elastic limits and oriented residual stresses 


Plastic deformation state of specimen Elastic limits in 10° psi 
Axial 
residual 
Net } Cum. | In com stress in 
Description prestrain tic pre In tension pression 108 psi 
str I 
Annealed { ( 14 13 Reference 
zero). 
I cyel ia] deformation extension 
I 1 15 3 ‘ 
E-( ( 2 4 lf 2 
E-C-C l 5 yi) 5 
E-C-C-E ( i | 3 7 
re 
Negative evel nitial deformation compression 
( 1 lt rs 
C-] ( 17 3 3 
C-E-E +4 
C-E-E-( ( 1 +S 


E= Extended plastica 

C=Compressed plast ! 

Average uncertainty approximately +1X10° p 

Different specimens were used for the X-ray residual stress 


5 measurements 
ist column) than were used for the elastic limit determinations 


One can see in table 1 that essentially symmetrical 
results of elastic limits with respect to plastic exten- 
sion or compression were obtained when the direc- 
tions of limit determination and prestrain are taken 
into account. The effect of plastic deformation was 
to raise slightly the elastic limit measured in the same 
direction as the deformation and to lower by a large 
amount the limit measured in the opposite direction; 
this is in accord with the Bauschinger effect. When 
a second plastic deformation, reversed in direction, 
was made, the elastic limits were affected by roughly 
the same amounts relative to the values pertaining 
to the original material in the annealed state, but 
reversed in direction. This demonstrated that it is 
the direction of the immediately preceding deforma- 
tion that is of major importance with regard to the 
elastic limits The third step in the positive cycle 
of plastic deformation, at a minus 1 percent net 
strain, showed a small increase in both tensile and 
compressive elastic limits. And finally, the last 
step of the cycle, an extension of 1 percent to reverse 
the net strain to zero, resulted in a strong reversal 


| of the elastic limits, as expected. 
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Elastic limit determinations were made during 
only the first two steps in the negative cycle of de- 
formation, since the results were obviously svmmetri- 
cal with those of the positive cycle. 

Longitudinal and crosssections were cut from the 
reduced sections of new test specimens at quarter- 
cycle stages of one complete cycle each of positive 
and negative 1 percent plastic deformation. The 
sections were mounted without heat in acrylic cement 
and carefully mechanically polished, electro-polished, 
and etched until all surface disturbance was removed. 
The X-ray diffraction measurements were made on a 
commercial horizontal diffractometer with cobalt 
Ka, characteristic radiation, using a proportional 
counter as a detector. Each specimen was con- 
tinuously rotated about an axis normal to its surface 
in order to bring more grains into diffracting position. 
The position of the peak of the (400) line, near the 
20 diffraction angle of 153°, was determined by the 
analytical method of three-point parabola fitting of 
D. P. Koistinen and R. E. Marburger [9], as described 
in the SAE Information Report TR-—182 [10]. Cor- 
rections were made in the intensity at each point for 
the Lorentz and polarization factors, and all values 
of the calculated Bragg angle were corrected for 
specimen temperature variation. The precision of 
peak determination was approximatly + 0.02 
in 26. 

The value of the crystal lattice spacing observed 
normal to the surface of a given section compared to 
that of the annealed material enabled one to caleu- 
late the sum of the principal stresses parallel to the 
surface of the Because of the uniaxial 
manner of deformation, it was assumed that one 
principal stress o,, Was axial in the specimen and 





section. 


two equal principal stresses, each or, were at right 
angles to the axis. Hence, measurements on the 
cross section vielded a Value, 2o7, twice the trans- 
verse stress, and those on the longitudinal section 
vielded (o,+ or), the sum of the axial stress plus the 
transverse Therefore, the axial residual 
stress was calculated for each pre-strain condition 
by subtracting one-half the sum of the stresses found 
in the cross section from the sum of the stresses found 
on the longitudinal section. 


stress. 


That is, Ox—(Onq+ O07) 
The resulting values of axial residual 
stress are listed in table 1. Because of the involved 
nature of the calculation and the many possible 
sources of error, only an estimate of average un- 
certainty, +1 10° psi, is expressed as a measure of 
precision of the results. 

The behavior here again shows a kind of symmetry 
between the cvele of plastic deformation begun with 
an axial extension and the cycle begun with axial 
compression. After the first 1 percent plastic de- 
formation, a residual axial stress is observed that is 
opposite in sign from the deformation; that is to say, 
after extension, a residual compressive stress is 
observed, and, after compression, a residual tensile 
That is the common behavior 
that has often been reported. After the second 
quarter step in the deformation cycles, the observed 
stress, though smaller in magnitude, has reversed 


l ) 
, (Zoey) 


stress is observed. 


sign so that it continues to be directed opposite to the 
deformation. After the third and fourth quarter 
steps, however, the residual stress is tensile in all 
cases, regardless of direction of deformation. This 
behavior is shown graphically in figure 1. 

These residual stresses, resulting in a change in the 
X-ray diffraction peak position, are oriented micro- 
stresses, indicating a nonzero average component of 
stress present in the coherently diffracting material 
of many grains of the specimen in the X-ray beam. 
Another type of microstress, in which the residual 
stresses may be random in magnitude and sign as 
one goes from grain to grain, results in a broadening 
of the diffraction line. Two other factors in the 
specimen may increase line breadth: random micro- 
stress that varies within each grain (microstress of 
the ‘third kind,” as it is sometimes called) and very 
small particle size. Previous study of individual 
diffraction spots on stationary film patterns made 
with a monochromatic misorientation goniometer 
[6], using some of the same specimens, however, 
revealed that neither of these factors were significant 
under the conditions of this investigation. 

Line breadth measurements were made on the 
(400) diffraction line of the brass using cobalt Ka, 
radiation, employing the diffractometer with a 3° 
beam slit, a 0.1° receiving slit, Soller slits, scanning 
rate of 0.2° per min, rate meter time constant of 
8 sec, recording chart speed of 0.4 in. per min, and 
specimen spinner in use. Since the a,—a, doublet 
was not fully resolved in these patterns, one-half 
the breadth of the low-angle side of the a, line was 
measured at half-peak intensity and multiplied by 
two. Measurements on the cross and longitudinal 
sections of a specimen were very nearly equal and 
were averaged together. The observed breadth was 
corrected for instrumental broadening by the method 
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described by H. P. Klug and L. E. Alexander [11] 
for high-angle diffractometer lines. The maximum 
random microstress, without regard to sign, is given, 
according to B. D. Cullity [12], by 


KE 


) 
4tané (6 


omz— 
where %m: is the maximum random tensile or com- 
pressive microstress, / is Young’s modulus of the 
material, @ is the Bragg angle ‘of diffraction, and 
8 is the corrected line breadth. The values of 8 and 
ou: are to be found in table 2. The precision in 8 is 
estimated to be about +0.05°, which is reflected as 
an uncertainty of +0.7> 10° psi in om. 

It is apparent that the maximum random micro- 
stress is of the same order of magnitude as the 
residual directed microstresses. One may also see 
that there is no significant difference in the random 
stress level in the annealed specimen and those of 
the first two quarter steps of both the deformation 
cycles, after which the stresses rise markedly. 


3. Discussion 


The essential conclusion of Rovinskiy and Sinay 
skiy was that the Bauschinger effect can be explained 
by:summing the residual directed microstress, de- 
termined by X-ray diffraction, with the externally 
applied macrostress to account for the lowered elastic 
limit measured in the direction opposite to a preced- 
ing plastic deformation (which we shall call the “re- 
verse elastic limit”). The results of this study 
clearly demonstrate that this conclusion cannot with 
validity be extended to those cases where the im- 
mediately preceding deformation has itself been pre- 
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FicuRE 2. Decrease in elastic limit and arial oriented residual 
stress versus cumulative plastic strain through one complete 
cycle of defor mation. 
oa= Elastic limit of annea 
o.= Elastic limit observed 
or= Axial oriented residual 


led material 
in reverse direction after deformation. 
microstress measured by X-ray diffraction 


ceded by a plastic deformation in the opposite direc- 
tion. This inequality of the residual stress and the 
decrease of the reverse elastic limit, Ao, beyond the 
first quarter-cycle is presented graphically in figure 2, 
where both quantities are plotted against the cumula- 
tive strain. Throughout at least the first four 
quarter-cvcle steps of a cycle of deformation, the 
Bauschinger effect, as can be seen in table 1, depends 
to a very Iligh degree upon only the immediately 
preceding deformation. The effect of earlier de- 
formations seems be merely to increase slightly 
the difference between the forward and reverse elastic 
limits. An entirely different situation exists with re- 
gard to the X-ray determined residual axial stresses. 
As can be seen in table 1 and figure 2, only the first 
quarter-cycle shows an axial residual stress that ap- 
proximates the Bauschinger decrease in reverse elas 
tic limit. The residual stress at the second quarter- 
evele, +2> 10° psi, agrees in sign but not in magni- 
tude with the decrease, +10 10° psi, in the limit. 
Clearly something other than the immediately pre- 
ceding deformation is beginning to have an important 
effect. After the third and fourth quarter-cycles of 
deformation, the direction effect of the immediately 
preceding deformation is increasingly masked by 
this cumulative factor. At this point, specimens in 
the positive cvcle and the negative cycle are indis- 
tinguishable on the basis of the axial residual stresses 
observed, whereas the Bauschinger effect, that is, 
the directionality in elastic limit, is undiminished. 

The validity of the various methods of calculating 
residual stresses from X-ray measurements continues 
to be open to some question. A new approach has 
recently been used by D. M. Vasil’ev [13,14] to com- 
pute the internal stress system from measurements 
on sections. He assumes as usual that the stress 
components lying in the plane of the cut are not af- 
fected by the cut. He does not make the usual as- 
sumption, however, that the component normal 
the surface of the cut is relaxed to zero; but rather 
that, in the equations relating strain to stresses in 
this region reached by X-rays, it appears multiplied 
py a factor k, to be evaluated, which lies between 
zero and one. He computes k from data obtained 
from additional X-ray measurements. His resulting 
values of k vary typically from 0.2 in the 
aluminum (99.99 percent pure) examined with 
Ka radiation to 0.7 in the ease of Fe,C examined 
with Cr Ka radiation. <A surprising result of his 
measurements and calculations is that all principal 
both axial and transverse, after either uni- 
axial extension or compression, are always compres- 
sive. This is obviously in conflict with the usual re- 
sult that the residual axial stress is opposite to the 
direction of deformation, as was used by Rovinskiy 
and Sinayskiy for the basis of their explanation of 
the Bauschinger effect. It is also in conflict with 
our result that, after cyclic axial deformation, the 
residual axial stress becomes tensile regardless of the 
direction of the last axial deformation. 

No attempt was made in this study to evaluate k 
by the somewhat complex procedure of multiple in- 
clined incidence X-ray exposures of Vasil’ev. How- 
ever, calculations were made of the axial stresses 


case of 


stresses, 
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using a k of 0.4, Vasil’ev’s value for copper with 
cobalt radiation. These calculations yielded results 
containing both tensile and compressive stresses, at 
least qualitatively similar to those obtained by the 
accepted procedure and listed in table 1. These 
newer values were not added to the table, however, 
because they did not seem to add appreciably to the 
results already presented, and, moreover, there does 
not appear to be adequate evidence that the value of 
0.4 used for k was more nearly correct than the usual 
value zero. 

While the directed residual stresses seem gradually 
to lose any relationship with the elastic hmits as the 
plastic strain cycle passes the three-quarter mark, 
the values of the maximum random microstress cal- 
culated from diffraction line breadth measurements 
and listed in table 2 show an abrupt change after the 
third quarter step in the deformation cycle. Since 
no such marked change in the Bauschinger effect is 
observed at this point, it seems clear that the latter, 
at least in the cyclic case, is the result of some factor 
or factors other than either of these types of micro- 
stresses revealed by X-ray diffraction. 


TABLE 2. Results of line breadth 
measurements 
Plastic Corrected Max. random 
condition line breadth microstress 
specimen B(26 a Mr (108 psi 
rN 
. J 
Annealed 0. 35 5.2 
Positive cycle of deformation 
} 0.30 15 
E-C 35 2 
E-C-( () 7.5 
E-C-C-E 8.2 
Negative cycle of deformation 
Cc 0.3 2 
G-j 35 5.2 
C-E-F 7-1) 7.5 
C-E-E-C ) 7.5 
E=extended 1 ( compressed 1°% 
Average uncertainty in B approximately 
+-1).6), 
Average uncertainty in oMz approximately 


It may be that the Rovinskiy and Sinayskiy 


model of the part plaved by the oriented microstress 


is correct for the first deformation of an annealed 
material; but that, as the material is subjected to 
further strain of a ecyelic nature, another type of 
residual stress arises that is always tensile in char- 
acter and = increasing, perhaps ‘asympotically, in 
magnitude, while that associated with the Bau- 
schinger effect is outweighed or masked, after the 
third quarter stage, altogether. The rise of this 


conjectural new type of residual stress is also indi- 
cated at this point in the cyele by the increase in 


line broadening. The Bauschinger effect, in the 
meantime, is not affected at all; the elastic limits 
are as strongly directional at the end of a evele of 
deformation as at the beginning. One must con- 


| 
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clude, therefore, that an explanation of the Bau- 
schinger effect that is valid in all must be 
based on some model other than the additive action 
of those types of residual microstress that are 
revealed by X-ray diffraction effects of peak shift 
or broadening. 

The plastic properties of metals in general are 
being explained with increasing success on the basis 
of the dislocation theory of behavior as exemplified 
in the work of A. H. Cottrell [15], N. F. Mott [16], 
and many others. Some application of the theory 
has been made to the problem of the mechanism 
of the Bauschinger effect, as for example in the work 
of R. L. Wooley [17] and that of S. N. Buckley and 
K. M. Entwistle [18]. <A critical review of these 
theories is beyond the scope of this paper; it appears, 
however, that there does not yet seem to be a theory 
sufficiently complete to enable one to predict quanti- 
tatively the behavior of the forward and reverse 
elastic limits of a cyclically deformed polycrystalline 
specimen such as used in this study. No doubt the 
complexity of the physical situation in experiments 
such as this with polycrystalline brass is such that 


cases 


it would be difficult, if not impossible, to apply 
dislocation theory quantitatively and with rigor. 


It has been suggested that the Bauschinger effect 
is a result of the nonsymmetry of the stress potential 
associated with a dislocation arrangement after 
flow, as, for example, arising from a pile-up against 
obstacles on the slip plane or from the interaction 
of other dislocations on parallel slip planes. If 
this is true, the elastic stresses in the crystal lattice 
would probably be quite localized, on the scale of 
perhaps 10? or 10° interatomic distances; whereas 
the stresses associated with X-ray measurements are 


averaged over much larger dimensions. Even if 
the tool of X-ray microstress measurement were 


sufficiently sensitive to reveal the nature of a localized 
dislocation configuration after an initial deformation 
of the order of 1 percent, it should probably not 
be surprising that this information would become 
masked by more gross structural effects and would 
cease to be a m: Lor contributor to the X- ray micro- 
stress after the cumulative effects of reversed plastic 
strain past the middle of a strain cycle became 
important in the specimen. 


4. Summary and Conclusions 


The principal experimental facts reported in this 
study may be summarized as follows: 
(a) The elastic limits indicated a strong Bau- 


schinger effect in alpha brass throughout a complete 
evele of axial plastic deformation of 1 percent 
amplitude. 


(b) There was no marked change in the effect as 
the cycle progresses; only a slight hardening was 


apparent, 

The axial oriented residual microstress was 
not additive to the applied stress to account for 
the Bauschinger effect in the latter part of the 
cycle; in fact, after the third-quarter stage, the 
microstress appears always to become positive in 
sign (tensile). 
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(d) The random maximum microstress (X-ray 
line broadening) remains near the value found in 
the annealed material until it exhibits a marked 
change at the three-quarter cycle condition of the 
specimen. 

These facts point to the following conclusions: 

(a) The effect of the cyclic deformation on the 
elastic limits seems to require two mechanisms (at 
least): one to account for the forward limit hardening; 
a second to account for the large reduction of the 
reverse limit. Neither mechanism would appear 
to be strongly affected by the cumulative strain of 
one complete cycle of 1 percent amplitude of plastic 
deformation. 

(b) Since both the oriented and the maximum ran- 
dom microstresses measured by X-ray techniques 
show a marked change in character after three- 
quarters of a cycle of deformation while the Bausch- 
inger effect does not, they are not directly related to 
the Bauschinger effect, at least under these conditions 
of reversed plastic deformation. 

(c) The eyelic uniaxial working of this materia! 
leads to a stress in the axial direction. 
Determined by X-ray diffraction, this is a positive 
(tensile) stress. 


residual 
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A Study by Polarization Techniques of the Corrosion 


Rates of Aluminum and Steel Underground for Six- 
teen Months 
W. J. Schwerdtfeger 
(July 6, 1961) 
One aluminum and one steel specimen were exposed underground for 16 months in 


order to make a running survey of their corrosion rates as influenced by weather and time. 
Cathodie and anodic data were obtained periodically as a basis for calculating rates of 


corrosion. Automatie polarizing and recording equipment, housed indoors, was used for 
making the measurements. The method is adaptable for field use in making corrosion 
rate or soil corrosivity studies where commercial power sources are not available. 

1. Introduction Before capping, the ends of the specimens were 
rounded with a file and smoothed with fine abrasive 
The present investigation was carried out to | cloth. The outer and inner surfaces of the steel 
evaluate polarization techniques for measuring | tube were cleaned with a motor-driven wire brush 
rates of corrosion of aluminum and_ steel under- | to remove any trace of seale or superficial rust. 
ground. An exposure period of at least 1 year was | The surfaces of the aluminum tube were polished 
decided upon in order to study the effects of time | with fine abrasive cloth, and scrubbed under hot 
and climatic changes on the corrosion rates of these | water. Each specimen was then weighed to the 
materials in one environment. For convenience, | nearest milligram. An insulated stranded-copper 


both materials were placed underground outside | wire was then soldered to an inside edge of each 
of the laboratory which housed the measuring | tube and the soldered joint coated with paint. 
apparatus, Next, the inside surfaces of the tubes were coated 

The theory upon which the measurements were | with a film of rust-preventive oil and the rubber 
based has been applied to the measurement of local caps pressed on and sealed. The lead wire came 
action corrosion currents on steel exposed to aqueous | through a hole bored in the center of one of the caps. 


environments in the laboratory [1].! Briefly stated, | The area where the wire protruded was sealed with 
the method consists of noting the value of cathodic | bitumen. 

and anodie current applied in increasing amounts The two specimens were buried underground 
to a corroding specimen which coincides with a | outside of the laboratory, Washington, D.C., in 
change-in-slope (break) in the respective polariza- | soil having a resistivity of about 7500 ohm-em. 
tion curve and which is indicative of the current | They were spaced 3 ft apart in a vertical position m 
necessary to stop local-action corrosion. holes dug 30 in. deep with a post-hole auger. An 


A logarithmic relationship between the value of | auxiliary electrode (steel channel 1.5 in. by 4 in. by 
current at the break in the controlling polarization | 19 in.) was similarly buried about 20 ft away from 
curve and the rate of polarization (based on the | the specimens. The rubber-covered wires from the 
curve) was observed to be reasonably applicable | specimens, auxiliary electrode and one for a _refer- 
to both the steel and the aluminum. A similar | ence electrode, led to the laboratory where polariza- 
relationship was described previously [2] together | tion measurements were made. 
with related work of other investigators [3]. } 

The technique described here could be put to 2.2. Electrical Measurements 
practical use in the field for measuring corrosion 
rates (based on weight loss), screening metals and 
alloys psa long-time exposure in aqueous environ- 
ments, or for measuring soil corrosivity. 


Potential and polarization measurements were 
made periodically (on 40 occasions) during the 16 
months that the specimens were exposed. The 
separate potentials were measured with a 200,000 
ohm/v voltmeter and the polarization data were 
2. Experimental Procedure obtained with a two-pen (current and poteatial) 
strip-chart recorder. The circuit used for balanc- 
ing out IR drop, automatically applying polarizing 
Two specimens were used, one of low carbon steel current, and recording current and potential, has 
and one of high purity aluminum. Each specimen been previously described : [4]. A eh cag ae 
an de. the: deees of & Walon. Gabent 2.40 ie. ok. be | sulfate half-cell was used for the reference electrode. 
12 in. long) with rubber caps sealed on the ends, When measurements were to be made, the reference 
leaving about 0.4 ft? of outer surface for exposure electrode was placed on the surface of the earth 

about 3 ft from the specimens and between them 
and the auxiliary electrode. The auxiliary electrode 
Figures in brackets indicate the léerature references at the eisd of this pape: served either as an anode or a cathode. 
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2.1. Preparation and Exposure of Specimens 


to the soil. 








Apparatus for obtaining polarization data has | —4 T ™ ae ee i i a ee ae ee Ser 





also been built for field use by using dry batteries | L ; 
and controls for manually varying the polarizing ~~ a 
) . . : A: ; : 4 e-e-— ° —e Os -O—__o——_o- 
current. Portable indicating-type instruments con- t - ae 4 
sisting of a potentiometer (galvanometer-type null - 2 J 
indicator) and a milliammeter are used. An earth- 
auger (1 in. diam) appropriately fitted with a cross- Ps 
bar for turning, is used as the auxiliary electrode. sa ow 
When obtaining polarization data in the field equal i : an 1 
When obtaining polarization data i ld « ie « * ot" & onoee J 
increments of polarizing current are applied at ° ° 
. : 
approximately equal time intervals (1 or 2 min ' — 
intervals). —8 * : 
2.3. Removal of Corrosion Products E } 
. P . 2 r Re) 
After 483 davs of underground exposure, the q 
. 3 . . oo 
aluminum and steel specimens were removed for eo OFF 


: : ‘ : , 56 Fe fe) ° 4 10-22-59 “| 
cleaning, weight loss determinations and _ pit-depth * tee, 

measurements. The aluminum specimen was sub- 
merged in hot water and brushed with a stiff-bristle 


brush. After the rubber caps were taken off, the r 
wire and solder were removed by applying heat a) 7] 
from a soldermg iron. The specimen was then | g ' or ; 
: Ate? : : os ; Q- 
placed in concentrated nitric acid for 5 min, washed G6 > © — 4 
: : : Z . as 3 e se 1-26-60 
: a ° 1 ° rTN . oO *-e a 
in running water, dried, and weighed. The effee- | j } “n 4 
. “. . - > 
tiveness of the cleaning procedure was checked by | © —g b ies - 
repeating the nitric acid treatment and reweighing. | $ q : a J 
There was no significant difference between the two | >= _4 Lb | 
weighings, indicating that the specimen was free of = o® 
° 7 ox 7 
corrosion products. z S > ODP” 
: w—6 fe . 3-92-60 + 
Rust which had formed on the steel specimen was : laa * ee, 
very adherent. First, the steel was cleaned cathod- | &@ ; tee, 
ically in a 3 percent NaCl solution for 5 hr at a oe 7 
current density of 1 amp/ft.2. The specimen was q 1 
then immersed for 2 hr in a 10 percent ammonium 345 oo” 4 
citrate solution (at 150 °F) made alkaline with . 200 4 
; ‘ a | e ° eo—g-—2 Sera I 6-15-60 
ammonium hydroxide. Finally the steel was hand ss S a “dl 


brushed (wire brush), washed. dried, and weighed. 


3. Results and Discussion 











Some of the polarization data covering the 16 <4 : ia 7 
months of exposure are shown for steel and alumi- a viata err, B-I-60 
num, respectively, in figures 1 and 2. The dates “ih re | 
are given so that the data can be related to the time & 4 
of year. The initial curves were obtained after the ~g + _] 
specimens were in the ground for 3 days and the ' | 
final curves shown were run on the day preceding q id 
their removal. All data were taken from the ad . ad “| 
recorder charts and plotted on semi-leg scale. The re +o——6- a “ied 12-5-60 
corrosion current controls the magnitude of the 36 i ~] 
polarizing currents 7p and -/q [i], indicated by the r ". 
breaks in the curves (fig. 1 and 2). It will be =8 “| 
observed in figures 1 and 2, that the current Jq | a a 
approaches in value the current Ip with the passing | 2 a 622845 2 4 6 810 20 40 


of time, particularly for the steel. Thus, the 


' ; APPLIED CURRENT, ma 
tendency is a change from cathodic to mixed control. 


- ° . ‘ 7 1] l¢ 2 1 CUrVES ? i 8 
The polarizing currents, Jp and Jq, used in ealeu- Ficure 1. Typical polarization curves of the underground steel 
. ‘ . . . ° specimen during the 16 month exposure pe riod.@ Cathodic, 
lating the values of corrosion current were obtained {nodic 
as shown in figures 1 and 2. Polarization rates a 
= : - lo convert applied current to current density (ma/ft?), multiply by 2.5 
(AV/AZ) were calculated at the slightly larger values 
of current indicated by arrows on the curves. These | the present data are shown in figure 3. The present 


rates seem to bear a logarithmic relation to the cur- | data for aluminum also appear to fit’ reasonably 
rents Jp as shown by previous work [2]. Data from | well. Thus, knowing the polarization rate (AV/AJ) 
the previous work on steel in a NaCl solution and | it is possible to estimate the current Jp from the 
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Figure 2. Typical polarization curves of the underground 
aluminum specimen during the 16 month exposure period. 


r of ‘athodic, 


ro convert applied current to current density (ma/ft?) 


Anodic. 
, multiply by 2.5. 
curve (fig. 3). The empirical equation for data 
(fig. 3) may be written: 
Ip=c R or log Ip=log e—n log R., where 
c=100 (y—intercept), R=AV/AT in mv/ma/ft’, 
n 1.1 (slope of curve) and Jp is expressed in 
ma/ft.? 
This relationship would be useful where the break 
in the curve is not readily apparent. The author 
is currently investigating the scope of this relation- 
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FicureE 3. Logarithmic relationship (cathodic polarization), for 
steel and aluminum, between the polarizing current, Ip, and 
the polarization rate, AV/Al. 

ABare steel underground (covering 16 months exposure); A bare aluminum 
underground (covering 16 months exposure); @ bare steel in 3 percent sodium 
chloride solution (covering 14 days exposure); coated steel (carbonate coating) 
in 3 percent sodium chloride solution (covering 29 days exposure); @ coated steel 
(carbonate, progressively damaged) in 3 percent sodium chloride solution (cover- 
ing 36 days exposure 


ship with possible application to ferrous alloys such 
as stainless steels. When the type of control is 
predominantly anodic, the anodic polarization data 
should fit the relationship. 

All data obtained for the steel and aluminum from 
the polarization curves (semi-log scale) used for 
calculating weight losses are shown respectively in 
tables 1 and 2. The calculated and actual weight 
loss for steel agree within 6 percent. On a percentage 
basis the same cannot be said for the aluminum but, 
considering the relatively small weight loss involved, 
the calculated weight loss is considered quite signifi- 
cant. Actually, the difference (about 50 mg) between 
the actual (as weighed) and calculated corrosion 
losses could be largely cleaning error for a soft metal 
such as aluminum. 

For the first 3 months of exposure, polarization 
data were obtained more frequently than later, 
During this initial period, changes in the corrosion 
rate of both steel and aluminum appeared to be 
occurring frequently. Also, the corrosion rates were 
higher and seemed to be influenced more by changes 
in the weather and the passing of time. The effects 
of time, temperature, and rainfall on the corrosion 
current for steel and aluminum in this particular 
underground environment are shown in figure 4. 
Potentials of both the steel and alu ninum are also 
plotted with soil resistance and climatological data. 
For the steel, the corrosion current (%) varied con- 
siderably for the first 2 months of exposure. The 
rainfall early in September 1959 and again early in 
October lowered the soil resistance and apparently 
caused the two peaks (fig. 4) in the corrosion current. 
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TABLE 1. Weight loss calculated from polarization data on a 
steel specimen erposed unde rground for 16 months 





Polarizing current ® Weight loss 
at break in curve Corro- | Exposure 
Time of sion time 
year current Calculated Actual 
Cathodic Anodic i umulative 
Ip I 
m me ma Days mg 
&-—12-59 2 10.0 2.0 1 A) 
14 2.3 9.2 1.8 3 14 
19 2.4 1.8 1.¢ S 357 
25 2.3 4.6 1. § 14 2 
94-59 6.6 13.2 1.4 24 1332 
x 5.9 R.4 3. 5 28 1712 
1s 3.5 4.2 1.9 3S 2387 
20 30 5.0 1.9 44 2009 
10-1-59 7.9 3. 5 5.0 51 3084 
9 4.1 6.7 9 5 ] 324 
22 4.4 6.4 2.¢ 72 4609 
11-4-59 2.4 7 2.3 RS h4n2 
13 2:3 7.7 La 94 5O7S 
30 - 6.5 1.Y lil 6828 
12-22-59 2.{ 1.9 1.4 133 7708 
1-13-it 1.9 2.8 1.1 l 8423 
i 1.4 1.7 0.77 168 872. 
2-19-60 $ 1.4 0. 67 1¥ Q1LA7 
34440 1.0 1.5 0. 60 206 Y381 
Y Lo 1.3 0.62 211 9457 
2 1.7 1.9 0.90 227 9761 
1-0 2.2 2.3 1.1 234 10061 
28 3.7 2.7 13 261 10721 
20-90 2 2 1.3 283 11435 
23 3.0 ( 1 286 11540 
6-2-0 1.8 206 11940 
15 3.2 1.0 1.8 309 12526 
30) 2.7 3.0 1.4 324 13088 
7-18-60 o-& 2.6 1.2 342 13674 
26 2.3 2.¢ 1.2 350 131g 
Qn 1.9 1.8 0.92 352 13969 
8-140 1.9 7 0.90 356 14060 
19 3.6 3.4 1.7 374 14645 
9-1-60 2. 2 2.3 1.1 387 15100 
14 3.0 3.1 1.5 400 15522 
10-11-60 2.9 2.¢ 1.4 427 16467 
26 2.@ 2.3 1.1 442 16954 
11-23-60 2.0 2.1 1.0 470 17723 
12-5-60 1.6 2.0 0.89 4182 18005 
6 483 18027 19055 
4 From curves in figure 1 and other curves not shown, 
bj Ip Iq Ip+ Iq). 
e Weight loss (g K ti, where K=2.8938X10-4 g per coulomb, i=i A verage 
current (amp) for the period (¢ in seconds) between successive readings. The 


value of i, at the instant of exposure and at the end of exposure is taken as the in 
itial and final values, respectively, as calculated 

i No anodic run; value based on other Jp/Iq ratios around that time 

e No apparent break; value based inversely on relative cathodic and anodic 
polarization rates AV/AZ and Ip 


During the next 5 months, the corrosion current con- 
tinued to decrease, presumably due to the drop in 
air temperature and resultant lowering of soil tem- 
perature. The lowest value of corrosion current was 
observed in March 1960, the coldest time of the 
vear when the soil resistance was maximum. From 
this time until the beginning of June, air tempera- 
tures continued to rise, soil resistance became lower, 
and the corrosion current again increased. During 
the summer 1960, fluctuations in the corrosion cur- 
rent apparently were caused by rainfall. With the 
coming of colder weather, the soil resistance again 
increased and the corrosion current diminished. The 
major fluctuations in potential of the steel occurred 


“ABLE 2 


alun 


Wetght loss calculated from polarization data on an 


minum specimen exposed underground for 16 months 


Polar 1 Weight loss 
tl i Corro Exposurt 
e of slo! time 
urrent Calculated Actua! 
Cathe Al cumulative 
| ] 
wa Days mu 

8-12-59 0 3 0.10 1 Os 

14 l ; 1] ; 2 6 

19 12 2 OSI s 6.4 

2 iT x 067 14 10.0 
~ 4-59 11 O71 24 15. ¢ 

is { x O28 3s 21.1 

2u (4 OZ 10 23. 4 
( 1-59 039 ( (V2 | 23.8 

) ( ( Teal ; 2 

2 03 is Olu 27.2 

4-59 030 ( O19 s 20 3 

; 024 j ( a4 30 

$4 03. Hs O02 111 $3. 1 
12-22-50 O34 y ( 133 37. 2 
l i () ( ] 10), ¢ 

2 01 13 009 168 $1.9 
2-19-41 ( } 2 19 13.8 

1400 Ol. { awe 1M) 5.0 

9 01 03 00S 211 15.3 

25 O14 033 on 227 4}. 5 
1- 6-60 Ole Ont » O34 17.6 

2s Olf OA Oll 261 19.7 
5-20-60 O14 042 O11 283 51.7 

23 O11 $4 009 2St} 1.9 
15-60 012 2 008 $004 53. 6 

30) 009 a 006 $24 4.4 
7-18-60 OO 01 OO. 342 55. 1 

26 OOS ( OO. 3H 55.4 
S- 1-60 009 12 in $e 55. 6 

19 O11 O13 Oot 74 AO. 5 

1-60 OO O11 rT 387 57.0 

14 O11 O12 006 400) 57.6 
10-11-60 007 O11 O04 127 5S. 7 

# 006 mH O03 442 AY. 1 
11-23-00 OO7 TI O04 i70 60.0 
2 »-H0 O06 (nit 003 {2 60.4 

+) iS 60.7 109 
From curves in figure 2 and other curves not shown 

Ip Iq/Ip+Iq 
Same as footnote it t ept A=0.9316X 10-4 g per coulomb 
No anodic run; value based on other Ipil tic iround that time 


during the early months of exposure. After the 
corrosion products began to build up, there was a 
tendency for the potential to become gradually more 
electropositive. 

The corrosion rate of the aluminum seemed to be 
less affected by changes in weather than that of the 
The current for the aluminum (fig. 4) dropped 
considerably during the first 2 months and gradually 
thereafter. In contrast with the steel, the aluminum 
corrosion current decreased about 30-fold while that 
of the steel decreased around 8-fold during the 16 
months. 

After the specimens had been cleaned, the pit 
depths were also measured. On the steel, the maxi- 
mum pit depth was 0.035 in. and the six deepest 
pits averaged 0.03 in. For the aluminum, the 
maximum pit depth was 0.008 in. The six deepest 
pits averaged 0.007 in 


steel. 
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Figure 4. Effect of time and weather on the corrosion rate and potential of steel and aluminum 
exposed underground for 16 months 


A Potential, steel corrosion current, steel; @ air temperature; & potential, aluminum; @ corrosion current, alumi 
num; ™ soil resistance, between steel specimen and reference electrode; | precipitation (water) in area—U.S. Weather 
Bureau 
4. Summary discussed in a previous publication [2], which is 


helpful in anticipating rates of corrosion. 

Both materials corroded at a faster rate during 
the early months of exposure than later. The cor- 
rosion rate of the aluminum diminished more 
rapidly and to a greater degree than that of the steel. 
The corrosion rate of the steel was affected more by 
temperature and rainfall than was that of the 
aluminum. 

This investigation was not intended to evaluate 
the relative merits of the metals exposed to corrosion. 


A steel and an aluminum specimen were exposed 
underground for 16 months to a soil having an 
average resistivity of about 7500 ohm-em. During 
this time, eathodie and anodic polarization curves 
were automatically recorded on 40 occasions with 
the idea of evaluating a technique for measuring 
rates of corrosion and noting the effeets of the 
weather on such rates. 

The weight loss of each material caused by cor- 
rosion agreed quite well with the calculated losses taps . ‘ 2 : : 
based on the polarization data. The cathodic While their relative behavior is of interest, it must 
polarization data for both the steel and aluminum | be remembered that only one soil environment was 
seem to fit reasonably well a logarithmic relationship, | involved. 
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The weather data were copied from the local 


climatological data sheets, U.S. Department of Com- 
merce, Weather Bureau. 
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On the spectrum of terrestrial radio noise at extremely low 
frequencies, H. R. Raemer. ./. Research NBS 65D (Radio 
Propagation) No. 6, (November— December 1961) 70 cents. 

A theory of the frequency spectrum of radio noise at extremely 
low frequency (ELF) is presented and the results compared 
with recent measurements of the first five “Schumann’’ reso- 
nant modes (between 8 and 34 ¢/s) made by Balser and 
Wagner [1960]. The source of this noise is assumed to be 
return strokes in vertical cloud-ground lightning flashes dis- 
tributed randomly in time, uniformly in angular displacement 
along the earth relative to the observer, and with statistics 
of stroke duration, inter-stroke intervals, and strokes per 
flash taken from studies of thunderstorms reported by J. C. 
Williams. Thus, the mathematical model for the noise 
sources is an extremely simple one, being analogous to the 
shot effect in electron devices. The electromagnetic model 
employs the familiar waveguide mode theory, assumes a 
sharply bounded homogeneous ionosphere, and neglects the 
earth’s magnetic field 

Agreement between the shape of the theoretical and observed 
spectrum is good for the first three modes and rather poor 
for the higher modes. 

It is found by matching the theoretical resonant frequencies 
to the observed resonances that the product of effective 
ionosphere height A and the square root of effective con- 
ductivity yo; is a decreasing function of frequency. The 
functional dependence of this quantity on frequency is deter- 
mined and used in the calculation of the mode spectrum 
Discrepancies between the theory and experimental results 
are believed to be partially due to the artificiality of the 
sharply bounded homogeneous ionosphere model and to fail- 
ure to give sufficient probability weighting to equatorial 
regions of abnormally high thunderstorm activity. 

These last items are the subjects of continuing work on the 
extension of the theory ‘ 
Bibliography of temperature measurement 
June 1960, C. Halpern and R. J. 
(1961) 14 cents 

There are presented more than 500 references to the field of 
temperature measurement. These references were collected 
from two general sources: Scientific and technical literature, 
and government reports. The period covered is from 1953 
to June 1960, with some from earlier dates. For convenience 
of the user, the references are divided inte a number of 
categories based on the type of instrument used. Some 
references to calibration of instruments and to scientific 
theories, on which temperature measurement is based, are 
also presented 

Corrected optical pyrometer readings, 1). I. Poland, J. W. 
Green, and J. L. Margrave, NBS Mono. 30 (1961) 54 cents. 
The table of corrected optical pyrometer readings enables 
optical pyrometer users to convert observed temperature 
immediately to the true temperature, if the effective emis- 
sivity of the material being observed is known. The table 
gives observed temperatures from 1,000 to 3,000 °K in in- 
crements of 5 degrees, from 3,000 to 5,000 -K in increments 
of 10 degrees, from 5,000 to 10,000 °K in inerements of 50 
degrees, and from 10,000 to 39,900 °K in increments of 100 
degrees. For these, true temperatures are tabulated for 49 
emissivities ranging from 0.02 to 0.98 in increments of 0.02. 
These calculations were made on a 650 electronic computer 
using Planck’s law, and the value C,=1.438 em deg for the 
radiation constant. 

Publications of the National Bureau of Standards July 1, 
1957, to June 30, 1960 (Includes titles of papers published in 
outside journals 1950 to 1959), B. L. Arnold, NBS Mise. 
Publ. 240 (1961) $2.24 

This Supplement lists the publications of the National Bureau 
of Standards from July 1, 1957, to June 30, 1960. A _ brief 
abstract for each publication is given. It also includes a list 
of those papers by Bureau authors that have appeared in non- 
National Bureau of Standards Journals from 1950 through 
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1959. Subject and author indexes are included, and general 
publication information is furnished. This catalog and pre- 
vious lists (see page 2) give reference to the publications 
printed by the Bureau during the period 1901 to June 30, 
1957. 

NBS—-source of American standards, W. A. Wildhack, /SA 
J. 8, No. 2, 45-50 (February 1961). 

A general description of the National Bureau of Standards 
and its program is given, with particular emphasis on devel- 
opment of standards and provision of calibration services. 
Calculation of properties of magnetic deflection systems, S. 
Penner, Rev. Sci. Instr. 32, No. 2, 150-160 (February 1961). 
A convenient matrix method for calculating properties of 
magnetic deflection systems is presented. This method is 
applicable to particle beams of small spatial and angular 
extent, and small energy spread. Equations for quadrupole 
lenses and for deflecting magnets are given. Examples are 
given to show the procedure for calculating the parameters 
of magnet systems. 

Stepless variable resistor for high currents, C. R. Yokley and 
J. B. Shumaker, Jr., Rev. Sci. Instr. 32, No. 1, 6-8 (January 
1961). 

A variable-area water rheostat has been applied to the prob- 
lem of current control in electric arc and plasma generating 
equipment. Resistance between 0.5 and 30 ohms can be 
obtained by adjusting the height of circulating water which 
surrounds pairs of fixed vertical aluminum plates. The elec- 
trolyte is ordinary city, water which is circulated through the 
tank and discharged. No special filtering, pumping, or other 
treatment is used. The unit is located in the open air out- 
side the laboratory and all control functions are operated 
from a remote switching panel. The unit has been operated 
using a maximum of 125 amp at 400 v de and 1000 amp at 
230 v de. 

Physical quantities proposed for radiation measurements, 
L. S. Taylor, Proc. [Xth Intern. Congress of Radiology, July 
23—30, 1959 (Miinchen, Germany), 1301—1302 (1960). 

This is the introductory statement for a symposium on radi- 
ation measurements. -It was designed to outline the back- 
ground in one phase of the problem of radiation quantities 
units. 
Improving rectifier circuits, G. F. 
p. 86-87 (April 7, 1961). 

The conventional half-wave rectifier and voltage doubler 
circuits can be modified to use an inductor-input filter ad- 
vantageously. The results of these modifications are reduced 
peak rectifier current and improved load-voltage regulation 
and smoothing. The first modification requires an additional 
rectifier; the second requires a two-winding input inductor. 
The robustness of life testing procedures derived from the 
exponential distribution, M. Zelen and M. C. Dannemiller, 
Technometrics 3, No. 1, 29-49 (February 1961). 

Almost all the statistical procedures in current use for eval- 
uating the reliability of components or equipment rest on the 
assumption that the failure times follow the exponential 
distribution. However, in practical situations one rarely has 
enough data to determine whether failure times are actually 
exponential. This paper studies the behavior of several 
statistical life testing procedures based on the exponential 
failure law if. the true failure law is the Weibull distribution. 
It is found that these statistical techniques, which are widely 
used, are very sensitive to departures from initial assumptions. 
Applying these techniques to life test data when the exponen- 
tial failure law is not satisfied may result in substantially 
increasing the probability of accepting components or 
equipments having poor mean-time-to-failure. 

How to evaluate accuracy, W. J. Youden, Materials Research 
and Standards (Am. Soc. Testing Materials, Philadelphia, Pa.) 
1, No. 4, 268-271 (April 1961). 

This paper is one of a series of four integrated papers on errors 
of measurement under preparation by four authors for the 
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June 1960 ASTM meeting. It presents a logical breakdown 
of the error in a measurement into (a) the systematic error 
inherent in the procedure, (b) the local systematic error of 
the laboratory using the procedure and (c) the random error 
(precision). This breakdown should facilitate efforts 
attain better accuracy. Several methods are given for identi- 
fying sources of error in measurements. . 

Improved standard for the calibration of vibration pickups, 
R. R. Bouche, Exp. Mech., 1-6 (April 1961 

As part of a continuing program on the development of 
primary vibration standards, an electrodynamie standard was 
modified so that it is free of transverse motion throughout 
almost all of its frequency range. In the vibration-pickup 
calibration service recently established at the National Bureau 
of Standards, the standard is used at frequencies from 10 cps to 
2,000 cps at accelerations up to 10 g. 

Velocity and acceleration pickups in common use calibrated 
on the standard provide reasonably accurate secondary 
vibration standards for use by government and industrial 
research laboratories. However an assembly of a small coil 
used with a magnet has been found to be a accurate 
secondary vibration standard. 

The development of more stable gage blocks, \l. R. Meyerson, 


to 


more 


T. R. Young, and W. R. Ney, ASTM Bull. 1, No. 5, 368-374 
(May 1961). 
American precision tool manufacturers have requested a 


higher degree of accuracy in gage blocksjand their calibration 
than is currently available in order to keep pace with the needs 
of the missile, aircraft, and other precision industries. To 
meet this need, a long-range research program aimed at an 
accuracy of one part in 10 million in gage-block length calibra- 
tions is under way at the National Bureau of Standards. 
Significant advances have been made recently with the 
development of highly stable blocks and ultraprecise measure- 
ment techniques. This paper describes the highlights and 
presents the data obtained in the initial phases of the pro- 
gram. 

The power spectrum and its importance in precise frequency 
measurements, J. A. Barnes and R. C. Mockler, JRE Trans 
Instrumentation 1-9, No. 2, 149-155 (September 1960). 

The power spectral density functions of a frequency multiplier 
chain, driven by several different crystal oscillators, were 
obtained by comparing the output with a second chain which 
was stabilized with an ammonia maser. The frequency of 
the maser stabilized chain was demonstrated to be relatively 
fixed; the power spectrum of the other chain was determined 
by two different methods. The results are compared.  Pos- 
sible errors and uncertainties introduced by the methods are 
discussed. An analysis is made that relates the instantaneous 
frequency fluctuations of a signal with the power spectral 
density function 

Analysis predicts that when frequency modulation occurs in 
the first stages of frequency multiplication or in the primary 
frequency oscillator, the output power spectrum is, in general, 
not symmetrical. Furthermore, the sidebands are increased 
in intensity by the multiplication process. This is, in fact, 
observed to be the case. It is shown that a frequency counter 
will measure the frequency of the center of gravity of the 
power spectrum. 

If signals having a complex power spectrum are used in precise 
frequency measurements, errors may result. 

The nation’s electronic standards program: where do we now 
stand? H. Lance, [RE Trans. Instrumentation I-9, Vo. 
- 94 100 (Se ple mber 1960). 

Recently increased standards activities throughout the nation 
have resulted in a large increase in the calibration services 
requested from NBS and in a number of related problems 
Several of these are discussed, including desirable criteria for 
the staffing and instrumentation of electronic standards 
laboratories, the choice of interlaboratory standards to be 
submitted to NBS, and the use of these standards before and 
after calibration. An attempt is made to clarify the concept 
of “traceability of calibrations to N BS,’’ and conditions under 
which a calibration may be called traceable are discussed. 
Suggestions made regarding the use of available calibra- 
tion services and the action to be taken when needed calibra- 
tion services art available. The need for an 
of standards laboratories is pointed out, and many areas are 
noted in which I could make valuable 
contributions. 
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Trends in the technology of automatic data processing, 
S. N. Alexander, AMA Report 48, 38-42 (1960). 

This is a summary of a talk presented to the Annual Office 
Management Conference of the American Management 
Association, Inc., in whichthe four elementary functions of 
automatic data systems are discussed as they 
relate to management. The progressive development of input 
and output devices and data storage and memory systems as 
well as improved logical organization of the central processor 
described. In conclusion, the point is made that an evolu- 
tion of appropriate automatic program procedures now in 
process is of comparable importance 

Reproducibilities of carbon and germanium thermometers at 
4.2 °K, M. H. Edlow and H. H. Plumb, Proc. Cyrogenic Eng. 
( onf., Boulder, Colo., Aug. 1960. 6, 542 (1961). 

Some new results are presented on the characteristics of 
secondary thermometers at temperatures in the liquid helium 
region. These results are for the reproducibility of the resist- 
ances of carbon and germanium resistors at 4.2 °K when they 
were eycled from 300 °K to 4.2 °K. The resistances of these 
thermometers were determined functions of the vapor 
pressures of liquid helium in a liquid helium bath at nearly 
constant temperature The of Allen-Bradley 
resistors (10°, 4 watt which have been continuously in 
liquid helium for many weeks at constant temperature have 
been found to drift. Similar drifts have not been detected for 
germanium resistors. 


Other NBS Publications 


Journal of Research 65A 
August 1961). 75 cents. 
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Phys. and Chem.) No. 4 (July 


Electrical properties and kinetics of electrode reactions, R. J. 


srodd. 
Effect of hydrostatic pressure upon the relaxation of bire- 
frigence in amorphous solids, R. M. Waxler and L. H 


Adams. 

Vapor pressures of platinum, iridium, and rhodium, R. F 
Hampson, Jr., and R. F. Walker 

Crystallization of bulk polymers with chain folding: 
of growth of lamellar spherulites, J. D 
Lauitzen, Jr. 

Phase equilibrium relations in 
oxidersniobium pentoxide, R. 8S. Roth and J. L. Waring. 

Solid state reactions involving oxides of trivalent 
S. J. Schneider, R. S. Roth, and J. L. Waring. 

Gamma irradiation of fluorocarbon polymers, R. E. 
and L. A. Wall 

Inhibition of diffusion flames of methyl bromide and trifluoro- 
methyl bromide applied to the fuel and oxygen sides of the 
reaction zone, E. C. Creitz. 
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cations, 


Florin 
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ber-October 1961 


Phys. and Chem. 
75 cents. 
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of a monitor for use in bremsstrahlung beams, 
E. G. Fuller and E. Hayward 
Mass spectrometric study of NF», NF:, NoFs, and NoF,, 
J.T. Herron and V. H. Dibeler 
tate of the reaction NO+N, and some heterogeneous reac- 
tions observed in the ion source of a mass spectrometer, 


J s. Herron 
Synthesis of the humites nMg2S8iO,-Mg(F, OH)., A. Van Val- 
kenburg 
1ase equilibria in systems involving the rare earth oxides 
Part III. The Eu,.0,;-In,O; system, 8. J. Schneider. 
Heats of hydrolysis and formation of dimethoxychloroborane, 
M. V. Kildav, W. H. Johnson, and E. J. Prosen 
Tritium-labeled compounds VII Isotope effects in the oxida- 
tion of p-mannitols-C’ and pb-mannitols-t to p-fructoses, 
L. T. Sniegoski, H. L. Frush, and H. S. Isbell. 
Franck-Condon factors to high vibrational quantum numbers 
I: No and Nt, R. W. Nicholls 
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Note on the ‘‘baffled piston” 
Some results on non-negative 
and B. N. Moyls 


Probability inequalities of the Tehebycheff type, I. 


problem, F. Oberhettinger. 
matrices, M. Marcus, H. Mine, 


R. Savage 


Journal of Research 65D (Radio Prop.) No. 5 (September 
October 1961 75 cents. 

Frequency dependence of D-region scattering at VHF, J. C 
Biair, R. N. Davis, Jr., and R. C. Kirby. 

Theoretical scattering coefficient for near vertical incidence 
from contour maps, H. 8. Havre and R. K. Moore 

Mutual interference between surface and satellite communi- 
cation systems, W. J. Hartman and M. T. Decker. 

VHF and UHF signal characteristics observed on a long 
knife-edge diffraction path, A. P. Barsis and R. 8. Kirby 

Iixperimental study of inverted L-, T-, and related transmis- 
sion-line antennas, S. Prasad and R. W. P. King. 

Reflection from a sharply bounded ionosphere for VLF 
propagation perpendicular to the magnetic meridian, D 
1) Crombie 

Resonance of the 
Poeverlein 

Observed attenuation 
radio waves, A. G. Jean, 
D. F. Wasmundt 

A note concerning 
waves, J. R. Wait. 

Computation of whistler ray paths, I. Yabroff. 

On the analysis of LF ionospheric radio propagation phe- 

J. R. Johler 


space between earth and ionosphere, H 
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Journal of Research 65D (Radio Prop.) No. 6 (November 
December 1961). 75 cents. 


The solar wind, E. N. Parker. 


Attentuation coefficients for propagation at very low fre- 


quencies (VLF) during a sudden ionospheric disturbance 
SID), E. T. Pierce 
Dipole radiation in a conducting half space, R. K. Moore 


and W. EF. Blair. 
Reliability of atmospheric 
Herman 
kffeets of the 
Blackband 
On the spectrum of terrestrial radio noise at extremely low 
frequencies, H. R. Raemer. See above abstract) 
The nonsingular embedding of transition processes within a 
more general framework of coupled variables, J. Heading 
Worldwide VLF standard frequency and time signal broad- 
casting, A. D. Watt, R. W. Plush, W. W. Brown, and A. H 


radio noise predictions, J. R 


ionosphere on VLF navigational aids, W. T. 


Morgan 
Design of panoramic ionospheric resources, L. H. Heisler and 
L.. D Wilson 


\ quick method for estimating the stage of the sunspot evele. 
W. B. Chadwick 
Measurements of low-angle radiation from a monopole, A. C. 


Wilson 


Ideal thermodynamic functions and isotope exchange 
functions for the diatomic hydrides, deuterides, and tritides, 
L. Haar, A. S. Friedman, and C. W. Beckett, NBS Mono. 
20) (1961) $2.75. 

Causes of variation in chemical analysis and physical tests of 
portland cement, B. L. Bean and J. R. Dise, NBS Mono. 28 
1961 

Thermal 


25 cents 

expansion of technical solids at low temperatures 
A compilation from the literature, R. J. Corruecini and 
J. J. Gniewek, NBS Mono. 29 (1961) 20 cents 

Mean electron density variations of the quiet ionosphere 
Mav 3. 1959, J. W. Wright, L. R. Weseott, and D. J. 
Brown, NBS TN40-3 (PB151399-3) (1961) $1.50. 

Bibliography on ionospheric propagation of radio 
1923-1960), W. Nupen, NBS TN84 (PB161585 
$7.00 

The absorption spectra of magnesium and manganese atoms 
in solid rare matrices, O. Schnepp. J. Phys 
Solids 17, Nos. 3/4, 188-195 (1961 
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Runge band system, R. W 

1705-1711 (1960). 

Ionospheric mapping by 
and R. M. Gallet, 
1960 

The activation energy for hydroge n atom addition to propyl- 


Phys. 38, 


methods, W. B. Jones 
12, 260-264 (December 


numerical 
Telecommun. J 


ene, M. D. Scheer and R. Klein, J. Phys. Chem. 65, 
No. 375, 375-377 (1961 

Intramolecular rearrangements. I. sec-butyl acetate and 
sec-butyl formate, R. Borkowski and P. Ausloos, J. Am. 
Chem. Soc. 83, No. 5. 1053-1056 (March 1961 


Titanium (1) chloride and titanium (11) bromide (titanium 
trichloride and titanium tribromide), J. M. Sherfey, 
Chapter IVB, Sec. 17, p. 57-61, Book, Inorganie Syntheses, 
by E. G. Rochow (MeGraw-Hill New York, 
N.Y., 1960) 

Science news writing, D. M. Gates and J. M 
133, No. 3447, 211-214 (January 1961 
Integrated starlight over the sky, F. E. Roach and L. R. 
Megill, Astrophys, J. 138, No. 1, 228-242 (January 1961). 
Elastomers for static seals at eryogenic temperatures, D. H. 
Weitzel, R. F. Robbins, G. R. Bopp, and W. R. Bjorklund, 
tev. Sei. Insts. 31, No. 12, 1350-1351 (December 1960). 

The electromagnetic fields of a dipole in the presence of a 
thin plasma sheet, J. R. Wait, Appl. Sci. Research 8, Sec. 
B, 397-417 (1960 

New standards for the space age, A. T. McPherson, 
nauti¢s 6, No. 1, 24-25; 50-54 (January 1961 

Near infrared atmospheric transmission to solar radiation, 


S00k Co., 


Parker, Science 


Astro- 


D. M. Gates, J. Opt. Soc. Am. 50, No. 12, 1299-1304 
December 1960 
Determination of crystallite size distributions from X-ray 


line broadening, A. Bienenstock, J 
187-189 (February 1961 

Correlation of visual and subvisual auroras with changes in 
the outer Van Allen radiation zone, B. J. O’Brien, J. A. 
Van Allen, F. E. Roach, and C. W. Gartlein, IGY Bull. 
No. 45, 1-16 (March 1961 

Are source for high temperature gas studies, J. B. Shumaker, 
Jr., Rev Instr. 32, No. 1, 65-67 (January 1961). 

The dependence of the melting temperature of bulk homo- 
polymers on the crystallization temperature, L. Mandel- 


Appl. Phys. 32, No. 2, 


Sc. 


kern, J. Polymer Sci. XLVIT, issue 149, 494-496 (Julv 
1960). 
Spiral patterns in geophysies, V. Agy, J. Atmospheric and 
Terrest. Phys. 19, 136-140 (1960) 
Iron (99.9+-), G. A. Moore and T. R. Shives, Metals Handb. 
1, 1206-1212 (1961) 


The interfacial properties of polyesters at 
surfaces, R. R. Stromberg, 
October 1959) 

Ionospherie absorption at 
pulsations, W. H. Campbell and H. Leinback, J. 
Research 66, No. 1, 25-34 (January 1961 

Some properties of new or modified excitation sources, M. 
Margoshes, Am. Soc. Testing Materials, Spec. Tech. 
Publ. No. 259, 46-58 (1959). 

Seasonal and day-to-day changes of the central position of 
the S, overhead current system, 8. Matsushita, J. Geophys. 
Research 65, No. 11, 3835-3839 (November 1960). 

Hydrogen sulfide precipitation of the elements from 0.2—0,5 
normal hydrochloric acid, J. I. Hoffman, Chemist-Analyst 
50, No. 1, 30 and 32 (March 1961 

The mechanism of electrolytic deposition of titanium from 
fused salt media, W. E. Reid, J. Electrochem. 108, 
No. 4, 393-394 (April 1961 

Studies of elevated temperature corrosion of type 310 stainless 
steel by vanadium compounds, H. L. Logan, Corrosion 17, 
109-111 (April 1961 

Deep penetration of radiation, U. Fano and M. J. 
Proc. Symp. Appl. Math. XI, 43-59 (1961 

\ simple low-temperature specimen holder for an 
diffractometer, D K. Smith, Norelco Reporter 
No. J, 11-12 (January-February 1961 

Microbalance techniques for high temperature application, 


and water 
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Soc. 


Jerger, 
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R. F. Walker, Book, Vacuum Microbalance Techniques, 
edited by M. J. Katz. I, 87-110 (Plenum Press Ine. 
New York, N.Y 1961 


Investigation of bond in beam and pull-out specimens with 
high-yield-strength deformed bars, R. G. Mathey and D. 























































Watstein, J. 

(March 1961). 

Vacuum ultraviolet photolysis of ethane: Molecular detach- 
ment of hydrogen, H. Okabe and J. R. MeNesby, J. Chem. 
Phys. 34, No. 2, 668-669 (February 1961). 

Shielding calculations for civil defense, C. Eisenhauer, Health 
Phys. 4, No. 2, 129-132 (1960). 

Hydrogen atom reactions with propene at 77 °K. Dispropor- 
tionation and recombination, R. Klein and M. D. Scheer, 
J. Phys. Chem. 65, 324-325 (1961). 

The effect of solvents on the y-ray 
acetate and acetone, P. Ausloos, J. 
No. 5, 1056-1060 (March 1961). 

On the nature of the crystal field approximation, C. M. 
Herzfeld and H. Goldberg, J. Chem. Phys. 34, No. 2, 
643-651 (February 1961 

A study of auroral coruscations, W. H. Campbell and M. H. 
Rees, J. Geophys. Research 66, No. 1, 41-55 (January 1961). 

The nature, cause and effect of the porosity in electrodeposits. 
VII. A microscopic examination of nickel-chromium coat- 
ings after atmospheric corrosion, D. W. Ernst and F. 
Ogburn, Plating 48, No. 4, 491-497 (May 1961). 

Fracture characteristics of notched tensile specimens of tita- 
nium and a titanium alloy, G. W. Geil and N. L. Carwile, 
Materials Research and Standards (Am. Soc. Testing Ma- 
terials, Philadelphia, Pa.) 1, No. 1, 16-21 (Jan. 1961). 

Absence of an isotope effect in the fractional recrystallization 
of alpha-p-glucose-1-t, H. 8. Isbell, H. L. Frush, and N. B. 
Holt, Anal. Chem. 33, 225-226 (Feb. 1961). 

Applications of graphs and Boolean matrices to computer 
programming, R. B. Marimont, SIAM Rev.. 2, No. 4, 
259-268 (Oct. 1960 

Comparison of national standards for roentgen measurement, 
H. O. Wyckoff, IXth Intern. Congress of Radiology, July 
23-30, 1959 (Miinchen, Germany), 1315-1318 (1960). 

Environmental factors in a family-size underground shelter, 
R. P. Achenbach, Proc. Meeting on Environmental Engi- 
neering in Protective Shelters, Feb. 8-10, 1960, Natl. Acad. 
Sci.-Natl. Research Council, Wash., D.C., 69-118 (1961). 

Absorption and scattering of photons by holmium and er- 
bium, E. G. Fuller and E. Hayward, Proc. Intern. Conf. 
Nuclear Structure, 760-763 (Kingston, Canada, 1960). 

Introductory remarks, U. Fano, Conf. Coherence Properties 
of Electromagnetic Radiation, AFOSR-583, Report on 
Conf., 31-38 (University of Rochester, The Institute of 
Opties, Rochester, N.Y., April 1961). 

Available heat sinks for protected underground installations, 
B. A. Peavy, Proc. Meeting on Environmental Engineering 
in Protective Shelters, Feb. 8-10, 1960, Natl. Acad. Sci.- 
Natl. Research Council, Wash., D.C., 69-118 (1961 

Interpretation of the appearance potentials of secondary ions, 
M. B. Wallenstein and M. Krauss, J. Chem. Phys. 34, 
No. 3, 929-936 (Mar. 1961) 

Observations on the chemiluminescent reaction of nitric oxide 
with atomic oxygen, H. P. Broida, H. I. Schiff, and T. M. 
Sugden, Trans. Faraday Soc 57, No. 458, 259-265 (Feb. 
1961). 

Studies of the low-temperature distillation of hydrogen iso- 
topes, T. M. Flynn, Cryogenics 1, No. 2, 1-5 (Dec. 1960); 
Advances Cryogenic Engr. 6, 236-244 (1961). 

Inert enclosed pump for shaped flow of ultraclean solutions, 
M. W. Barnes and R. H. Noyce, Rev. Sci. Instr. 32, No. 
3, 353 (Mar. 1961 

On some partial differential equations of Brownian motion of 
a free particle, A. Ghaffari, (Abstract) Proc. Intern. Conf. 
on Partial Differential Equations and Continuum Mechan- 
ies (Univ. of Wisconsin, Madison, Wis., June 7-15, 1960), 
348-350 (1961). 

tadiation protection standards, L. 
824-831 (May 1960). 

Atomic beam frequency standards, R. C. 
Beehler, and C. 8. Snider, IRE Trans. 
I-9, No. 2, 120-132 (Sept. 1960). 

VLF phase perturbation associated with meteor shower ion- 
ization, C. J. Chilton, J. Geophys. Research 66, No. 2, 
379-383 (Feb. 1961 

Departures from the Saha equation for ionized helium. | 
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N. Thomas and J. B 
588-595 (Mar. 1961 
Tests for regression coefficients when errors are 
M. M. Siddiqui, Annals Math. Stat. 31, 

(Dec. 1960). 

Electron diffraction studies on solid a-nitrogen, FE. M. 
and L. Marton, Acta Cryst. 14, Pt. I, 11-19 
1961). 

A comment on the NRL solar Lyman-Alpha results, J. T. 
Jefferies and R. N. Thomas, Astrophys. J. 133, No. 2, 
606-607 (Mar. 1961 

Hard gallium alloys for use as low contact resistance electrodes 
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